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ABSTRACT 

We have used HST imaging of the central regions of the globular cluster 47 Tucanae (= NGC 104), taken 
with the WFPC2 and ACS cameras between 1995 and 2002, to derive proper motions and U- and V-band 
magnitudes for 14,366 stars within 100" (about 5 core radii) of the cluster center. This represents the largest set 
of member velocities collected for any globular cluster. The stars involved range in brightness from just fainter 
than the horizontal branch of the cluster, to more than 2.5 mag below the main-sequence turn-off. In the course 
of obtaining these kinematical data, we also use a recent set of ACS images to define a list of astrometrically 
calibrated positions (and F475W magnitudes) for nearly 130,000 stars in a larger, ~ 3' x 3' central area. We 
describe our data-reduction procedures in some detail and provide the full position, photometry, and velocity 
data in the form of downloadable electronic tables. We have used the star counts to obtain a new estimate for the 
position of the cluster center and to define the density profile of main-sequence turn-off and giant-branch stars 
into essentially zero radius, thus constraining the global spatial structure of the cluster better than before. A 
single-mass, isotropic King-model fit to it is then used as a rough point of reference against which to compare 
the gross characteristics of our proper-motion data. We search in particular for any evidence of very fast- 
moving stars, in significantly greater numbers than expected for the extreme tails of the velocity distribution in 
a sample of our size. We find that likely fewer than 0. 1 %, and no more than about 0.3%, of stars with measured 
proper motions have total speeds above the nominal central escape velocity of the cluster. At lower speeds, the 
proper-motion velocity distribution very closely matches that of a regular King model (which is itself nearly 
Gaussian given the high stellar density) at all observed radii. Considerations of only the velocity dispersion 
then lead to a number of results. (1) Blue stragglers in the core of 47 Tuc have a velocity dispersion smaller 
than that of the cluster giants by a factor of \/2, consistent with the former being on average twice as massive 
as normal, main-sequence turn-off stars. (2) The velocity distribution in the inner five core radii of the cluster is 
essentially isotropic, and the detailed dependence of cr^ on R for the brighter stars suggests that heavy remnants 
contribute only a fraction of a percent to the total cluster mass. Both of these results are in keeping with earlier, 
more realistic multimass and anisotropic models of 47 Tuc. (3) Using a sample of 419 line-of-sight velocities 
measured for bright giants within R < 105", we obtain a kinematic distance to the cluster: D = 4.0 ± 0.35 kpc, 
formally some 10%-20% lower than recent estimates based on standard CMD fitting, and more consistent with 
the value implied by fitting to the white-dwarf cooling sequence. And (4) by fitting simple models of isotropic, 
single-mass stellar clusters with central point masses to our observed cr^iR) profile, we infer al-a upper limit 
of M, < 1OOO-15OOM0 for any intermediate-mass black hole in 47 Tuc. The formal best-fit hole mass ranges 
from if only the kinematics of stars near the main-sequence turn-off mass are modeled, to ~7OO-8OOM0 if 
fainter, less massive stars are also used. We can neither confirm nor refute the hypothesis that 47 Tuc might lie 
on an extension of the M, - a relation observed for galaxy bulges. 

Note: all online material is also available at http://www.astro.le.ac.uk/^dml31/47tuc.html| 

Subject headings: globular clusters: individual (NGC 104) — astrometry — stellar dynamics 
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1. INTRODUCTION 

Galactic globular clusters, which are ancient building 
blocks of the halo, represent an interesting family of "hot" 
stellar systems in which some fundamental dynamical pro- 
cesses have taken place on time scales comparable to the age 
of the universe. Intermediate in mass between galaxies and 
open clusters, globulars are unique laboratories for learning 
about two-body relaxation, mass segregation and equiparti- 
tion of energy, stellar collisions and mergers, and core col- 
lapse. 

The whole concept of core collapse, linked to the gravother- 
mal instability which may develop due to the negative specific 
heat of self-gravitating systems, was first investigated theoret- 
ically in the 1960s and observed indirectly in the 1980s (see 
jMevlan & Heagie 1997 for a review). The stellar density in 
the core may increase by up to six orders of magnitudes dur- 
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ing the phases of deep collapse. This significantly increases 
the frequency of interactions and collisions between stars. 

Binary stars play an essential role during these late phases 
of the dynamical evolution of a globular cluster, as they trans- 
fer energy to passing stars and can thus strongly influence 
the cluster evolution — enough to delay, halt, and even re- 
verse core collapse. At the same time, stellar collisions are 
effective in destroying binaries (the outcome of most binary- 
binary interactions being the destruction of one participant), 
in hardening those that remain, and in ejecting stars towards 
the outer parts of the cluster Observational evidence of pos- 
sible products of stellar encounters include blue stragglers. 
X-ray sources, pulsars, and high-velocity stars. 

In their pio neering radial- velocit y study of the globular M3 
= NGC 5272. lGunn & Griffinl (l!979) noted the puzzling pres- 
ence of two stars that they called "interlopers". These are 
two high-velocity stars located in the core of the cluster, both 
about 20" from the centre. They have radial velocities rela- 
tive to the cluster mean of H-17.0 km s"' and -22.9 km s"\ 
corresponding to 3.5 and 4.5 times the velocity dispersion in 
the core. These radial velocities are still close enough to the 
mean radial velocity of the cluster ((vr) ~ -147 kms~', high 
enough to make contamination by field stars very unlikely) to 
carry a stron g implication of membership . 

Similarly, iMevlan. Dubath. & Mavor (1991) discovered 
two high-velocity stars in the core of the globular cluster 47 
Tucanae. Located respectively at about 3" and 38" from the 
centre, they have line-of-sight velocities relative to the cluster 
of -36. 7 km s"' and +32.4 km s"', corresponding to > 3 times 
the core velocity dispersion but appearing in a total sample 
of only 50 radial velocities. Repeated observations over 1.5 
years indicated that neither of these two stars is a binary or a 
pulsating star, and high-resolution echelle spectra confirmed 
their luminosity classes and, consequently, their membership 
in 47 Tucanae. 

Prompted in 1995 by the presence of these four unusually 
fast-moving stars in two different globulars, and their poten- 
tial link to the extreme dynamical processes in high-density 
environments, we decided to investigate further the core of 47 
Tucanae, the closest of the two clusters. The capability limit 
of the radial velocity observations which could be obtained 
from the ground in the crowded core of 47 Tuc having already 
been reached, we concluded that if more progress was to be 
made in the search for high-velocity stars, it had to be made 
by obtaining proper motions — a task for which only HST is 
suitable. 

We thus used WFPC2 to obtain images of the core of 
47 Tucanae, at three different epochs over four years be- 
tween 1995-1999, in order to perform precise astrometry and 
obtain a complete census of high-velocity stars. Choosing 
the F300W (w t/-band) filter allowed stars to be measured 
over the whole color-magnitude diagram, from the red-giant 
branch to well down the main sequence, ultimately yielding a 
velocity database of unprecedented size for a globular cluster 
Meanwhile, subsequent observations of the center of 47 Tuc, 
by unrelated WFPC2 and ACS imaging programs between 
1999-2002, have provided extremely useful supplements to 
our original dataset. 

In this paper, we present our analysis of these HST data. 
Not surprisingly, we have found it possible to address a num- 
ber of issues beyond simply characterizing the stellar veloc- 
ity distribution (e.g., see iMinniti et an il997>. But the latter 
does remain our primary focus here, and our look at other 



questions (estimating the distance to 47 Tuc; defining the in- 
ternal velocity-dispersion profile as a function of stellar lu- 
minosity/mass; assessing the possibility of a compact cen- 
tral mass concentration) is not as comprehensive. However, 
we are also providing full details of the data themselves, in- 
cluding extensive tables of star-by-star astrometry, photome- 
try, and proper-motion solutions. This thorough census of the 
stellar distribution and kinematics in 47 Tuc, used together 
with sophisticated modeling techniques, will ultimately allow 
for unique and precise constraints to be placed on relaxation 
processes, stellar collision and ejection rates, and many other 
aspects of the dynamical structure and evolution of globular 
cluster cores. 

It is perhaps worth noting that previous studies of inter- 
nal globular-cluster dynamics using HST-based proper mo- 
tions (Drukier et al. 2003; McNamara. Harrison. & Anderso^ 
2003) have employed samples of ~1000 member stars and 
tended to focus on deriving the stellar velocity disper- 
sion very near the cluster centers. The largest sample of 
ground-based proper motions comes from the analysis by 
^n Leeuwen et al. (2000) of 9847 stars in NGC 5139 = 
uj Centauri, using obse rvations over a ~50-year baseline. 
Ivan de Ven et alJ (l200(j) have used a high-quality subset of 
2295 of these stars to explore the internal dynamics of this 
large cluster and estimate a distance to it. Our full velocity 
sample for 47 Tuc includes 14,366 stars, and the majority of 
these prove useful for a variety of precise kinematics analyses. 

1.1. Outline of the Paper 

We begin in ^2.1l bv giving the basic details of the WFPC2 
and ACS image sets that we ha ve us ed to derive proper mo- 
tions for stars in 47 Tuc. Section l2T2l then focuses on the con- 
struction and astrometric calibration of a comprehensive cata- 
logue of positions and F475W magnitudes for nearly 130,000 
stars in one central ACS field measuring about 3' on a side. 
This "master" star list is presented in Table|3 and both it and 
an associated image that we have made of the field are avail- 
able from the online edition of the Astrophysical Journal. We 
then use thi s list to re-evaluate the coordinates of the center of 
47 Tuc. In ^2.3l we describe our procedures for performing lo- 
cal coordinate transformations of the data at other epochs into 
the master reference frame, in order to obtain relative proper 
motions for as many stars as possible in a rather smaller area 
{R < 100") of the sky. 

Section |3l discusses our derivation of the velocities them- 
selves, focusing on statistical properties such as goodness-of- 
fit and error distributions in order to identify a useful working 
sample for kinematics analyses. A catalogue of U, V, and 
F475W photometry, epoch-by-epoch displacements, and as- 
sociated proper motions for 14,366 stars is given in Table |5l 
This is also available electronically, along with an SM code 
which extrac ts an d plots the data for any given star's position 
vs. time. In ^3.21 we also describe a set of line-of-sight veloc- 
ities that we ultimately use to estimate a kinematic distance to 
47 Tuc. 

In ^4. II we use our master star list to construct the number- 
density profile at R < 100" for stars brighter than the main- 
sequence turn-off. This provides a direct extension of a wider- 
field, ground-based V-band surface-brightness profile already 
in the literature. Combining these data, we fit a standard 
single-mass and isotropic King ( 1966) model to the cluster, 
to give a rough framework for the ph ysica l interpretation of 
some of our results. In particular, in ^4.21 (supplemented by 
Appendix 15} we describe the calculation of projected, two- 
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TABLE 1 

Basic Data on 47 Tucanae = NGC 104 



Property 



Reference 



Cluster Center (J20()0) 
Galactic Coordinates 
Apparent Magnitude 
Integrated Colors 
Main-Sequence Turn-off 
Metallicity 

Central Surface Brightness 

King (1966) Core Radius (V < Vjo) 

King (1966) concentration 

Foreground Reddening 

Field Contamination (V < 21) 

Heliocentric Distance: 

Main-Sequence Fitting 
Main-Sequence Fitting 
White Dwaif 
Kinematic 



a = 0i>24'"05?67, 5 = -72°04'52" 
£ = 305?9, fo = -44?9 
Vi„i = 3.95 

(B-y) = 0.88. (f/-V)=1.25 
Vto = 17.65 
[Fe/H] = -0.76 

Hvo = 14.26 ±0.26 mag arcsec"^ 
ro = 20"84±5"05 
c = log(r,/ro) = 2.01 ±0.12 (Wo 
E{B-V) = QM 

Sfoic = (0.8 ± 0.2) stars arcmin"^ 



Z) = 4.85±0.18kpc 
D = 4.45 ±0.15 kpc 
Z) = 4.15±0.27kpc 
D = 4.02 ± 0.35 kpc 



62 



8.6 ±0.4) 



Central Velocity' Dispersion fm* ~ 0.85 Mq): 

Line-of-sight aSR = 0) = 

Plane-of-Sky a^(R = 0)-- 



11.6±0.8kms-' 

: 0.609 ±0.010 mas yr"' 



this paper, i|2.2.2l 
HaiTis (1996) 
HaiTis (1996) 
Harris (1996) 

Zoccah et al. (200n : Percival et al. (2002) 
HaiTis (1996) 
this paper, ! 4. 1 
this paper, ! ~ 
this paper, i|4. 1 



Harris (1996) 

Ratnatunga & Bahcall (1985) 



Gratton et al. (2003^ 
Percival et al. (2002) 
ZoccaUetal (2001') 
this paper. 



this paper, il6.5l 
this paper jl6.5l 



dimensional proper-motion velocity distributions for generic 
iKing models. 

Sections|5]and|6lthen examine various aspects of the stellar 
kinematics in the central 5 core radii of 47 Tuc. Some prelim- 
inary results from earlier stages of this work have been pre- 
sen ted in conference proceed ings by King & And erson ( 200 11) 
and lMcLaughlin et alJ J2003I) . which naturally are superseded 
here. 

In § ^5.11 and 15.21 we construct the one- and two- 
dimensional distributions of proper motion and compare them 
both to Gaussians and to King models with finite escape ve- 
locities. We look especially for evidence of stars with total 
speeds on the plane of the sky exceeding the nominal central 
escape velocity of 4 7 Tu c, but find only a few dozen potential 
candidates. Section l531 summarizes the overall properties of 
these high -velocity stars. 

In ^6. II we go on to compare the velocity dispersion of blue 
stragglers in our field to that of similarly bright stars on the 
cluster's giant branch. Section considers the run of ve- 
locity dispersion with clustercentric radius, as a function of 
stellar magnitude, and obtains an estimate of t he av erage ve- 
locity anisotropy in the central regions. Sec tion l63l then com- 
pares the velocity dispersion profile of the brighter stars in 
our proper-motion sample to that of our much smaller radial- 
velocity sampl e, to derive a kinematic estimate of the distance 
to 47 Tuc. In ^6.5l we focus on the kinematics at the smallest 
projected radii, to fit the proper-motion velocity dispersions 
there with models based on those o f .Kingi til 966) but allowing 
for the possible presence of a dark central point mass. 

We should emphasize from the start that, although 
47 Tuc is kno wn to be rotating ( Mevlan & Mavor 1986; 
I Anderson & Kin g 2003a), we do not attempt to include this 
complication in any of our kinematics analyses. The justifi- 
cation for this is essentially that we are working here only on 
relatively small scales, in regions of the cluster for which ro- 
tation is indeed dynamically dominated to a large extent by 
random stellar motions. This point has been made previously 
bv iMevlan & Mavor Cl 986.) . and we illustrate it again, quan- 



titatively, in ^6.4l of this paper. 

For reference throughout the paper, some basic data on 47 
Tuc are provided in Table Some of the numbers there rely 
on new analyses of the HST data that we have collected. Note 
up front the small field contamination predicted by the Galaxy 
model of Ratnatunga & Bahcal] (1985), which for the area 
covered by our proper- motion sample (very roughly, « 3.5- 
4 arcmin^) amounts to of order 3 (±2) interloping field stars 
brighter than V < 21. For the majority of our work, this is 
evidently a negligible effect. 

2. HST ASTROMETRY AND PHOTOMETRY 
2.1. The Available Data and General Approach 

This project began with a series of WFPC2 expo- 
sures of the center of 47Tuc in 1995, 1997, and 1999 
(GO-5912, GO-6467, and GO-7503, PI Meylan). The 
goal of these observations was to search for the proper- 
motion analogues of the high-velo city "cannonball" stars that 
iMevlan. Dubath. & Mayoi]lll99 1') had found with CORAVEL 
radial velocities from the ground (and which have well- 
known cou nterparts in the Galactic globular cluster M3; 
iGunn & Griffia.1979.) . These original exposures are confined 
within the inner ~4-5 core radii {R < 100") of the cluster. 
They were taken with the F300W U -band) filter in order 
to suppress the background from the red giants and thus allow 
better position measurements of the more numerous stars at 
the main-sequence turnoff and fainter. As a result, we found 
it possible to measure accurate motions not only for the fast- 
moving stars in the cluster, but for many thousands of the av- 
erage members as well. Furthermore, subsequent (unrelated) 
HST observations of the core of 47 Tuc have nearly doubled 
our original four-year time baseline and allowed the deriva- 
tion of more precise proper motions. 

Additional WFPC2 images of the same central field as 
the Meylan pointings were obtained in 1999 and 2001 (GO- 
8267 and GO-9266, PI Gilliland) through the similarly short- 
wavelength filter, F336W. More recently, ACS images that 
cover the same area were obtained in 2002 for various call- 
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TABLE 2 

WFPC2 AND ACS Observations of 47 Tucanae 



Program 



Data set 


ID 


Kb. 


Filter 


Date 




MEYLANel 


5912 


15 


F300W 


25 Oct 1995 = 


1995.82 


MEYLANe2 


6467 


16 


F300W 


03 Nov 1997 = 


1997.84 


GILLILUl 


8267 


28 


F336W 


05 Jul 1999 = 


1999.51 


MEYLANe3 


7503 


16 


F300W 


28 Oct 1999 = 


1999.82 


GILLILU2 


9266 


11 


F336W 


13 Jul 2001 = 


2001.53 


WFC-MEUR 


9028 


20 


F475W 


05 Apr 2002 = 


2002.26 


HRC-MEUR 


9028 


40 


F475W 


05 Apr 2002 = 


2002.26 


HRC-BOHL 


9019 


10 


F475W 


13 Apr 2002 = 


2002.28 


WFC-KING 


9443 


6 


F475W 


07 Jul 2002 = 


2002.52 


HRC-KING 


9443 


20 


F475W 


24 Jul 2002 = 


2002.56 



bration programs (Pis Meurer, King, and Bohlin), all through 
the slightly redder F475W filter We have reduced all of these 
images from the HST archive and included them in our analy- 
sis. The ACS data in particular have proven extremely useful 
in providing much higher-precision position measurements 
than are possible with the lower-resolution WF chips of the 
WFPC2. In addition, the even coverage of these ACS-WFC 
data (as contrasted with the lop-sided WFPC2 footprint) al- 
lows us to construct an accurate, uniform, and nearly com- 
plete census of stars — independently of any proper-motion 
goals — within about ±1.5 arcminutes of the cluster center. 

All these data sets and their basic attributes are listed in Ta- 
ble|2l There are 182 independent exposures taken as parts of 
ten distinct sets, which we refer to loosely as ten "epochs" 
spanning nearly seven years in total. Our general approach 
to collating these for analysis is first to combine all the ex- 
posures for each epoch so that we have a single position and 
flux for each star measured in a frame natural to that epoch. 
This intra-epoch averaging also gives an empirical estimate of 
the error in position and flux for each star at each epoch. We 
then compare the positions of stars measured at the different 
epochs to derive proper motions. 

The proper motions we measure are, of course, simply 
changes in the relative positions of stars over time. The many 
observations are taken at different times, at different point- 
ings and orientations, through different filters, and with dif- 
ferent instruments. Each observation therefore has a different 
(and a priori unknown) mapping of the chip coordinates to the 
sky. Before we can compare relative positions of stars mea- 
sured in different images, we must transform all our positions 
into a common reference frame. We have chosen to use the 
WFC images of the GO-9028 data set to construct this "master 
frame," since these images have a large and very even spatial 
coverage. The large majority of stars found in any of our data 
sets w ill b e found in the GO-9028 data set. 

In ^2.21 then, we define this master frame and discuss its 
astrometric and photometric calibration. We also use it to find 
a new estimate for the cluster center, which will be useful for 
our later analyses. After this, we briefly describe the process 
by which we transform multi-epoch observations into the ref- 
erence frame ( ^2.3> . In ^we detail our derivation of the 
proper motions themselves and define a working sample for 
investigation of the cluster kinematics in the rest of the paper 

2.2. The Master Star List 

Special care is required in constructing the star list for the 
reference frame, since even if a star is not optimally measured 



in the GO-9028 data set, we still want to allow for the possi- 
bility that it might be found well in other epochs. Our primary 
goal is therefore to make the master list as complete as possi- 
ble. 

The GO-9028 data set has some small dithers and some 
large dithers about a central pointing. We do not want the 
master star list to be affected by the location of the inter- 
chip gap in the centeral pointing, so we made use only of the 
central-pointing image and the pointings that had ditherings 
larger than the inter-chip gap to generate an unbiased master 
list. This amounts to 13 images out of the 20 in GO-9028. We 
first fit a PSF to every peak in each of these 13 images, and 
then corrected each peak's raw pixel position (from the _f It 
images) for distortion according to the prescription in Ander- 
son (2006, in preparation). We next determined the transfor- 
mation from each image into the frame of the central pointing 
and identified a star wherever a peak at the same master-frame 
location was found in 7 or more of the 13 individual images. 

This gave us a list of 129,733 coincident peaks, covering 
an area of about 202" x 202", with positions in the distortion- 
corrected frame of the central-pointing image (j8cd01a9q). 
Plotting the star list on a stacked image that we made of the 
field shows that no obvious stars are missing from the list. 
There are, however, a (relatively small) number of PSF arti- 
facts that occurred in the same place in all images and which 
were therefore misidentified as stars. 

2.2.1. Purging of Non-stellar Artifacts 

In order to isolate PSF artifacts in the master list, we went 
through the list star by star For each star we found every 
fainter neighboring "star" within 25 pixels. In Figure [J we 
show the distribution of neighbors for all saturated stars only 
(about 3000 altogether). On the horizontal axis we plot the 
distance in WFC pixels from the bright star to the fainter 
neighbor, and on the vertical axis we show the magnitude dif- 
ference. PSF artifacts occupy a clearly recognizable region in 
this parameter space (e.g., the clump of points at a ^10-pixel 
distance and '^8-magnitude flux difference). 

We therefore define the discriminating line shown in Figure 
^to distinguish real stars from (possible) artifacts. For every 
star in turn in the full master list, we have flagged all neigh- 
bors that fall above this line. These correspond to objects that 
are too close to brighter stars to be considered trustworthy, and 
they should not be used in any detailed analyses. This proce- 
dure is bound to reject some real stars along with true artifacts, 
but it does so in a quantifiable (and therefore correctable) way. 
The end result is a more robust and better-defined star list, 
consisting ultimatel y of 1 14,973 reliable stars. The full star 
list is presented in ^2.2.6l below. after we have described the 
astrometric and photometric calibration of the data. 

2.2.2. Finding the Cluster Center 

The ACSAVFC images provide us the widest, the most uni- 
form, and the deepest survey to date of the central regions of 
47 Tuc. Previous images taken by WFPC2 have a very asym- 
metric and non-uniform coverage, due to the lop-sided shape 
of the WFPC2 footprint and the gaps between the chips. By 
contrast, the dithered set of WFC images have essentially uni- 
form coverage out to a radius of about 100", which is nearly 
five core radii. These data therefore permit us the best deter- 
mination to date of the cluster center Indeed, given the high 
degree of completeness in the present star counts, it is difficult 
to imagine any significant improvement in the near future. 
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Fig. 1. — Distribution of neiglibors found around saturated stars in tlie 
GO-9028 F475W exposures used to construct the master frame. Eacli point 
represents a fainter neiglibor of a saturated star and is plotted at the appropri- 
ate distance and magnitude offset from the main star. PSF artifacts show up 
as very faint objects relatively close to the star. All points above this line are 
considered to be potential artifacts. Some real stars are necessarily rejected 
as well, but in a very regular way so that completeness can be defined well. 

Note that a determination of the center from star counts 
should not depend on incompleteness corrections, provided of 
course that incompleteness is a function only of stellar mag- 
nitude and clustercentric radius. Nor should mass segregation 
enter the problem, so long as the stellar distribution is radially 
symmetric. Thus, in determining the center we work directly 
with our master star list, making no attempt to correct for ei- 
ther of these effects. 

We begin with our list of 1 14,973 stars after culling faint 
neighbors (possible PSF artifacts). For each star we have a 
position {x,y) in the reference-frame coordinate system, ro- 
tated in order to align the y axis roughly with North. Next, an 
array of trial centers (xcjo) is defined, and for each trial cen- 
ter in turn we find all the stars within a radius of 1500 WFC 
pixels (75") from that center We then divide these stars into 
sixteen, 22.5°-wide pie wedges — shown schematically in the 
left panel of Figure|2] — and use two tests to compare the distri- 
butions of stars in the eight distinct pairs of opposing wedges. 

In the first test, we look at the difference in the total number 
of stars between the members of the opposing-wedge pairs de- 
fined by each ixQ,yo) in our pre-defined array of trial centers, 
and we form a "quality-of-fit" statistic as the sum of the dif- 
ferences over the eight wedge pairs. The coordinates which 
minimize this statistic then define a "median" estimate of the 
cluster center. The middle panel of Figure |5] shows a contour 
plot of the sum of differences for a 3"x3" subgrid of trial 
centers, with (0, 0) corresponding to the best center ultimately 
implied by this method (the conversion to calibrated right as- 
cension and declination is discussed below). 

The second test is to generate the cumulative radial distribu- 
tion for stars within each of the sixteen wedges for any speci- 
fied {xo,yo). We then find the absolute value of the integrated 
difference between the radial distributions in any two oppos- 
ing wedges, and define a quality-of-fit statistic as the sum of 



these absolute differences over the eight wedge pairs. The 
right panel of Figure |2l shows a contour plot of this statistic 
over a grid of trial centers, with (0,0) again corresponding to 
the best-estimate coordinates which minimize the statistic. 

These two approaches give quite consistent positions for 
the cluster center. To get an idea of the accuracy of our center, 
we again take eight pairs of opposing wedges. For each pair, 
we find the location along the wedge axis that minimizes the 
difference between the cumulative radial distributions for the 
two wedges independently of any others. This yields eight 
estimates of the center along different axes using indepen- 
dent samples of stars. From the scatter among these inde- 
pendent estimates, we find that our final center is good (in the 
master-frame coordinate system) to about ±5 WFC pixels (or 
~ 0."25) in both directions. This uncertainty is indicated by 
the errorbars at (0,0) in the middle and right-hand panels of 
Figure |2] 

2.2.3. Astrometric Calibration 

We now have a position for the cluster center in the refer- 
ence frame, which is based on the distortion-corrected and 
rotated frame of the first image of GO-9028. In order to 
transform the Master-frame positions into absolute RA and 
Dec, we used the image header information from several 
WFPC2 images (u2ty0201t, u2vo0101t, u4f 40101r, 
and u5jml2 0dr) to obtain absolute positions for seven 
stars — five stars at the center and two stars in the outskirts. 
These four images were taken at different pointings and ori- 
entations, so they should all use different guide stars and give 
independent estimates of the absolute coordinates. 

Table |3l gives details of these seven stars. First are their 
IDs in our master star list and calibrated F475W magnitudes 
(both of which items are described in general below). Follow- 
ing this are the stars' locations in the master frame, both in 
terms of pixel positions and in terms of relati ve RA and Dec 
offsets from the cluster center determined in ^2.2.21 Then we 
list the average absolute RA and Dec (J2000) obtained from 
the header information in the four WFPC2 images. Combin- 
ing the absolute positions of the five central stars with their 
relative offsets in the reference frame then sets the absolute 
astrometric zeropoint of our master-frame system. The two 
outer stars Q and 77 are used to fix the orientation angle. As is 
also stated in the bottom line of Table |3l the absolute position 
of the cluster center is 

RA(J2000) = 00^24™05?67 ±0?07 
Dec(J2000) = -72°04'52."62±0f'26 , ^ ' 

where the uncertainties come from the averaging of the five 
stars and essentially reflect the ±0"25 internal uncertainty in 
the uncalibrated master-frame coordinates of the cluster cen- 
ter. This absolute calibration should be good to about 0. 1 arc- 
second throughout our ^^3' x 3' master field, although we note 
that it may ultimately suffer from a small (~1") inaccuracy 
if the positional errors of the HST guide stars for the WFPC2 
frames we have used are correlated (see Taff et al. 1990). 

The position of our adopted cent er is interm ediate to those 
determined by Guhathakurta et all (1992) and ICalzetti et alJ 
11993). This is illustrated in Figure|3l the left panel of which 
shows a 20" x 20" region about our adopted center with the 
5 centra l reference stars in Table |3l marked. Our star e is the 
star that Guhath akurta et al . used as a reference position (their 
star E). Our absolute coordinate for this star differs from theirs 
by about 1.5 arcseconds. The center position as estimated by 
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QUALITY OF FIT FOR CENTER QUALITY OF FIT FOR CENTER 

FROM MEDIAN ANALYSIS FROM CUMULATIVE DISTRIBUTIONS 
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Fig. 2. — Left: a shematic illustrating the wedges used to home-in on a cluster center; see text. Middle: contours of equal quality-of-fit statistic for determining 
the center using a "median" analysis, i.e., by comparing the total numbers of stars in pairs of opposing wedges. Right: contours of equal quality-of-fit statistic for 
determining the center using the cumulative radial distributions of stars in opposing wedges. In the middle and right panels, dotted contours represent changes of 
20% in the quantity plotted. EiTorbars represent the estimated uncertainty in the center. 



TABLE 3 

Stars Used for Astrometric Calibration, and the Resulting Cluster Center 



Point/Star 


Master Frame ID 


F475W 




)'raw 


ARA 


ADec 


RA (J2000) 


Dec (J2000) 








[pixe 


s] 


[arcsec] 


[hh:mm:ss] 


[dd:mm:ss] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


a 


M052296 


13.6 


1991 


415 


-4.532 


-5.983 


00:24:04.69±0.10 


-72:04:58.58±0.12 


P 


M060833 


14.3 


1961 


239 


-4.924 


+2.935 


00:24:04.61±0.10 


-72:04:49.69±0.11 


7 


M061604 


13.8 


2163 


182 


+5.180 


+3.720 


00:24:06.79±0.10 


-72:04:48.90±0.11 


5 


M056400 


13.6 


2208 


281 


+6.809 


-1.665 


00:24:07. 15±0. 11 


-72:04:54.25±0.11 


e 


M056630 


12.9 


2013 


319 


-2.903 


-1.413 


00:24:05.04±0.11 


-72:04:54.08±0.11 


c 


M004529 


14.2 


2113 


812 


-6.731 


-74.669 


00:24:04.21 ±0.11 


-72:06:07.29±0.I0 


»? 


M034527 


13.4 


3380 


545 


+62.662 


-26.247 


00:24: 19.25±0. 10 


-72:05:18.77±0.12 


Center 






2066 


277 


0.000 


0.000 


00:24:05.67±0.07 


-72:04:52.62±0.26 



Note. — Key to columns: 

Column (1) — Label for the star in Fig.|3] See text for the distinction between stars a, /3, 7, 5, and e vs. stars C, and 77. 
Column (2) — Stellar ID on the sequential numbering system for the full master-frame star list of Tablel^ 
Column (3)— Calibrated F475W magnitude. 

Column (4) — "Raw" .v position, in pixels, in the top chip (WFCl) of the j 8cd0 la9q_f It frame from program GO-902S. 
Column (5) — "Raw"^' position, in pixels, in the top chip (WFCl) of the j8cd01a9q_f It frame from program GO-9028. 
Column (6) — RA offset in arcsec (positive Eastward) from the cluster center in the master-frame system. 
Column (7) — Dec offset in arcsec (positive Northward) from the cluster center in the master-frame system. 
Column (8) — Calibrated, absolute light ascension. 
Column (9) — Calibrated, absolute declination. 



o 

QJ 
Q 




R. A. 




1 -1 

ARA (arcsec) 



Fig . 3 . — Left: stars in a 20" x 20" field about our cluster center. Point size is correlated with F475W magnitude. The 5 stars used in the absolute calibration 
are indicated (see Table [sl. as are the cluster centers estimated bv Guhathakurta et al^ lI992) and .Calzetti et aL lI993i) . Right: a close up of the inner A" x4". 
The uncertainty in the center is about =bO"25 in each coordinate. 
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iGuhathakurta et aP is labeled with a "G," placed at the posi- 
tion they report relative to star e (not at the absolute coordi- 
nates given in their paper). We also mark the center found by 
[Calzetti et al. with a "C." Finally, the right panel of Figure 13 
shows a close-up of the 4" x4" box around our center. 

2.2.4. F475W Photometry 

In the course of fitting PSFs to find positions, we also com- 
puted the average F475W fluxes of stars in the master frame. 
These fluxes include a spatially-dependent correction, of or- 
der ~ 3%, for the fact that the ACS flat fields were designed to 
preserve surface-brightness rather than flux (Anderson 2006, 
ACS ISR in preparation). However, they were calculated from 
only the inner 5x5 pixels around each stellar peak, and they 
correspond to a 60-second exposure. To calibrate the photom- 
etry, we must turn this flux into that which would be measured 
through the standard "infinite" aperture of 5 arcseconds (100 
WFC pixels) in 1 second. We first compute the ratio of our 
measured flux to the flux contained within 10 pixels, or 0."5. 
This ratio is 1.238. The encircled energy curves in the ACS 
Handbook ("Pav lovskv et al.ll2005h then show that 92.5% of 
the light should be contained within this radius. Thus, we 
scaled all our fluxes by a net factor of 1.339 upward and then 
divided by 60 (seconds) befo re adding th e VEGAMAG ze- 
ropoint of 26.168 (De Marchi et al.1 12004 to obtain a final, 
calibrated F475W magnitude for each star in the master list. 

2.2.5. Completeness Fractions 

The broad, uniform coverage of our WFC master frame also 
makes this a useful data set for constructing an accurate sur- 
face density profile for the cluster We can detect almost all 
the stars there are, so it should be possible to come up with 
the definitive radial profile from the center out to nearly 5 core 
radii. There are, however, two issues that complicate the con- 
struction of any density profile from star counts in globular 
clusters: incompleteness and mass segregation. Mass seg- 
regation really just means that the density profile can differ, 
at least in principle, for stars in different mass (magnitude) 
ranges. It is a physical effect, separate from any instrumen- 
tation or data-reduction issues, and we discuss it briefly in ^ 
below, where we actually derive a density profile for the inner- 
most parts of the cluster Incompleteness, on the other hand, 
is a technical limitation of the observations themselves. 

Incompleteness is always a joint function of both the im- 
ages and the algorithm used to find stars in the images. It 
can arise from several sources. A particular star might not be 
found because (1) it is too close to a brighter neighbor and is 
not bright enough to generate its own peak in the image; (2) it 
could land on a defect in the chip or it could be hit by a cosmic 
ray; or (3) it could be close to the background and not bright 
enough to generate a peak above the noise. The fact that we 
have generated our list from many observations all at different 
pointings saves us from (2). And since the stars we analyze 
in this paper are always several magnitudes brighter than the 
faintest stars that can be detected in the master frame, issue 
(3) is not important for us. The first source of incompleteness 
(bright-star crowding) is the only thing we need to concern 
ourselves with here. 

The usual strategy for evaluating incompleteness involves 
applying a set data-reduction procedure to a large number 
of artificially generated data sets. This is a rather daunting 
prospect in our case, and we have instead addressed the prob- 
lem directly from the algorithm we used in ^2.2.1l to identify 



(likely) non-stellar artifacts in our original list of 129,733 PSF 
peaks in the master frame. 

First, we used the ACS images from program GO-9028 to 
create a circular image of the region R < 150", with a uni- 
form pixel size of 0."05 pixel"'. (Note that this reference im- 
age completely contains the master field itself, which is non- 
circular and measures only ~ 3' on a side. This is because 
we required a point to be covered by at least 7 ACS pointings 
to contribute to the master star list, but only one was suffi- 
cient to build the circular mosaic.) We then went pixel by 
pixel through the intersection of this circular frame with our 
master field, made a list of all stars found by our PSF-fitting 
within a 25-pixel (1 ."25) radius from each point, and recorded 
the faintest magnitude, m^ax, of the stars which could have 
survived the artifact-purging procedure of ^2.2.11 and Figure 
n This yields an estimate of the limiting magnitude at every 
pixel in our master frame. 

Given the limiting mmax at every point in our field, we cal- 
culated a "completeness fraction" for every star in the master 
list individually. Knowing the position (x^,,y^,) and magnitude 
of each star in the list, we looked at the values of m^ax at 
all points within 100 pixels (= 5") of (;c, ,y*) and counted the 
number of pixels for which m^ax < m^,. That is, we computed 
the fraction / S [0, 1] of the local area around each star where 
an identical star could fall but be discounted as unreliable by 
our criterion in ^2.2. II (or not be detected at all). The local 
completeness fraction at (x*,y*,OT*) is then just c = (1 -/). 

With c defined pointwise in this way, every star in our mas- 
ter list is interpreted as representing 1 /c actual stars, and cor- 
rected density profiles follow simply from the sums of 1 /c 
over all points within specified areas on the sky. Again, we 
actually construct such a profile in ^ 

2.2.6. Image of the Central Regions and the Final Star List 

As was mentioned just above, as part of our estimation of 
completeness fractions we constructed a meta-image of the 
innermost R < 150" from the center of 47 Tuc. This con- 
venient reference image, 47TucMaster . fits, oriented in 




Fig. 4. — ACS-based image of the central R < 150" in 47 Tuc, inside of 
which we identify stars for the unbiased master hst of Table |4| The image 
is available in .fits format from the online edition of the Astrophysical 
Journal. Pixel size in the online image is 0"05 per pixel, and coordinates 
are defined with the cluster center at {x,y) = (3001 , 3001) px. North is up and 
East is to the left. 
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the usual North-up, East-left way, is available through the on- 
line edition of the Astrophysical Journal. A low-resolution 
version of it is shown in Figure |3 

Table |3 finally presents all important information on the 
129,733 coincident PSF peaks in our master field, which is 
fully contained in the circular area of the meta-image. (A sam- 
ple of Table H can be found at the end of this preprint.) Only 
1 14,973 of these peaks can be said confidently to be bona fide 
stars, but to be comprehensive, we have here retained (and 
flagged) th e 14,7 60 peaks which could be PSF artifacts ac- 
cording to ^2.2.11 The table gives the offset of each detection 
from the cluster center in arcseconds; the calibrated F475W 
magnitude; the absolute RA and Dec in real and celestial for- 
mats; the faintest magnitude a star at each position could have 
and still be found; a flag indicating whether the star survives 
the artifact purging; the local completeness fraction for the 
brightness and position of each star; a serial ID number; and 
{x^y) coordinates in both our meta-image and the original im- 
age of the ACS/GO-9028 centi-al pointing. 

The absolute astrometric calibration of Table |4] should be 
good to about 0"1, but the relative positions should be much 
better than this — particularly for bright, unsaturated stars and 
those with small separations, for which we estimate an accu- 
racy of '^O^'OOl. Of course, all the positions refer specifically 
to the epoch (2002.26) of the GO-9028 data set. 

2.3. Reducing Images from Multiple Epochs 

With a well-defined reference frame in hand, the next step 
toward deriving proper motions is to measure the positions 
of stars in individual images taken at different epochs and 
transform these into the master coordinate system. We be- 
gan this task by measuring each star in each data set in Table 
121 with appropriate PSF, derived according to the method 
of [Anderson & King (2000). For the ACS-WFC observa- 
tions, we used a single PSF to treat the entire 2-chip associ- 
ation. We also used a single PSF for the entire ACS-HRC 
chip. The PSF does vary significantly with position over 
the ACS, but our data were reduced before the methods of 
[Anderson & Kine (2006) were developed to deal with this. 
Nevertheless, Anderson (2002) shows that the biggest effect 
of this variation on astrometry is a small bias of 0.01 pixel 
in the positions. Since all our images are well-dithered, this 
error averages out and is included in the internal uncertain- 
ties for the positions. The constant-PSF assumption can also 
introduce systematic errors of up to 0.03 magnitudes in the 
photometry, but this is not of concern to us here. 

All the raw measured posit ions were corrected for d istor- 
tion using the prescriptions in lAnderson & Kin^ i2003 b) for 
WFPC2 data and Anderson (2006, in preparation) for ACS 
data. These corrections include the global and fine-scale dis- 
tortion corrections, as well as a correction for the 68th-row 
defect. Charge-transfer inefficiency should not be an issue, as 
the background is relatively high in the images that we have 
used. 

It was then necessary to combine all the observations of 
each star in the multiple pointings at each epoch. In all cases 
we adopted the centermost pointing as the "fiducial" frame for 
each epoch, and used general 6-parameter linear transforma- 
tions (along the lines of those described in eq. O below) to 
transform positions from the individual, distortion-corrected 
pointings for that epoch into the central frame. In finding 
these transformations for WFPC2 images, we treated the PC 
chip and the three WF chips independently of the others; for 
the ACS images, the HRC chip and two WF chips were like- 
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Fig. 5. — Illustration of the criteria defining stars, on the various chips 
of different detectors, with positions that are "well-measured" enough to be 
used in deriving proper motions. Examples are shown at one epoch for each 
of the PC and WF chips on the WFPC2, and the WFC chips and HRC chip 
of the ACS. The discriminating fines drawn are applied to the list of positions 
and fluxes of all stars found at every epoch on each chip, with any detections 
falling above the appropriate line discarded as too uncertain to contribute 
meaningfully to the detemiination of a proper motion. The horizontal axes 
are given here in instrumental magnitudes to indicate the relationship be- 
tween positional error and raw S/ N. Th e location of the subgiant branch in 47 
Tuc (cafibrated V ~ 17.2; cf. Fig. 1101 is clear from the sharp increase in the 
density of points in each panel. 



wise considered independently. In this way, we found an av- 
erage position for each star in each of 27 chip-epoch com- 
binations. The standard deviation of a star's position in the 
independent pointings at each epoch defines the uncertainty 
in X and y locations. 

At this point we examined the errors in the stars' positions 
as a function of their flux at each epoch in order to define a 
criterion to define which stars were well measured and which 
were not well measured within each data set. Figure|5lshows 
a plot of the position errors as a function of the instrumental 
magnitude -2.51og[o[A^DA'] for some sample chip-epoch com- 
binations. We drew the fiducial lines shown in Figure [S] to 
distinguish well measured from poorly measured stars; only 
stars falling below these lines in any of the data sets were 
retained for any further analysis. In addition to this, we dis- 
carded all stars found within the pyramid-affected region of 
any of the WFPC2 chips (ie., any at chip coordinates x < 100 
or y < 100). We finally rejected any saturated bright stars 
from the star list for every data set. 

This left us with an average position, including uncertain- 
ties, and a flux for every star that we deemed well measured 
in each chip of each data set. The positions were still in 
the distortion-corrected coordinate systems of the individual 
epochs, however, and needed next to be transformed into the 
system defined by our master frame. 

We should note explicitly here that the GO-9028 WFC data 
(epoch 2002.26) were subjected to the above analyses (and to 
the coordinate transformations described next) in exactly the 
same way as were any of th e oth er epochs in Table |2l — even 
though the master frame in ^2.2l was itself defined using data 
from this epoch. This is because in constructing the master 
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frame we used only 13 images from the 2002.26 WFC set of 
20 and made no effort to impose any "quality control" on the 
stellar positions or magnitudes. However, any given star in 
the master list could, in principle, be observed in one or more 
of the seven other pointings at this epoch; or alternately its po- 
sition might be so poorly constrained (relative to the criterion 
illustrated in Figure |5} that it is not useful for proper-motion 
measurements. Thus, to make optimal use of the GO-9028 
WFC data requires that they be re-analyzed in parallel with 
the other epochs, with the master frame viewed simply as an 
externally imposed construct. 

2.3.1. Linking with the Master System 

The transformation into our master coordinate system is 
non-trivial. Given the small proper-motion dispersion in 47 
Tuc (see Tabled and our short (< 7 year) time baseline, we 
need to resolve displacements of order 0."001 and less for the 
stars in our sample. The pointing of HST is good to only 
0."1 at best, however, and thus we cannot transform our po- 
sitions to an absolute frame with anywhere near the accuracy 
required. 

The only information we have about how the coordinate 
system of one frame is related to that of any other comes from 
the positions of stars that are common to both frames. Each 
observation typically has thousands of stars in common with 
the master frame. Once these stars are matched up, the two 
sets of positions can be used to define a linear transformation 
from one coordinate system to the other. Specifically, if we 
have a list of positions {xi^yc / = 1 , . . . , A/"} for A/" stars at some 
epoch, and a list {X, , Yj] for the same stars in the master frame, 
then we specify 

Xi = A-{xi-XQ) + B-(yi-yQ) + XQ 
Yi = C-(xi-XQ) + D-(yi-yo) + Yo, 

for constants A, B, C, D, and (x(),yo), (Xq, Fo) to be determined. 
It may look like there are eight free parameters in this trans- 
formation, but in fact we have the freedom to pick the zero- 
point offset in one system. In particular, if we choose (jccyo) 
to be the centroid of the Af matched stars in the non-master 
frame, then (Xo,Yo) is necessarily the centroid of the group 
in the Master system. With these set, then, standard least- 
squares regression can be applied to solve for the linear terms 
A, B, C, and D from the 2J\f observed pairs of coordinates. 
Our last concern is the question of how many stars should be 
used to do this. 

2.3.2. Local Transformations 

If we were to use all the stars in common between the mas- 
ter system and the data from any other epoch, the solution 
of equation (|2j would correspond to the best global transfor- 
mation between the two coordinate systems. Unfortunately, 
such a chip-wide solution may not give the best transforma- 
tion for a given point. In particular, although all of our frames 
have been distortion-corrected, any residual distortions will 
introduce systematic errors in a global transformation. Such 
distortion errors tend to be cumulative, i.e., they are larger 
for stars that are farther apart. As a result, the distance be- 
tween two nearby stars can be measured much more accu- 
rately than the distance between stars at different corners of a 
chip. We therefore decided to perform more local transforma- 
tions, using smaller groups of relatively nearby stars to define 
different transformations at different positions in an observed 
frame. Such a scheme may sound computer intensive, but it 



is straightforward enough to implement and is not, in fact, ex- 
ceedingly slow. Given that we will never be able to remove 
distortion perfectly, this is a useful way to minimize its effect 
on our results. 

Two sources of error influence the choice of the number 
A/lrans of stars to usc in solving the system of equations (|2j for 
local transformations. First, the positions and {X,-,y,} 

themselves are obviously subject to some uncertainty. If A 
denotes a representative value of this uncertainty, then the 
average transformation is fundamentally uncertain at a level 
~ A/VA/trans, which suggests that we would like A/trans to be 
as large as possible. Second, the uncorrected residual distor- 
tion introduces appreciable systematic error if the "nearby" 
stars used to define the local transformation at any point come 
from too far away. This implies that we would like A/Irans to 
be as small as possible. 

After some experimentation, we found that a reasonable 
compromise between these opposing tendencies was reached 
with A/trans = 45 for our data. In practice, we took a single 
star from the list of well-measured positions at one epoch; 
found this target star's nearest 55 neighbors; and matched 
these neighbors to positions in the master frame. We then 
used these 55 pairs of coordinates {not including the target 
star itself) to solve for the coefficients in equation (|2ji; dis- 
carded the 10 stars which deviated most from the solution; 
and re-solved for the transformation using the 45 remaining 
stars. This was then used to transform the original, target star 
only into the reference frame. These steps were repeated for 
every well-measured star in the combined frame for each of 
our ten epochs, until we had measurements of RA and Dec 
position (relative to the cluster center) vs. time for all stars in 
a single, unified coordinate system. 

We also performed this procedure using A/Irans = 20 and 
A/Irans =100 and did not find significant differences, in gen- 
eral, from our results with A/Irans = 45. However, the lower 
number is approaching the limit of what we feel comfortable 
with in terms of vulnerability to small-number statistics, while 
the higher number is coming close to bringing in too-distant 
neighbors that could well be affected by uncorrected distor- 
tion errors. 

There is another, somewhat more subtle source of error in 
our local transformations, which comes from the fact that all 
the stars are physically moving with respect to each other 
from one epoch to the next. Thus, even if a star's position 
could be measured perfectly in every image, it is impossible 
to associate its location {x,y) at some epoch perfectly with its 
position iX^Y) in the master frame using a finite network of 
A/Jrans moving neighbors. Put another way, any one network of 
neighbor stars might in fact have a real net motion relative to 
the cluster center, but our approach assumes that all such net- 
work motions are identically zero. This error gets averaged 
away over many networks, of course, so that the mean veloc- 
ity we estimate for any group of stars will be unbiased. How- 
ever, any estimate of the velocity dispersion, a, is affected. 
Details of how we correct for this are given in Appendix IaI 
which also discusses the correction of velocity dispersion for 
unequal measurement errors. Here we simply note that the 
net effect of the local-transformation error is an artificial in- 
flation of the intrinsic stellar cr by a factor of (1 + 1 /A/ii-ans)'"'^, 
or about 1 .01 1 for our chosen A/Irans = 45. 

2.3.3. V- and U -Band Photometry 

In addition to calibrated F475W photometry for all stars in 
our master field ( ^2.2.4t . we have obtained calibrated V-band 
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and F300W (roughly U -band) magnitudes for the subset of 
stars falling in the WFPC2 field of the original Meylan expo- 
sures (GO-5912, GO-6467, GO-7503; see Table |3 and note 
that V exposures were also taken as part of these programs for 
photometric purposes only). As we will describe below, all 
stars for which we have derived proper motions were required 
to be detected in at least one of these three early epochs, and 
thus the full proper-motion sample has multicolor photome- 
try. Although the standard V bandpass is not substantially 
different from F475W, it is useful for relating our analyses to 
various ground-based observations, and to such things as the- 
oretical stellar mass-luminosity relations. The shorter U -band 
photometry is useful for constructing a color-magnitude dia- 
gram of our velocity sample in order to investigate how CMD 
position might influence stellar kinematics. 

V photometry of the WFPC2 fields was calibrated 
against an HRC image of the core of 47 Tu c, us- 
ing the VEGAMAG zeropoints in iDe Marchi et alJ j2004 . 
Star-by-star comparison with th e calibrated photometry in 
|Bica. Ortolani. & Barbuv (1994) shows agreement at the ~ 
0.05-mag level on average. The U magnitudes were calibrated 
using the VEGAMA G zeropoint in the WFPC2 Data Hand- 
book dBa^get et alJ l2002i) . The main-sequence turn -off in 
our V, ({/ - V) CMDs agrees well with that in Edmond s et alJ 
(|2003). We rep ort th e V and U photometry for our proper- 
motion stars in ^3.11 where we now define the velocity sam- 
ple. 

3. VELOCITY SAMPLES 

Our multi-epoch astrometry data are most certainly not ho- 
mogeneous. The 2002 epoch is relatively uniform, thanks to 
the even ACSAVFC coverage; but the earlier WFPC2 obser- 
vations come from a different camera with different pointings, 
dithering patterns, fields of view, and resolutions. Thus, not 
every star in any one of the ten data sets of Table |2] can be 
foun d in a ll of the other nine. 

In ^3.11 we describe how we go from a heterogeneous sam- 
ple of multi-epoch positions to a homogeneous sample of 
proper motions. We also examine the quality of the fits of 
straight lines to the position-vs.-tim e dat a, and the uncertain- 
ties in the final proper motions. In ^3.21 we discuss the third 
component of velocity, that along the line of sight. Specifi- 
cally, we have a large sample of ground-based radial veloci- 
ties for stars in 47 Tuc, and we are interested in using these in 
conjunction with our proper-motion sample t o ob tain a kine- 
matic estimate of the distance to the cluster ( ^6.3> . 



3.1. Proper Motions 

To define a working sample of plane-of-sky velocities for 
kinematic analyses, we have chosen to work only with stars 
which are observed in at least three separate data sets prior to 
year 2000 (i.e., in three or four of the MEYLANel, MEY- 
LANe2, MEYLANe3, and GILLILUI epochs in Table |2}, 
and in at least one of the ACS data sets (any of the WFC or 
HRC fields) from 2002. Thus, we only derive proper motions 
for stars that have a minimum of four separate {x,y,t) mea- 
surements spanning a minimum of about 4.4 years (2002.26- 
1997.84 = 4.42), and all derived velocities are effectively "tied 
down" by a precise ACS position measurement. In practice, 
most of our stars are in fact measured 6 or more times over 
the full 6.7-year timespan (1995.82-2002.56) of our epochs, 
and a good number do in fact appear in all 10 of the data sets 
listed in Table |2| 
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Fig. 6. — Outline of the fields of view of our master star list (large square 
area; Table derived from ACS program GO-9028, and of the original 
WFPC2 observations of GO-5912, 6467, and 7503. Only stars in the inter- 
section of the two fields are retained for inclusion in the full proper-motion 
sample of Tablelsl 



We further decided to include stars in the proper- motio n 
sample only if they also appear in the master star list of ^2.21 — 
even if another ACS "epoch" might contain the star Ulti- 
mately, this has some effect on the spatial distribution of the 
velocity sample. Figure|6lshows the RA and Dec offsets from 
our estimated cluster center for a subset of stars in our mas- 
ter list (the larger, rotated square field), and the positions of 
stars in the original MEYLAN epochs (the four smaller, more 
densely filled squares). Evidently, the combination of all our 
requirements gives us proper motions over a patch on the sky 
with roughly the familiar WFPC2 pattern, but with slightly 
larger gaps than normal between the chips (cf. 32. 3> and with 
the outermost corners shaved off. 

We consider a star to be "observed" at any pre-2002 
(WFPC2/MEYLAN or GILLIL) epoch only if it is detected 
in at least ten of the individual pointings/ditherings that went 
into that observation, and if it survives the culling based on 
position error vs. instrumental flux discussed around Figure|5] 
above. For any of the higher-precision, 2002 ACS observa- 
tions, we take a star to be detected only if it is identified in at 
least four individual pointings/ditherings and if it again sur- 
vives the appropriate culling by position error. These criteria 
are imposed to ensure that, whenever we include the mea- 
surement of a star's position at any one epoch in our proper- 
motion determination, we also have enough information to 
estimate accurately the position uncertainty at that epoch. 

The estimation of position uncertainties is of critical im- 
portance to the derivation of proper-motion uncertainties, and 
thus to the ultimate inference of intrinsi c, err or-corrected stel- 
lar kinematics. As was suggested in 32.31 we take a star's 
position in any epoch to be the mean of the positions mea- 
sured in the separate ditherings that were combined for that 
epoch and transformed into the master reference frame. The 
uncertainty is then simply the standard sample deviation of 
those independent multiple measurements. Note that we do 
not work with the uncertainty in the mean position (which 
would involve dividing the standard deviation by \/Nd for Nd 
ditherings), because the theory of Unear regression — which 



HST Proper Motions in the Core of 47 Tuc 



11 



we use to derive velocities — actually requires that the square 
of the errorbars on the data be unbiased estimates of the vari- 
ance in the position measurements. 

Thus, given (x,y) = (aaooo, (^2000) positions and uncertainties 
as functions of time for any star that satisfies the minimum 
criteria just set out, we derive a plane-of-sky velocity in each 
direction by a sta ndard, error-weigh ted least-squares fit of a 
straight line (e.g.. iPress et alJll992l Section 15.2), allowing 
both the slope and the intercept to vary. This procedure auto- 
matically yields uncertainties in each component of velocity. 
The measurement-error distribution in each component of ve- 
locity for any single star will be Gaussian, with a mean of 
zero and a dispersion equal to the fitted least-squares error- 
bar, if the position measurement errors are similarly Gaussian 
distributed. We have assumed that this is the case. 

The most important advantage of performing a weighted 
least-squares regression here is that it returns a value for 
each straight-line fit. Knowing the number of degrees of free- 
dom in the fit (= Ne-2, where 4 < A^,, < 10 is the number 
of independent epochs at which the star's position was mea- 
sured), it is then possible to calculate the probability that the 
of the fit could have occurred by chance if the tru e mo- 
tion of the star were really linear (see lPress et al.llT992 ). We 
can then make use of this to exercise some quality control 
over the proper motions, by excluding from kinematics anal- 
yses any stars with fitted velocities whose probabilities are 
lower than some specified threshold. 

By fitting straight lines to our position-vs.-time data we 
are obviously assuming that the stars are moving at effec- 
tively constant velocity. To justify this in general, note that 
the characteristic gravitational acceleration in the core of 47 
Tuc is of order (see Table 1, and §^ and |6l below) oq ^ 
o-Q/ro ~ (0.6mas yr"')^/20arcsec ~ 2 x 10~^mas yr"^. Even 
integrated over 7 years, the velocity change induced by such 
an acceleration is only ~ 10"'* mas yr"', which is, as we shall 
see, a small fraction of the velocity uncertainties we infer Of 
course, this does not exclude the possibility that some stars 
could be significantly accelerated by "nonthermal" processes 
such as stellar or black-hole encounters, or by virtue of being 
in tight, face-on binaries, or in some other way. Such stars will 
simply not be described well by straight-line motion, and the 
probabilities inferred from our weighted linear regression 
will reflect this fact. Thus, a low P(x^) can reflect legitimately 
nonlinear data as well as simple "bad" measurements. 

In all, we have 14,366 stars with RA and Dec positions 
measured to satisfactory precision in at least three pre-2000 
epochs and at least one 2002 ACS epo ch, an d with measured 
V and U magnitudes (as described in ^2.3.3> . Table|5] which 
is published in its entirety in the electronic version of the As- 
trophysical Journal, contains the instantaneous J2000 RA and 
Dec positions in our master frame (epoch 2002.26) for all of 
these stars, expressed both in arcseconds relative to the cluster 
center determined in ^2.2.2l (see also Table 1) and in absolute 
celestial coordinates. The F475W, V , and U magnitudes of 
each star are also reported. The RA and Dec offsets (also in 
arcsec) from the nominal master-frame position, and their un- 
certainties, are then given for every epoch in which the star 
was detected. The proper-motion velocities and uncertainties 
implied by the weighted straight-line fitting to the offsets vs. 
time then follow, along with the values for the fits and the 
probabilities P(x^) that these values could occur by chance 
if the motion is truly linear. We have not tabulated the inter- 
cepts of the linear fits, as these only reflect the choice of an 
arbitrary zeropoint in time and are of no physical interest in 



the constant-velocity case. (A sample of Table|5lcan be found 
at the end of this preprint.) 

Table |5] includes all stars for which we have estimated ve- 
locities, regardless of whether or not the x^ values of the fits 
have "good" probabilities. We stress again that a reliable sam- 
ple for kinematics work should exclude stars with very low 
P(X^), but that some such stars could still be of interest for in- 
vestigations of non-constant-velocity phenomena (which we 
do not pursue in this paper). 

Two points should be noted regarding the connection be- 
tween Table |5] and the larger master list of stars in Table |3 
First, at the 2002.26 WFC-MEUR epoch in Table El the RA 
and Dec offsets from the absolute master-frame positions are 
consistent with within the uncertainties for most stars, but 
they do not vanish exactly — even though this is the epoch that 
was used to construct t he m aster frame. The reason for this is 
that the master list in ^2.2l was defined using a very specific 
subset of the 20 individual exposures comprising th e WF C- 
MEUR data set, while the relative positions found in ^2.3l and 
used here to compute proper motions were allowed to come 
from different combinations of the 20 pointings. As a result, 
the transformations in the latter case cannot be expected in 
general to give positions identical to those in the master list, 
and the offsets for this epoch in Table|5]essentially reflect sta- 
tistical noise. Second, there are some stars in Table [S] for 
which no offset at all is given at the WFC-MEUR epoch — 
an entry of "n/a" appears instead — even though all stars in the 
proper-motion table are guaranteed by construction to appear 
in the master list. In these cases, the uncertainties in the stars' 
positions at the master-frame epoch were larger than accept- 
able according to our flux-based criterion defined in Figure|5] 
above, and thus these data were excluded from our fitting for 
proper motions. 

Having derived proper motions for the best-measured stars 
in our field, we now describe the selection of a velocity sub- 
sample which is best suited for kinematics analyses. In par- 
ticular, we can be more specific about how the x" probabil- 
ities in Table [S] are used to cull stars with particularly poor 
proper-motion measurements (or nonlinear motions), and we 
can quantify the distribution of the velocity uncertainties. 

3.1.1. "Good" and "Bad" Proper Motions 

Figure0shows plots of position (in milliarcseconds of RA 
and Dec offset from the master-frame position) vs. time (in 
years from the arbitrary zeropoint 1999.20) for two stars with 
relatively good proper-motion determinations. In the left pan- 
els is a bright star, near the cluster center, which has position 
measurements in Table|5]for 8 of our 10 data sets/epochs. In 
the right panels is a much fainter star, farther from the clus- 
ter center, which has fewer {x,y,t) datapoints. In both cases, 
the fitted velocities jia and /ia in the RA and Dec directions 
are given in units of mas/yr, with jia defined to be positive 
for Eastward motion. The x^ probabilities for the linear fits 
of the two proper-motion components are also listed, showing 
that these data are fully consistent, within measurement error, 
with the basic assumption of constant velocity. 

The much smaller position uncertainties in the most recent, 
ACS data relative to the pre-2002 (WFPC2) epochs are note- 
worthy. The right-hand panels particularly illustrate how the 
ACS epochs play a crucial role in defining the overall motions 
of our stars throughout the ~7-year baseline of the observa- 
tions. Indeed, if only the 4-year span of WFPC2 data had been 
fit to find a proper motion for the faint star in FigureQ the RA 
component would have had the opposite sign (though with a 
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Fig. 7. — Two examples of stars with good qualities of fit for straight- 
line proper motions: P(x^) ^ 0.001 in both velocity components. The ID 
numbers of the stars correspond to those in Tables |4] and |5] Plots such as 
these can be generated for any star in Tablelslusing an SM macro available in 
the electronic edition of the Astrophysical Journal. 

larger errorbar) from that obtained when the ACS data are 
included. Evidently, relying on only a few epochs of astrom- 
etry, even with the HST, can lead to spurious proper motions 
for some individual stars. 

Figure |8] next shows two stars with much less satisfactory 
proper- motion fits. Both stars here are fairly bright and near 
the cluster center, and both show more scatter about the best- 
fit linear velocity than is acceptable. On the left-hand side, the 
problem is primarily with the RA displacements in the upper 
panel, where the early (WFPC2) epochs do not match well 
onto the precise later (ACS) position measurements and Pix^) 
is uncomfortably small. On the right-hand side, the scatter in 
both the RA and Dec motions is so large (relative especially 
to the small errorbars on the ACS positions at year ~ 2002) 
that there is effectively no confidence that the best-fit constant 
velocity is an accurate representation of the data. 

Plots like those in Figures Q and |8] can be generated for 
any star from the data in Table |5j using an SM macro 
(pmdat . mon) that we have packaged and made available in 
the online edition of the Astrophysical Journal. After inspect- 
ing many such graphs and comparing the kinematics of sam- 
ples of stars defined by imposing various lower limits on the 
allowed value of P(x^) for the fitted velocities, we eventually 
decided to include in detailed analyses only stars which have 

Pix^) > 0.001 (3) 

for both RA and Dec components of proper motion. This cri- 
terion defines a velocity sample which cleanly excludes "bad" 
data such as those illustrated in Fig.|8l(as well as — again — 
potentially perfectly good data that are just nonlinear) and 
shows quite robust kinematics, in the sense that samples de- 
fined by revising the P(x^) threshold in equation (|3j moder- 
ately upward — even by as much as an order of magnitude — do 
not have significantly different statistical properties. 

3.1.2. Proper-Motion Uncertainties 
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Fig. 8. — Two examples of stars with poor qualities of fit for straight-line 
proper motions. With P(x") < 0.001 in one or both velocity components, 
stars such as these are not included in any working samples for kinematical 
analysis in this paper. They are, however, retained in Table Isl Plots such as 
these can be generated for any star in Tablelslusing an SM macro available in 
the electronic edition of the A.strophysical Journal. 

Figure |5] shows the distribution of uncertainties (least- 
squares errorbars) in the RA component of proper motion 
for stars in Table |5l which are brighter than V < 20 and have 
P(X^) > 0.001 for the fitted values of both and ^s- We 
have divided this subsample into three broad bins of projected 
clustercentric radius, R. It is immediately apparent that stars 
at /? > 20" have systematically higher velocity errorbars, typ- 
ically by factors of ^2, than those at < 20". This is because 
this radius is completely contained in the high-resolution, PC 
chip of the WFPC2 camera (see Figure |6|l, which affords 
higher precision in the measurement of (pre-2002) stellar po- 
sitions than the WF chips which cover the larger radii in our 
field. From this alone it is clear that our velocity uncertainties 
correlate with clustercentric position. 

Panels on the left-hand side of Fig.|5]plot the proper-motion 
uncertainty against the velocity itself, fi^. The colored 
lines, which correspond to different magnitude bins within 
each radial bin, connect the median errorbar in each of 
a number of discrete /i^ bins. The right-hand panels then 
show the corresponding normalized histogram of velocity un- 
certainties in each radius/magnitude bin. These show that, for 
stars of a given magnitude, the distribution of A^ not only 
peaks at smaller values for R < 20", but it is more sharply 
peaked there as well. Put another way, stars at the larger clus- 
tercentric radii have a stronger tail toward high velocity er- 
rorbars. Conversely — and quite naturally — at any fixed clus- 
tercentric radius, fainter stars always have larger average un- 
certainties and broader distributions. Indeed, the stars in the 
faintest magnitude bin illustrated here (19.5 <V < 20) have, 
at radii R > 25", "typical" velocity errorbars of order 0.6 mas 
yr"' or more, which is comparable to the intrinsic velocity 
dispersion at the center of 47 Tuc (Table Thus, stars of 
such faint magnitude will be of limited use for the statistical 
characterization of intrinsic cluster kinematics. Stars that are 
fainter still have even larger velocity uncertainties and are en- 
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Fig. 9. — Characterization of velocity uncertainties Aq in the RA com- 
ponent of proper motion, fia- Left-hand panels are scatter plots of Aa, with 
colored lines connecting median values in bins of fia, for stars of different 
magnitudes. Right-hand panels ai'e normalized histograms of the distri- 
butions, showing the tendency for the uncertainties to increase with V magni- 
tude and to be smallest in general at S < 20", corresponding to the region of 
the proper-motion field covered by the PC chip in the early WFPC2 images. 

tirely useless in this context. Thus, we impose the magnitude 
limit 

y < 20 (4) 

when choosing stars for any kinematics analysis. This is not 
to say that the velocities of individual faint stars are always 
unreliable — we shall see evidence to the contrary in ^5.31 — 
but only that their group properties are poorly constrained. 

The distribution of velocity uncertainties for the Dec com- 
ponent of motion, fig, is essentially identical, in all respects, 
to that shown for the RA component in Figure |9] The basic 
message of these plots is that, due to strong correlations with 
stellar magnitude and position, the errorbars on our fitted ve- 
locities cover a wide range of values and cannot be considered 
even approximately equal. Nor, in general, are they negligible 
relative to the intrinsic stellar motions. It is therefore impor- 
tant to account properly for the effects of measurement error 
on the observed proper-motion distribution and its moments 
(i.e., velocity dispersion). The details of this are discussed in 
Appendix El where we also consider the statistical implica- 
tions of the local-transformation approach to obtaining rela- 
tive proper motions (cf. ^2.3.2> . 

3.1.3. The Sample for Kinematics Analyses 

Figure[lO|shows color-magnitude diagrams and spatial dis- 
tributions of the 14,366 stars with proper motions listed in 
Table |5] split into a sample in the left-hand panels which we 
consider useful for kinematics analyses (i.e., which comprises 
only stars satisfying both criteria in eqs. 13 and 0) and a 
sample in the right-hand panels that we exclude from such 
work (i.e., which consists of stars with Pix^) < 0.001 in either 
component of proper-motion velocity, and/or with magnitude 
V > 20). 

The CMD of the upper left-hand panel exhibits a very well 
defined cluster sequence, including a red giant branch, a main- 
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Fig. 10. — Color-magnitude diagrams and spatial distributions of proper- 
motion stars which are used for subsequent kinematics work (left-hand pan- 
els), and those which are not (right-hand side). See text for details of the 
criteria by which stars are rejected. 



sequence turn-off at about the appropriate Vtq = 17.65 (Table 
[!), and even a small blue-straggler population. This is, in fact, 
not unexpected given that the foreground contamination in our 
small field (with an area of « 3.5-4 arcmin^) is expected to 
amount to perhaps 3 ± 2 brighter than V = 21 (from the fore- 
ground star density estimated by Ratnatunaa & Bahcall 19851 
see see Table Note that we have no stars brighter than 
y = 14 with measured velocities at all. Such stars were gener- 
ally saturated in the 2002 ACS frames, and thus did not meet 
our selection criteria for the proper-motion sample. 

The rejected stars in the right-hand panels of Figure ^| 
amount to only ^ 10% of the full sample in Table|5] 30% of 
them (417/1392) are rejected simply because they are fainter 
than y = 20. The 975 others are rejected solely because they 
have a low probability for the linear fit to one or both com- 
ponents of their velocity. Of these, as Figure^Jsuggests, the 
majority are brighter than y = 1 8 .5 and located within R < 20" 
of the cluster center, where the stellar density is highest and 
crowding is most problematic. 

It is important to check whether the rejection of these 
stars leads to significantly different kinematics for the "good" 
proper-motion sample than would have resulted if all stars 
were retained. Thus, in Figure ^2 we show the cumulative 
distributions of RA and Dec velocities for stars with V < 18.5 
and R < 20, split into samples with P(x^) > 0.001 (accepted 
into a final kinematics sample) and P(x^) < 0.001 (rejected). 
Kolmogorov-Smirnov tests applied to these indicate that the 
velocity distributions of the rejected stars are consistent with 
their having been drawn from the same parent distribution as 
the accepted stars. Thus, excluding the former from analysis 
does not give biased kinematics in the end but simply reduces 
the net uncertainty in our final results. 

Before going on to use our sample of 12,974 good proper 
motions to investigate these kinematics in detail, we first 
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Fig. 11. — Comparison of the cumulative RA and Dec velocity distri- 
butions for bright stars in the inner 20" of our proper motion field, for the 
sample of stars with P(x^) > 0.001 in both proper-motion components, vs. 
those stars with P(x^) < 0.001 in at least one component. In each panel, Pks 
is the probability (from a Kolmogorov-Smirnov test) that the proper-motion 
distributions of the "good" and "bad" stars are in fact drawn from the same 
parent distribution. The high values of this statistic show that our results are 
not biased by excluding the poorest velocity fits from kinematical analysis. 

briefly describe so me c omplementary line-of-sight velocity 
data that we use in ^6.3l (onlv) to derive a kinematic distance 
to 47 Tuc. 

3.2. Radial Velocities 
3.2.1. The Data 

The radial velocities used in this study come from two 
different instruments: the Geneva Observatory's photo- 
electric spectrometer CORAVEL ( Baranne, Mavor, & Poncet 
Il979t) mounted on the 1.5 m Danish telescope at Cerro 
La Silla, Chile, and the Rutg ers Fabry -Perot spectrom- 
eter (e.g., Gebhardt etai] 11994 at CTIO Details of flie 
CORAVEL data are given bv 'Mav or et all ( fT983l fT984 >) and 
[Mevlan, Dubath, & Mavor (1991); typically, the measure- 
ments for relatively bright giants and subgiants have uncer- 
tainties of ~ 0.6 km s"'. The Fabry-Perot velocities include 
some previously published (Gebhardt et al. 1995), and some 
newly measured (Gebhardt et al. 2006, in preparation) from 
two runs at CTIO in 1995: one using the 1.5-meter telescope 
on 16-22 June 1995, and one with the 4-meter on 6-7 July 
1995. The velocities normally have precisions of 1 km s"' or 
better, depending on the stellar magnitude and crowding. 

All together, there are neai-ly 5,600 stars with CORAVEL 
and/or Fabry-Perot radial velocities, with Fabry-Perot con- 
tributing a large majority of the data. However, each of 
these samples covers a much larger area of 47 Tuc than does 
our proper-motion sample. We work here only with a much 



smaller subset of stars which lie within a radial distance of 
105" of the cluster center and thus overlap the proper- motion 
field (cf. Figure I10> . In addition, we found it necessary to 
exclude from our analysis any Fabry-Perot stars fainter than 
the horizontal branch in 47 Tuc (which is also roughly the 
limiting magnitude of the CORAVEL data). This drastically 
reduces the sample size and thus requires some explanation. 

The basic strategy of Fabry-Perot observations is described 
in Gebhai-dt et al. ( 1994, 1995, 1997). The idea is to scan an 
etalon across a strong absorption line (Ha in our case) to build 
a small spectrum. The instrumental resolution of the Rutgers 
system is around 5000, and the FWHM of the Ha line is larger 
than this. For the 4-meter observations of 47 Tuc, 21 steps 
across the Ha line were used, while the 1.5-meter run ob- 
served two fields using 17 and 41 steps. The full range in all 
three cases was about 5A. The exact wavelength coverage of 
any particular star depends on its location in the field, due to a 
wavelength gradient introduced by the Fabry-Perot. Given the 
full set of scans, DAOPHOT ( Stetson 1987) and ALLFRAME 
JStetso n' 1998) were used to determine the brightness of each 
star in the field. The scans were then combined to make a 
small spectrum, which was fitted for the Ha velocity centroid. 

There are particular advantages and disadvantages of this 
technique compared to the traditional slit echelle spec- 
troscopy employed by CORAVEL. The latter provides large 
wavelength coverage, but it is obviously necessary to choose 
specific stars on which to place the slits. With this configura- 
tion, the observer has no control over what light goes down the 
slit, and neighboring stars can be be a severe contaminant in 
crowded regions. By contrast, Fabry-Perot spectroscopy pro- 
vides only a small wavelength coverage but once the scans are 
completed, every star in the field produces a spectrum. This 
provides an enormous advantage in terms of the number of 
velocities that can be obtained. In principle, difficulties due to 
crowding should also be reduced since the use of DAOPHOT 
allows one to apportion the light appropriately in crowded re- 
gions, i.e., accurate photometry can be obtained for faint stars 
near bright ones. 

Crowding inevitably remains an issue, however, with its 
main effect being a progressively stronger bias in the mea- 
sured velocity dispersion at fainter magnitudes. If the spec- 
trum of a faint star is affected by either a bright neighbor or 
the collective light of the bulk of unresolved cluster stars, the 
measured velocity is pulled closer to the cluster mean, and the 
velocity dispersion th at fo llows is spuriously low. In the con- 
text of our work in ^6.31 where we compare proper-motion 
to radial-velocity dispersions to estimate the distance to 47 
Tuc, this effect would ultimately lead to a short value. There 
are two options for avoiding this. One is to remove the effect 
statistically, by using detailed simulations to define empiri- 
cal, magnitude-dependent corrections for the stellar veloci- 
ties. This approach is pursued by Gebhardt et al. (2006, in 
preparation). The other option is simply to make a magnitude 
cut in the velocity sample, keeping only those stars which are 
bright enough not to be significantly contaminated in the first 
place. This is what we have done here. 

Figure shows the distribution of the combined 
CORAVEL and Fabry-Perot radial velocities for stars in three 
bins of clustercentric radius. The left-hand panels include 
only stars with 11 <V < 14, a magnitude range that has 
been well explored in previous radial-velocity studies, and 
which we have confirmed does not suffer significantly from 
the contamination problems just discussed. The limit V > 1 1 
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Fig. 12. — Histograms of CORAVEL and Fabry-Perot radial veloci- 
ties in the inner 105" of 47 Tuc, comparable to the area covered by our 
proper-motion survey. Left-hand panels are for stars blighter than the clus- 
ter horizontal branch, and right-hand panels are for stars just 0.5 mag fainter. 
The sharp drop in the (error-corrected) dispersion a- from the bright to the 
faint magnitudes — especially at the smaller clustercentric radii — signals a 
crowding-induced bias in the Fabry-Perot velocities. 

excludes only a handful of measured stars in this region of 
the cluster but is imposed to guard against rad ial-velocity "jit- 
ter" in the brightest cluster giants ( Gunn & Griffin 1979) . The 
faint limit V < 14 corresponds to the magnitude of the hori- 
zontal branch. The right-hand panels then include stars just 
slightly fainter than this, in the range 14 < V < 14.5. 

We have calculated the error-corrected velocity dispersion 
(7, (see eq. IA7I ) for the stars in each of the magnitude ranges 
and radial bins defined in Figure [21 These are listed in the 
appropriate panels of the figure, where we have also drawn 
(strictly for illustrative purposes) Gaussians with dispersions 
(7; on top of the velocity histograms. Comparing these re- 
sults for the two magnitude bins at any given radius immedi- 
ately shows the downward biasing of ct, caused by crowding. 
Comparing the different radii confirms this interpretation of 
the situation, since ( 1 ) the difference between cr, for the bright 
and faint samples decreases toward larger clustercentric radii, 
where the stellar densities are lower; and (2) the velocity dis- 
persion at the smallest radii R < 25" for the fainter stars es- 
pecially is dramatically (and unphysically) lower than that at 
larger R. Although cr. for the stars with V < 14 is formally 
lower at R < 25" than in the range 25" < R < 50", by some 
0.7 km s~', the uncertainty in the calculated dispersion is also 
0.7-0.8 km s"' in both radial ranges. Thus, the difference in 
this case is not highly significant and is likely a reflection of 
small-number statistics rather than serious contamination. 

Given the results in Figure we choose to include only 
stars brighter than the horizontal branch, 11 <V < 14, in the 
radial-velocity sample used in this paper There are only 419 
such stars in the inner 105" of 47 Tuc, and thus any kine- 
matic estimate of the distance is fundamentally limited to a 
precision of no better than (2 x 419)"'/- ~ 3.5%. Clearly, it 
is desireable to do better than this, and in principle it is pos- 
sible to enlarge the radial-velocity sample by, say, applying 
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Fig. 13. — Velocity vs. magnitude for stars in our "good" proper-motion 
sample [V < 20 and P(x^) ^ 0.001 for both and fig; from Table Isl. and 
for all stars at R < 105" in the CORAVEL and Fabry-Perot radial-velocity 
sample. Since the biasing problem illustrated in Figure fT2l precludes the use 
of any line-of-sight velocities when V > 14, we do not cun'ently have a single 
star for which all three components of velocity have been reliably measured. 

position-dependent magnitude cuts to remove contaminated 
stars and/or using data from clustercentric radii beyond the 
proper-motion field. However, either of these things would 
require a sophisticated analysis of the Fabry-Perot data and an 
in-depth dynamical modeling that are well beyond the scope 
of this paper. 

3.2.2. Overlap with the Proper-Motion Sample 

Figure shows the RA and Dec components of proper 
motion vs. stellar magnitude for our "good^proper-motion 
sample of ^3.1.31 [i.e., those stars in Table |5] with P(x^) > 
0.001 and V < 20], and also for the combined CORAVEL 
and Fabry-Perot line-of-sight velocity sample in the inner R < 
105" of 47 Tuc. Again, we do not make use here of any data at 
y > 14 in the latter case, but Figure^jemphasizes that radial 
velocities can indeed be measured for many such faint stars. 
The great potential of these data, if their systematic problems 
can be solved, is clear. 

Figure^lalso illustrates the fact that there is currently not 
one star in 47 Tuc which has reliable measurements for all 
t hree c omponents of its velocity. As was also mentioned in 
^3.1.31 stars brighter than V = 14 were not measured in any 
ACS frames, and thus do not appear in our proper-motion 
sample; but fainter than this limit we cannot trust the radial- 
velocity measurements in general. As a result, only statistical 
comparisons of motions on the plane of the sky and along the 
line of sight are feasible at this point. 

Finally, Figure O compares the spatial distribution of the 
419 stars whose radial velocities can be used for a distance es- 
timate (large points), to that of our good proper-motion sam- 
ple (small dots, plotted for only a fraction of the 12,974 stars 
to outline the shape of the field). The circles drawn on this 
graph have radii of 25" (which completely encloses the PC- 
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Fig. 14. — Comparison of the spatial distribution of our good proper- 
motion sample and the useabl e radia l-velocity stars. All positions are plotted 
relative to the cluster center ( i|2.2.2i . and for reference the large, dotted cir- 
cles have been plotted with radii of 25", 50", 75", and 100". 

chip area of the proper -mot ion field), 50", 75", and 100". As 
we discuss further in ^6.31 the obvious differences between 
the non-uniform spatial coverages of the two kinematics sam- 
ples mean that some care must be used in properly comparing 
the measured velocity dispersions to derive a distance to the 
cluster 

4. SPATIAL STRUCTURE AND A MODEL VELOCITY 
DISTRIBUTION 

To optimally explore the implications of the velocity data 
that we have assembled, we would like to relate them to a 
dynamical model for 47 Tuc. Ideally, this should be con- 
strained by the spatial distribution of stars as a function of 
magnitude (mass) and allow for the possibilities of rotation, 
velocity anisotropy, axisymmetry (rather than sphericity), the 
presence of heavy stellar remnants, and so on. The proper 
construction of such a comprehensive, global model is clearly 
a large task that would involve bringing in more data than 
we have collected here (e.g., the spatial and velocity distribu- 
tions of stars at radii beyond the ~ 100" extent of our proper- 
motion sample; information on the stellar luminosity function 
as a function of position in the cluster; etc.). Our much more 
modest goal in this paper is to try and extract as much model- 
independent information as possible from our proper motions, 
and simply to check whether they are broadly consistent with 
a very low-level description of the cluster. 

The zeroth-order model we use as a reference point is the 
familiar single-mass, isotropic, modified isothermal sphere of 
|Kin2 (1966). Note that Mevlan (1988, 1989) has already 
fitted 47 Tuc with more realistic models of multimass and 
anisotropic clusters, but using only radial velocities and the 
surface-brightness profile of subgia nts and g iants as observa- 
tional constraints. Thus, applving iMevlanf s results as they 
are to assess fine details of our proper-motion data would 
not really be any more appropriate than using a single-mass 
and isotropic model (and re-fitting them would, as we've sug- 
gested, involve work beyond the scope of this paper). At any 



rate, iMeylanf s models imply that velocity anisotropy is only 
important at large radii R > 500", i.e., in regions of the cluster 
not probed by our proper motions. 

Similarly, Mevlan & Mavor ( 198^ and lAnderson & Kin3 
(|2()03a) have quantified the rotation of 47 Tuc, using both 
line-of-sight velocities and proper motions. However, in the 
innermost R < 100" covered by our new velocity data, rota- 
tion appears to be present only at the level of < 3 km s~', 
which is only a fraction of the velocity dispersion there (see 
iMevlan & Mavor 1986, and ^6.41 below). Thus, compari- 
son of the proper-motion kinematics with a non-rotating and 
spheri cal m odel can still be profitable in this region. 

In ^4.11 then, we use our master star list from 32.2.61 to 
construct a new density profile for the brightest stars in the 
core of 47 Tuc, and combine this with wider-field, ground- 
based surface photometry to fit a Kine ( 1966) model to the 
cluster's spatial structure. In 34.2l and AppendixlHl we discuss 
the proper-motion velocity distribution predicted by such a 
model. 

4. 1 . Density Profile of Turn-off Stars 

The most important concern for us in constructing a 
surface-density profile from our master star list is the ques- 
tion of mass segregation. This involves massive stars (or 
heavy remnants) sinking toward the cluster center and follow- 
ing a more concentrated radial distribution than lower-mass 
stars, so we cannot really speak of a single density profile 
for the cluster. In fact, the data in Table |3 above could be 
used to explore mass segregation in some detail, but once 
again it is beyo nd the scope this paper. Some early work 
has been done bv'Paresce. de Marchi. & Jedr zeiewskil ( 119951) . 
Anderson & King ( 1996), and Anderson ( 199^! 

The best constraint on the spatial structure of 47 Tuc outside 
the core regions observed by HST is the V-band surface pho- 
tometry collated from a variety of sources by Mevlan ( 1988t). 
The cluster light in this band is dominated by giants and sub- 
giants, all of which have roughly the same (main-sequence 
turn-off) mass. Our immediate interest, then, is in defining an 
improved profile in the central region for stars of this mass 
only. The brightness of the main-sequence turn-off in 47 Tuc 
is V = 17.65 (see Tabled, corresponding roughly to magni- 
tude 17.8 in the F475W bandpass for which our master star 
list has photometry. Thus, after first checking that the ra- 
dial distribution does not vary significantly as a function of 
magnitude for stars brighter than the turn-off, we have used 
Table 13 to construct a single density profile for all stars with 
'MF475W < 17.8. We emphasize that the full master star list has 
to be used for this, and not any sample of only proper-motion 
stars from Table |5j because the latter have a highly irregular 
distribution on the sky and a complicated (and unquantified) 
incompleteness as a function of magnitude and clustercentric 
position. 

We first define a series of concentric circular annuli {Rj}, 
all contained entirely within the area of the master frame. The 
estimates c, of completeness percentage for i ndivid ual stars 
in the master list (Column 10 of Table Bl 



see 32.2.51 1 are then 
used to calculate the density in each annulus: if there are A/la„ 
stars with mF475w < 17.8 in the range Rj <R < Rj+\, then the 
completeness-corrected number is A/'cor = Ym^T ^^^^/"^'^ ^^'^ 
the net density follows as E = Afcor / T^iR^j+i ~^^)- The results 
of this simple exercise are given in Table|6l Note, from com- 
parison of the listed A/'raw and A/^oi, that the corrections for 
incompleteness in this magnitude range are quite modest over 
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Density Profile of Cluster Stars With 
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734.5 ±23.3 



Note. 



Number densities in last column are 

,2 



E(R) = A/'co,/[vr(/?2-^i)] with Ri and R2 the inner 
and outer radii of an annulus (Column 1), 



in units 

of arcmin. Uncertainty in T,{R) is v^/V^/[7r(/?| - 
/?])]. To transform into V-band surface brightness 
on the scale of Meylan (1988) (uncorrected for red- 
dening), set /^v [mag arcsec"^] = (24.237 ± 0.021) - 
2.5 log(E/arcmin"^). 



our entire field. Also, no corrections are made for contami- 
nation by foreground or background stars, since the level of 
such contamination is negligible in this small area of the sky; 
see Table^again. 

It remains only to combine these star counts with the 
wider-field V-band surface photometry taken from iMevlanI 
We first define an effective radius for each an- 
nulus, ^eff; = The number density at each 
Rsff is then converted to a surface brightness as /i(.Reff.>) = 
C-2.5 log I](/?eff./), where C is a constant determined by the 
overlap of our data with the calibrated Mevlan photometry. 
In practice, we only use the ground-based surface bright- 
nesses fly at clustercentric radii R > 20" (to guard conser- 
vatively against any potential seeing-induced blurring of the 
central profile) and obtain C as the median value of [^y + 
2.5 logT^iReff)] for all of our data points with 20" < /?eff < 
70". We find that the new HST number counts are brought 
nicely onto the V-band surface-brightness scale of Meylan 
il988h (uncorrected for reddening) with the transformation 



= (24.237 ±0.021)-2.51og(E/arcmin"^) . (5) 

Figure [Tsl shows the composite surface-brightness profile 
for 47 Tuc. Open squares represent the data from Mevla3 
while the filled circles correspond to our star counts 
after the application of equation (|5J. We have also plot- 
ted, as the open triangles, earlier HST-based number densi- 
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Fig. 15. — V-band surface-brightness profile of 47 Tuc, con- 
structed by combining the ground-based surface photometi'y compiled by 
l^vlan 1 1988) at R > 20" wit h HST star counts at i? < 100" from 
IHowell. Guhathakurta. & Gilhland' f2060) and our own master list (TableHl. 
Only stars brighter than the main-sequence turn-off were included in the HST 
number-density profiles, which were converted to fiy as described in the text. 
The bold, dash-dot curve is the fit of a single-mass, isotropic, Kins 1 1 96fJ 
model, with the parameters indicated (also listed in Table [TT. This model 
only attempts to describe the cluster structure at small radii. The poor fit 
at S > 1000" is due to neglecting velocity anisotropy and multiple stellar 
masses, but the disagreement with fj,v in the extreme halo does not impact 
the core region R < 100" covered by our proper-motion data. 



ties of turn-off stars from 'How ell. Guhathakurta. & Gillilandl 
(H)00, their Figure 4), transformed by = (14.927 ±0.033)- 
2.5 log(I]HGG/arcsec"^) to give the best match to the ground- 
based surface photometry at 20" <R< 70". The agreement 
between the two HST da ta sets is very good overall, con- 
firming the conclusion of Howell. Guhathak urta. & Gillilandl 
( 2000) that 47 Tuc has a smooth, well-defined, and reasonably 
large core that shows no sign of the cluster having undergone 
core collapse. 

The do t-dash curve in FigureFTSlis our fit of a single-mass, 
isotropic. King (1966) model to the combination of our data 
points with the brightness profile of Mevlan ( 1988). Only data 
at R < 1100" were used to constrain the fit since, as we dis- 
cussed above, we only require a good description of the cen- 
tral regions covered by our proper-motion data. Our best esti- 
mates of the core radius ro, central surface brightness fiv.o, and 
concentration parameter c = login / rp) (for r, the tidal radius 
of the model) are written on Figure^]and also listed in Table 
^above. Note that the fitted c = 2.01 corresponds to a dimen- 
sionless central potential Wo = 8.6 (see Appendix I5t. We es- 
timate uncertainties of about ±25% in ro, ±0.25 mag in ^y.o, 
and ±0.1 dex in c (or ±0.4 in Wo). (Note that these uncer- 
tainties are correlated, such that fits with higher c tend to re- 
quire smaller ro and brighter fiv.o-) Thus, while the new values 
of ro and fiy.o here are formally slightly smaller and brighter 
than the "standard" values in the literature (e. g.. iHairHl 1 996t 
l^well. Guhathakurta. & Gilliland 2000), the differences are 
at the ~ l-cr level. Our core radius is in good agreement with 
the ro = 21" recently derived bv iMapelli et al.. (.2004) . also 
from HST star counts. 
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Give n our neglect of the data at large R> 1 100" when fit- 
ting the iKing I (119661) model in Figure its failure in these 
outermost parts of the cluster is not surprising. In fact, no 
single-mass and isotropic Kine model can satisfactorily fit the 
entire run of surface brightness throughout 47 Tuc. This is 
already well known and forms the basi s of the multimass an d 
anisotropic modeling of this cluster bv lMevla'nljl988lll989h . 

4.1.1. Effect of Mass Segregation 

As we suggested above and will show explicitly in ^ a 
model that assumes velocity isotropy in order to give a de- 
scription of the density distribution over radii R < 1000" is 
appropriate enough as a benchmark in the discussion of our 
proper-motion observations at < 100". More worrisome in 
principle is the assumption of a single stellar mass and the 
implication that stars in any magnitude range should have a 
common spatial structure. This is patently untrue even in the 
core of a relaxed cluster; but even so, it turns out that the 
proper-motion stars in our case do not, in fact, span a particu- 
larly wide range in mass. 

To assess this issue briefly. Table summarizes a theoret- 
ical mass-V magnitude relation for stars from 0.5-0.9Mo. 
In the second column of the table, the absolute magni- 
tudes associated with < 0.8 M0 are those calculated by 
Igaraffe et al. ( 1998) for a 12.6-Gyr old cluster with [M/H] = 
-0.5. For 0.8M(7, < < .9Mf^ we have reUed on the 
models of Bergbusc h & VandenB erg (2001), as illustrated by 
Uriley et al. (2004), for a 12-Gyr o'ld cluster with [Fe/H] = 
-0.83 and [a/Fe] = +0.3. The apparent V magnitudes in the 
third column of Table follow from My given a reddening 
of E(B-V) = 0.04 mag and an assum ed h eliocentric distance 
D = 4 kpc to 47 Tuc (see Table[2 and ^6.3l below). The simple 
formula on the bottom li ne, fo r as a function of V, is our 
own approximation. As ^3.1l discussed in detail, our proper- 
motion sample is limited by the constraint V < 20, which 
corresponds to stellar masses 0.65 < < O.9M0. The 
large majority of the stars for which we can meaningfully ad- 
dress kinematics questions differ in mass by < 20% from the 
main-sequence turn-off and giant-branch stars (V < 17.65, so 
m, ~ O.85-O.9M0) for which we fit the King ( 1966) model 
in Figure [T5I 

King models with high concentration such as the one in Fig- 
urell5lare not substantially different, in their inner parts, from 
bounded (finite central density) isothermal spheres. Crude 
expectations for the relative density profiles at small radii 
for stars of different masses in 47 Tuc can therefore be ob- 
tained by considering models of two -component clusters in 
this idealized limit. Taff et al. (Il975h show that two popula- 
tions of stars with masses 1112 > m\ in an equilibrium isother- 
mal sphere have density distributions related by ptir) / piiQ) = 
[p\{r)l P\(0)T-l"''- As expected, then, p2{r) fafls off more 
rapidly than p\{r) and the heavier stars are more strongly con- 
centrated toward th e cent e r. To quantify this, note that the 
scale ro defined by iKing I (11966 1 see our equation (IB5> be- 
low) is, by construction, nearly equal to the radius at which 
the surface density of an isothermal sphere falls to half its 
central value. As a result, it is also related by a unique nu- 
merical factor to the radius at which the volume density is 
halved. Thus, we have calculated p/p(0) and [p/ piO)]'"-^"^' 
vs. (r/ro) for an isothermal sphere, for a range of fixed mass 
ratios, and located the dimensionless radii at which the di- 
mensionless p and fall (separately) to 1 /2. For modest 

values of 1712/ nii, the core radii of the two mass classes turn 



TABLE 7 

Stellar Mass-Luminosity Relation 



m*/MQ Mv V = (Mv + 13.13) 



0.9 <3.9 <17.0 

0.8 "^5.01 ~18.14 

0.7 6.27 19.40 

0.6 7.58 20.71 

0.5 8.85 21.98 

m,/M0~ 0.9 V<n.O 



~0.9-0.081(V-17) y>17.0 



Note. — Based on the evolution- 
ary models of 'Baraffe et al. ( 1998) and 

IBersbusch & VandenBerEj 1200 IL see also 
Brilev et al. 2004). Apparent magnitudes in Col- 
umn (3) correspond to the absolute magnitudes 
in Column (2) for an assumed distance of £) = 4.0 
kpc to 47 Tuc and a reddening of E(B-V) = 0.04 
mag. 

out to be related roughly by 

ro.i/nu ^ (ma/mi) , (6) 

accurate to better than 1% for 1 < ma/mi < 3 and in good 
agreement with the alternative approximation developed by 
Grindlav et al. (2002, their equation 5a). 

Returning to 47 Tuc, we do not expect equation (|6jl to hold 
exactly; but in the absence of detailed modeling it is useful 
as an order-of-magnitude guide to the situation in the core 
region where all mass classes are expected to be nearly ther- 
malized (with the velocity dispersion of each species scaling 
as ~ m^'). For example, the core radius ro of the faintest 
stars in our proper-motion sample should be only slightly 
larger than that of the turn-off stars with ~ 0.85 Mq: 
(0.65/0.85)~"'^^ = 1.17. Even the "average" cluster stars with 
~ 0.5 Mq should be distributed with a scale radius differ- 
ing by only - 35% [(0.5/0.85)-"-5'* = 1.36] from that of the 
density profile shown in Figure[21 We will refer back to this 
approximate relation on occasion below, when discussing var- 
ious aspects of the velocity distribution of our proper-motion 
stars. First, we describe the exp ected f orm of this distribution 
for an isotropic and single-mass lKingl -model cluster. 

4.2. Projected Velocity Distribution 

As was noted above, the inner regions of King (7966^ mod- 
els with high concentrations (or central potentials) are dy- 
namically quite similar to isothermal spheres. Given that our 
proper-motion sample covers only (parts of) the innermost 
R < 5ro in 47 Tuc, we first ask how different from Gaus- 
sian we might expect the observed velocity distribution to be. 
Very simple considerations suggest immediately that even for 
a sample of ^10,000 stars the departure from Gaussianity is 
somewhat subtle. 

First, King models are defined by a finite escape velocity at 
every radius in the cluster. If W{r) = -<})(/)/ (Jq, where 0(r) 
is the gravitational potential and ao is a parameter closely 
related (though not equal) to the central velocity disper- 
sion, then the model escape velocity at radius r is Vmax/co = 
^2W{r), which is a monotonically decreasing function of ra- 
dius (see AppendixIBlfor more detail). A star moving outward 
from radius r with velocity Vm-iAr) would just reach the clus- 
ter tidal radius with zero velocity. Our proper-motion field 



HST Proper Motions in the Core of 47 Tuc 



19 



contains the center of 47 Tuc, where W(r) = Wo = 8.6 ± 0.4 
for the model fit in Figure Thus, this model implies that 
we should see no stars with total proper-motion speed greater 
than ~ 4.150-0. 

Second, if the velocity distribution were perfectly Gaussian 
and isotropic, then the probability of observing a star with 
relative speed /^tot = (M^ + My)''^^ on the plane of the sky would 
be 



^c(Mtot) = 



Mtot / 2 /--> 2 \ 

— exp(-^f„,/24„) 



2 



Mtot > (7) 



Integrating this, only a tiny fraction e"(4 i5) /2 ^ 1.8 x 10"'* 
of all stars are expected, on average, to have /itot > 4.15(Tiso. 
This shows that simply counting the numbers of stars with 
very high velocities in a sample of 10,000 at the center of a 
cluster like 47 Tuc cannot easily distinguish between the usual 
Gaussian approximation and regular Kine-model behavior 

Nevertheless, the possibility exists that rather more exotic 
physics — such as a central black hole, or a high rate of stel- 
lar collisions — could strongly influence the dynamics, and so 
we would like to assess the significance of any fast-moving 
stars in the core of 47 Tuc as carefully as we reasonably can. 
Appendix 151 therefore develops expressions defining the nor- 
malized probability distribution of observable veloc i ties as 
a function of proje cted c lustercentric radius in King ( 1966) 
models. Equation (IB14t gives Ni{v.\R/ ro), the distribution 
of line-of-sight velocity — or of any one component of proper 
motion — at fixed position R / ro (ro is the cluster core radius, 
determine d by fitting to the surface-density profile as above). 
Equation ( IB17> gives N2(fJ-x , ^^y\R / ro), the joint distribution of 
any two orthogonal components of velocity on the plane of 
the sky [e.g., we could equally well take (/i^, /i,) to be the RA 
and Dec proper motions, or the radial and tangential velocities 
relative to the cluster center]. We will work here mainly with 
the latter function, but note that the operations we perform on 
it can also be applied to the one-dimensional distribution. 

As Appendix E] discusses, isotropy ensures that the the 
proper-motion components /if and /i,, only appear in the 
(position-dependent) velocity distribution N2 in the combina- 
tion (fi^ + ii^.) = We would like to know the total number 
of stars, in a sample of size A/Jot distributed between projected 
radii R[ < R < R2, that lie within a velocity interval d^xd^y 
about the speed /^tot- If S(/?) is the observed projected num- 
ber density of the velocity sample, then this spatially-averaged 
velocity distribution is 



Ir, 



(8) 



which clearly requires the appropriate iKing l-model ro to be 
known. Note that N{^tot) has dimensions of inverse veloc- 
ity squared (A^2 contains a factor ctq in its denominator) and 
satisfies / / Nintoddn^dny = J N(fitot)2'Kfitotdfitot = Mot- It is 
effectively a "density" in the ^x-i^y velocity plane. 

Figure ^] shows the distribution x N(ptot) for A/Jot = 
10,000 stars distributed between Q < R < 5ro a ccording to 
the surface density profile S(/?) of a iKingI (119661) model with 
Wo = 8.6 — meant to recall our proper-motion sample. This is 
the solid curve in the plot. The vertical dotted line marks the 
location of the maximum possible total speed in this model, 
Vmax/co = \/2Wo = 4.15. For comparison, the dashed curve 
traces the isotropic bivariate Gaussian density 

Nc{^Mo^)=:^^^exp(-^liJ2al^ , (9) 
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Fig. 16. — Comparison of the projected, two-dimensional velocity distri- 
bution averaged over the inner 5 core radii of a Wq = 8.6 King 1 1966) model, 
vs. an isotropic, bivariate Gaussian distribution with the same average disper- 
sion. Both theoretical curves are normalized to a total population of 10,000 
stars, similar to our proper-motion sample. Significant differences between 
the two appear only at high speeds /^lot ^ 2(to, above which relatively few 
stars are expected in either case. The vertical, dotted line is at /imt = 4.15cro. 
the maximum velocity allowed at the center of a Wq = 8.6 King model. The 
Gaussian shown here would predict <1 star in 10,000 to be found above such 
a speed. 

for CTiso = 0. 923(To equal to the density-weighted avera ge o f 
the proiected lKind -model dispersion (cr?, defined in eq. IB9I ) 
over the inner five core radii. Evidently, the differences be- 
tween the two curves are not easily visible — let alone statis- 
tically significant — until /itot ^ 3(To- The Gaussian distribu- 
tion would predict only e"'*^ x A/Jot — 1 1 1 stars at such high 
speeds; the .King model, fewer still. 

4.2.1. Accounting for Measurement Errors and Irregular 
Spatial Distributions 

While informative, the idealized calculations just above still 
lack any recognition of the facts that observed velocities are 
subject to measurement error, and that a kinematical sample 
of stars need not trace the overall stellar density profile (as 
our proper-motion sample clearly does not; e.g., see Figure 
above). Both of these complications will affect the over- 
all shape of the King-model distribution N(fitot) and alter the 
relative strength of its wings. 

To deal with these issues, consider a set of stars with proper- 
motion measurements and uncertainties {fix,i,IJ-y.i,^x,i,\:i}, 
and clustercentric radii {/?,}, where ; = 1 , . . . , A/Jot- For each 
individual star we know the i deal ' King^-model distribution 
N2{lJ.x, IJ-y\Ri/ fo) from equation IB 171 and the bivariate Gaus- 
sian distribution Pi{6iix, Sfiy) of velocity measurement errors. 
The predicted observable probability density at is then the 
convolution of N2 and P,-. Strictly speaking, the velocity er- 
rorbars Ai , and A^,_,-, which appear in P,, are not exactly 
equal for all stars, but in general they are nearly so. We 
therefore define (for these purposes only) a single errorbar 

A; = ( A^ ; + A^ ,) /2 for both velocities of each star, so that 
the approximate measurement-error distributions 



Pi(Sfix,Sfiy) = 



1 



27rA, 



■ exp 



(J/i,)^ + ((5/i/ 
2A? 



(10) 
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are isotropic. Then the error-convolved velocity distribution 
at star /, 

N^{^JL„^iy\Ri|rQ,'Ei) = 

_ J 

(11) 

remains isotropic as well, i.e., it still depends on /ij and /i,, 
only through /itot- 

Although cumbersome, equation il l\ is straightforward to 
compute for every star once ro is known. (Note that the in- 
tegrals are not actually infinite, since the original N2 in the 
integrand vanishes for proper-motion speeds above the local 
maximum velocity.) Then, rather than assume that the proper- 
motion stars accurately follow the self-consistent Kins -model 
surface-density profile, we simply find the velocity distribu- 
tion for stars distributed on the sky exactly as the observed 
ones: 

A/tot 

A^(Mtot) = ^ A^2*(Ai.o Mvl/^z/ro, A,) . (12) 
1=1 

It is this last equation (rather than the idealized eq. |8l) that we 
use to produce models for comparison against the observed 
proper-motion distributions in the next Section. 

5. THE OBSERVED PROPER-MOTION VELOCITY 
DISTRIBUTION 

5.1. One-Dimensional Distributions 

The proper motions in Table |5] are given in terms of their 
RA (positive Eastward) and Dec (positive Northward) com- 
ponents, which we write as (/i^j/ia). In a dynamical context, 
however, it is more meaningful to refer the velocities to the 
cluster center Thus, if the coordinates of any point are writ- 
ten in terms of the projected clustercentric radius R and the 
azimuth 6 (measured in degrees North of East), it is a simple 
matter to make the transformation 



Mfi = Ma cos 8 + fis sin Q 
/iQ = -/ic sin 8 + /i5 cos 



(13) 



The uncertainties in and /ie then follow from the usual 
rules: 



A2 



Al cos2 8 + A|sin2 8 
sin^8 + A2-cos2 8 



(14) 



Note that /ig is positive for motion outward from the cluster 
center on the plane of the sky; /ie is positive for clockwise 
motion; and by the right-hand rule, the line-of-sight velocity 
V; is positive for motion away from the observer 

To investigate the distribution of these velocities, we first 
look at the entke sample of all "good" proper motions — the 
12,974 stars in Table|5]which have P(x^) > 0.001 for both the 
fitted and fis, and magnitudes V < 20 (see ^3.1> . Figure 
[TtI displays a number of histograms of the one-dimensional 
velocities for these stars. The top panels give the separate 
and lis distributions, the middle panels present those for 
the /iR and /ie components, and the bottom panel shows the 
histogram of speeds /itot = (Pa + A^s)'''^ = (Pr + /"e)'''^- For 
reference, speeds in mas yr"' are converted to km s"' by 

1 mas yr"' = 4.74(D/kpc) km s"' (15) 

for a distance D to 47 Tuc. With D ~ 4 kpc (Table 
1 mas yr~' ~ 19 km s~' . 



We have calculated the first two moments of each of the 
first four velocity distributions in Fig.[n] All of the means are 
consistent with zero — as they should be, given that we have 
used each star's nearest neighbors to derive velocities relative 
to the net cluster motion on the sky (see ^2.3l and Appendix 
IaI . We have estimated the true velocity dispersion in each 
component /i^, fig, and fiQ, using equation ( IA14> : 



1.1 
45 



obs " 



i=l 



— Va2 



(16) 



where A/Jot = 12,914-; A, is the uncertainty in the proper- 
motion component /i,; and the leading factor corrects an un- 
avoidable error introduced by our local transformations (as 
discussed in AppendixlAl. 

The true one-dimensional dispersions computed from this 
formula are listed in the appropriate panels of Figure ^\ 
The uncertainty in each is about ±0.004 mas yr"', and 
thus the four different proper-motion components are all 
consistent with a single, average (T,jue = 0.578 mas yr"' = 
11.0(Z)/4kpc) km s" . This is the first and simplest direct 
evidence that the velocity distribution in the core of 47 Tuc 
is (as expecte d) iso tropic. This conclusion is justified more 
rigorously in ^6.21 below, where we more carefully examine 
velocity dispersion as a function of stellar magnitude and clus- 
tercentric radius. 

Meanwhile, each of the top four panels in Figure ^\ has 
three curves drawn in it. The dotted curves are Gaussians with 
dispersion equal to (1 + 1 /45)'/^ x (Ttrae, while the much nar- 
rower, dashed curves are the measurement-er ror d istributions 
given by the sum of Gaussians in equation JAIL scaled up 
by the total number of stars. The thick, solid lines represent 
the convolution of the "intrinsic" dotted Gaussians with the 
normalized error distributions. These convolutions are them- 
selves the sum of Gaussians. They closely follow all of the 
observed histograms, showing that the intrinsic velocity dis- 
tribution in the core of 47 Tuc is in fact very nearly Gaussian. 
This is not assumed in the derivation of eq. jl6> in Appendix 
|A]but, from the discussion in ^4.21 it is just what we expect 
even if the cluster actually has the structure of a iKina til 96^ 
model. 

The curves in the bottom panel of Figure [n] assume the 
velocity isotropy suggested by the four upper panels. With 
the measurement-error distribution given by the sum of Gaus- 
sians, if the intrinsic velocity distribution were also Gaussian 
then the expected distribution of observed speeds would be 
the convolution (cf. eq. Q) 



Mo.Pa(/i,oO = V^^exp| 



Mtot > 
(17) 

where A,- = ( A^ + A| ,) /2 is the isotropized velocity uncer- 
tainty of each star. The dotted curve in the bottom panel of 
FigureFTTlis this equation evaluated with (tiso = (1 + 1/45)'/^ x 
0.578 but A/ = for all stars, i.e., it is the Gaussian approx- 
imation to the intrinsic distribution of speeds in our proper- 
motion sample (with the true dispersion inflated by trans- 
formation error). Conversely, the dashed curve is equation 
( I17> with CTiso = but using the set of {A,} actually mea- 
sured for the stars in our sample; this is the distribution of 
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Fig . 17. — One-dimensional distributions of velocities in the RA and Dec directions and in the projected radial and azimuthal directions relative to the cluster 
center, as well as the distribution of speeds ^lot = (At^ + Atl)'^^, all f or the sample of 12,974 "good" proper motions. Bold, dotted curves in all panels are Gaussians 
with the error-corrected dispersions crtme indicated (from eq. fT^ ). while nan'ower, long-dashed curves ai'e the velocity measurement-error distributions. Bold, 
solid curves are the convolutions of the two. The good agreement with the data in all cases indicates that the intrinsic velocity distribution in the core of 47 Tuc 
is nearly Gaussian, at least for proper-motion speeds /ito( ^ 2-2.5 mas yr"' . Isotropy is also indicated by the similarity of crirue in the four directions shown. 



speeds that would obtain if all apparent motion were due 
solely to measurement error The heavy solid curve run- 
ning through the histogram is the full equation ( I17> . with 
(7;,^ = (1 -I- 1 /45)'/2 X 0.578 and each A,- as observed. The 
main effect of velocity measurement error is perhaps most 
clearly seen by comparing the dotted and solid curves in this 
representation of the data: it works to take power from the 
peak of the intrinsic velocity distribution and move it out to 
the wings. 

5.2. Two-Dimensional Distributions 

The distributions in Figure^jare, of course, averages over 
the irregular spatial distribution of our velocity sample and 
a wide range of stellar magnitudes, 14 < V < 20. We expect 
the velocity distribution to vary with clustercentric radius and, 
if mass segregation and energy equipartition have taken hold 
to any extent, with stellar mass. It has to be emphasized in 
particular that the average dispersion crtme = 0.578 mas yr"' 



is only an indicative number and is not suitable, e.g., for di- 
rect comparison with the radial-velocity dispersion of giant- 
branch stars to estimate a distance to the cluster 

We therefore look now at the velocity distribution for stars 
in a number of distinct bins in magnitude and radius. To 
do this, we take advantage of t he isotropy indicated by Fig- 
ure [n] and confirmed in ^6.21 below, to construct the two- 
dimensional joint distribution of (/iQ,,/!^), or equivalently 
(MSiMe)^ rather than the one-dimensional histograms for in- 
dividual components. We also switch to compare the data 
explicitly with the velocity distributions expected in a iKind 
ri966) model cluster with a finite escape velocity. 

For these purposes we first restrict the sample of "good" 
proper motions from Table |5l further, by imposing the color 
criterion 



1.5 < (U-V)<2.6 , 
0.6 < (U-V) < 1.9 , 



V < 16.75 

V > 16.75 



(18) 
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Fig. 18. — Two-dimensional velocity distributions for stars at all observed 
clustercentric radii, divided into four magnitude bins. The total number of 
good proper motions and the oberved, error-corrected (Jtruc in each bin are 
given in each panel, as is the predicted single-mass, isotropic King -model 
dispersion accounting for the stars observed positions. The dotted curve in 
each case is the intrinsic mo del v elocity distribution convolved with the ob- 
served spatial distribution ( j|4.2t . The bold curve is this model convolved 
with the velocity measurement-eiTor distribution, and the dash-dot curves are 
its 68% confidence limits. Note the excess stars in the high-velocity tails, and 
the fact that the kinematics at V > 19.5 are error-dominated in general. 



which is designed to exclude blue-straggler stars and extreme 
outliers from the main cluster sequence in the CMD of Figure 
[Tol This leaves us 12,791 stars to work with. Given this sam- 
ple, we define a series of total proper-motion speeds {/itotj} 
and choose a range of stellar magnitude and clustercentric ra- 
dius to examine. We calculate the observed speed for each 
star in the specified range of V and R; count the number which 
fall in each interval [/itot.>,Mtot.>+i); and divide the population 
of each velocity interval by the "area" TriPi^tj+i - Mtotj)- The 
result is an isotropized "density" profile in projected veloc- 
ity space, Niiitot), which can be compared directly to model 
distributions of the type discussed in ^4.21 Uncertainty in the 
observed A^(/itot) is estimated by generating 1000 artificial sets 
of velocity data from the subsample of real stars in any magni- 
tude/radius bin (using a standard bootstrap with replacement; 
e.g., Ifress et al. 1992, Section 15.6) and calculating N(^tot) 
for each artificial set to find the 68% confidence intervals on 
the actual distribution. 

5.2.1. N(fitot) vs. Magnitude for Stars at all Radii R < 100" 

FigureFTSlshows the observed N(fitot) for the proper-motion 
sample of 12,791 stars over our entire field, divided into four 
broad magnitude bins. The curves in all panels of the figure 
are the theoretical distributions for a Kine ( 1966) model with 
Wo = 8. 6 an d ro = 20."84, computed as described in Appendix 
inland ^4.21 In order to normalize these curves to t he da ta, the 
model velocity scale co must be known (see eq. IB 171 ). and 
we evaluate this as follows. 

Given Wo and ro, we have computed the dimensionless pro- 
jected velocity-dispersion profile (J'^^^iR / ro) / ctq (given as cr? 
in eqs. IIB8I and IIB9I ). Then it is straightforward to calculate 
the predicted iKind -model velocity dispersion for any collec- 
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(19) 



At the same time, we can estimate the true projected velocity 
dispersion in any direction on the plane of the sky, corrected 
for both measurement and transformation errors, from equa- 
tion M6\ . Doing this separately for both the RA and Dec (or 
both the R and <d) components of proper motion yields the im- 
proved value cTf^yg = (Cn-ue.a + ^trae sV^^ which we have calcu- 
lated for the stars in each of the magnitude bins in FigurelTSl^ 
In particular, for the stars with V < 18.5 throughout our en- 
tire field (the subsample shown in the top left panel of Figure 
we find (Jtrue — 0.583 i 0.004 mas yr , while using their 
positions in equation jl9t gives (J^^g^/(jl = 0.849 ±0.015. Re- 
quiring that (Tpi-ed = cTtrae for thcsc stars then implies 

do = 0.633 ±0.010 mas yr"' = (12.0±0.2)(D/4kpc) km s"^ . 

(20) 

Note that this quantity really is only a normalization factor — 
nowhere in a King model does the velocity dispersion actually 
equal ctq — and that its value here is specific to the values W) = 
8.6andro = 20."84. 

Each panel in Figure ^| lists the error-corrected Ctrue for 
the stars in question, and dpied from equations ( I19l l and ( I20t . 
The errorbars on (Ttrue are 68% confidence intervals obtained 
numerically from the same bootstrap experiments used to 
estimate the uncertainties in the full velocity distrib utions. 
The dotted curve in each panel is the intrinsic iKind -model 

Ni^itoi) = I]^°'A^2(MA-,Mvl^i/'"o), which is free of velocity- 
measurement errors but is calculated using an inflated do = 
(1 + 1 /45)'/^ X 0.633 = 0.640 mas yr"' to account for the er- 
ror caused by local coordinate transformations (see ^2.3.2l and 
Appendix|A|l. The solid line in each case is the model curve 
after convolution with our measurement errors (see ^4.2> . and 
the dash-dot lines connect its one-sigma uncertainties (which 
reflect Poisson noise in the number of stars predicted in each 
of the velocity bins that we defined to construct the observed 
distribution). 

The main conclusion to be taken from these graphs is that, 
after convolution with the measurement errors, the isotropic 
and single-mass model velocity distributions describe the 
proper-motion data as a whole quite well for speeds /itot ^ 2- 
2.5 mas yr"'. This is essentially a re-statement of the good 
agreement with Gaussians at these velocities in Figure [n| 
The poorest agreement is shown by the faint stars in the lower 
right-hand panel of Figure ^| where the observed A^(/itot) 
points tend to fall above the model curve for /itot ^ 1 mas 
yr"' and slightly below it at larger speeds. This simply re- 
flects the discrepancy between the estimate of dtrae and the ex- 
pected (Tpred for these stars, but this is in turn mostly due to the 
fact that — as either Figure ^| or Figure |9| above shows — the 
observed velocity distribution for V > 19.5 is strongly domi- 
nated by measurement errors. 

From our estimate of ctq in equation ( I20t and the fact that 
the m aximiim total velocity for bound stars at the center of a 
iKind 119661) model is Vmax/co = \/2Wo, if our Wo = 8.6 model 

' To prevent contamination by extreme velocity outliers, we estimate uuue 
using only stars with total speeds satisfying /itoi < 3.0 mas yr"' if V < 19; 
/itot < 3.2 mas yr"' if 19 < V < 19.5; and fj.tot < 3.6 mas yr"' if 19.5 <V < 
20. The upper limit increases with magnitude in recognition of the larger 
measurement errors for fainter stars. 
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Fig. 19. — Same as Figure fTsl but only for stars in the innermost R < 20", 
where velocity uncertainties are smallest. Note the constancy in o-prej as a 
function of stellar magnitude, reflecting the single-mass assumption of the 
models. This contrasts with the rise in the observed crtmc toward fainter V, 
presumably as a result of mass segregation and equipartition of energy. 

for 47 Tuc were fully correct we would naively expect to see 
few or no cluster members with true speeds /itot ^2.6 mas 
yr"'. Indeed, all the dotted curves in Figure fall to zero 
at about this speed. Faster-moving stars are clearly found in 
our sample, however. At relatively faint magnitudes V >19, 
most of these are clearly attributable to measurement error, 
which inevitably scatters some slow stars to very high appar- 
ent speeds. Thus, although there are a few individual excep- 
tions, the majority of the NjjMot) measurements in the tails of 
the lower panels of Figure^]are statistically consistent with 
the error-convolution of the model curves with /imax — 2.6 
mas yr"'. The same is basically true for the brighter bin 
18.5 < y < 19. It is only for the brightest subsample of stars 
with V < 18.5 — for which the Kin2 model we have fit should 
in fact be most appropriate — that a statistically significant ex- 
cess of high-velocity stars is clearly evident. 

We discuss these high-velocity stars in more detail in ^5.31 
just below but note here that, their specific properties aside, 
there are very few of them indeed. Over all magnitudes V < 
20, only 46 of the 12,974 good proper-motion stars in Table|5l 
have observed /itot > 2.6 mas yr"', and fewer than 10 of these 
can currently be said unequivocally not to be due to velocity- 
measurement error There is, then, little to no evidence for a 
substantial nonthermal population of stars in the core of 47 
Tuc. 

5.2.2. Nifitot) vs. Magnitude in the Well-Measured 
Innermost 20" 

The lKinal -model CTpied for the four magnitude bins of Figure 
1 181 are very similar to each other because the model assumes 
that stars of all magnitudes have the same mass and thus the 
same spatial and velocity distributions (any slight differences 
in cTpied are solely due to differences in the spatial distributions 
of the observed stars in the subsamples). What is perhaps 
slightly surprising is that the observed au-ue are also, with the 
exception of the faintest, most error-prone stars, very similar 
both to each other and to the single-mass model dispersions. 



From Table Q the stellar mass decreases by ~ 0.04Mo for 
each 0.5-mag increase in V. Timescales in the dense core 
of 47 Tuc are short enough that equipartition of energy has 
almost certainly been established (Mevlan 1988, 1989), in 

-1 /2 

which case we should expect to find CTtrue ~ w^* and thus a 
systematic increase of roughly 3% in velocity dispersion from 
panel to panel in Figure The fact that this is not observed 
suggests that the effect may have simply been blurred out by 
velocity measurement errors. This is further indicated by the 
significantly — and unphysically — lower value of CTtrue in the 
faintest bin, 19.5 < V < 20. 

Thus, in Figure^]we show the observed and model N(^tod 
in the same four magnitude bins as Figure^] but now for 
3,097 stars in just the innermost region of our field, R < 20". 
As was discussed in ^3.1.21 this is where the velocity uncer- 
tainties are smallest for any stellar magnitude (see Figure|9jl. 
The error-corrected CTtrue in each magnitude bin is again shown 
in each panel of Figure ^] along with the model CTpred and 
Nintoi) curves for ctq = 0.633 mas yr"', all as above. The ex- 
pected increase of CTtrue toward fainter magnitudes is in better 
evidence now, and causes the observed A^(/itot) in the faint bins 
to be broader than the single-mass model dis trib utions. Com- 
parison of these results with those in Figure^]suggests that 
it is this central-most region, with its relatively precise proper 
motions, that will provide the most useful kinematical con- 
straints on detailed multimass dynamical models for 47 Tuc. 

5.2.3. N(^toi) vs. Radius for the Brightest Stars 

Finally, Figure |20| presents the velocity distribution of the 
6,342 brightest proper-motion stars in four concentric circular 
annuli covering our observed field. From Table0 the magni- 
tude range V < 18.5 chosen here corresponds to stellar masses 
OT* ~ 0.77-0.9Mq. The mean magnitude is (V) ~ 17.7, and 
thus (ot*) ~0.84Mq — essentially the main-sequence turn-off 
mass, and thus very similar to the stars used to constrain the 
iKing l model fit in ^4.11 These are also the stars for which 
the velocity errorbars are lowest (Fig.|9}. In fa ct, for clarity 
in Figure we only show the error-convolved iKind -model 
distributions and their one-sigma confidence bands, as these 
distributions are not grea tly different from the intrinsic ones 
(dotted curves in Fig.Tl8land Fig.ll9ll. 

Perhaps the most striking aspect of Figure |20| is that the 
estimates of CTtrue come within about 2% of the single-mass 
King -model CTpred at all radii. That is, the kinematics of stars 
with the main-sequence turn-off mass are described to rather 
high accuracy by a dynamical model which assumes that their 
spatial distribution traces the distribution of total mass in the 
cluster — even though this includes heavier neutron stars and 
white dwarfs and a large number of much less massive main- 
sequence stars, all stratified by mass segregation into differ- 
ent relative density profiles. At the same time, however, the 
slight disagreement between the observed and predicted ve- 
locity dispersions is informative. Looking at the top two pan- 
els in particular, CTtrue stays essentially constant as the cluster- 
centric radius increases out to 50", while CTpred falls by just 
over 4%, an amount which exceeds the observational error- 
bars by a factor of ^^3. Thus, the stars in 47 Tuc are mov- 
ing in a net gravitational potential which is less centrally con- 
centrated than that generated by the turn-off and giant-branch 
stars alone. 

This implies that the core of 47 Tuc cannothe dominated by 
heavy stellar remnants such as neutron stars or white dwarfs, 
which would have a significantly smaller core radius than 
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Fig. 20. — Similar to Figures [l8land[T9l but now showing the velocity 
distribution for bright stars only, in four intervals of clustercentric radius. 
Measurement errors are small at these magnitudes (Figure|9|, so for clarity 
only the fully eiTor-convolved model distributions and their i-cr confidence 
bands are drawn. The somewhat slower fall-off of observed crtme with radius, 
relative to the model cTpred implies that these main-sequence turn-off mass 
stars are moving in a potential that is more extended than the one generated 
self-consistently by their spatial distribution. Thus, the core of 47 Tuc appears 
to be dominated by average (fainter) main sequence stars, rather than by more 
spatially concentrated heavy remnants. 

the bright stars and lead to a steeper drop-off in velocity 
dispersion with clustercentric radius (for neutron stars with 
m* = 1.4M0, eq. |6l gives ro ~ 16" vs. the 20"84 measured 
for m, ~ 0.85-0.9Mo). This conclusion is consistent with 
other lines of argument. The multimass and anisotropic mod- 
els of Mevlan (1988, 1989) require, essentially to obtain a 
good fit to the detailed shape of the surface-brightness pro- 
file in Figure [21 above, that heavy stellar remnants make up 
only ~0.1% or less of the current total mass of 47 Tuc. The 
most recent theoretical and observational estimates of the to- 
tal neutron-star population in 47 T uc he in the range Mi-am ^ 
300-1500, easily compatibl e withiMevlanfs indirect in ference 
Jivanova. Fr egeau. & Rasi oll2005HHeinke et alJ2005l) . But in 
this case, even allowing for mass segregation to place most 
of the remnants in the core leads to a structure in which nor- 
mal stars with < O.85M0 dominate the gravity in the in- 
ner few core radi i by a t least an order of magnitude (see Fig- 
ure 5 of Mevlanl ll988l) . For an average stellar mass around 
m» ~ O.5M0, we expect the total mass distribution to have a 
large core with ro ~ 27 "4 (from eq. again), which is then 
consistent with the shallow radial dependence of the proper- 
motion fTtnie i n Fi gure 1201 This will also be seen, more quan- 
titatively, in ^6.5l below. 

A second point to be taken from Figure|20lis its demonstra- 
tion that the "high-velocity" stars — those which contribute to 
N(fi(ot) measurements which lie significantly above the tails 
of the .King -model curves — are not confined exclusively to 
the very center of 47 Tuc but can be found out to five core 
radii away. This is, in fact, expected for proper-motion high- 
velocity stars, whereas such stars observed in radial velocities 
are all expected to be observed close to the cluster center in 
projection. We now present more details on all of the very 
fast-moving stars in our full sample. 



5.3. High-Velocity Stars 

As was discussed abo ve, the nominal central escape veloc- 
ity for our ' Kind (1196 6') model of 47 Tuc is Vmax = V^Wo x 
(Jo — 2.63 mas yr"', or 49.8(D/4kpc) km s~'. The presence 
of stars in our sample with observed proper-motion speeds 
above this value is due in part to unavoidable measurement 
error (although Figure^] shows that this cannot account for 
all of them), and at some level it may also reflect the inevitable 
limitations of single-mass modeling (as in Figure^] fainter 
stars can attain broader velocity distributions and still remain 
bound to the cluster). Nevertheless, such stars are also of po- 
tential physical interest as the possible byproducts of close 
stellar encounters and strong scattering in the very dense core 
environment. 

Thus, Table |8] lists the basic properties of all the stars in 
Tabled with P(x^) > 0.001 in both the RA and Dec com- 
ponents of proper motion, V < 20, 0.2 < (U-V) < 2.6, and 
/itot > 2.6 mas yr~'. There are only 46 such stars in our full 
sample of 12,974 "good" proper motions. FigureOTIidentifies 
these stars relative to the rest of the velocity sample in plots 
of their distributions in the projected velocity planes (/ic/ia) 
and (/iR,/ie)i their {ctioooj^im)) positions relative to the clus- 
ter center, and a{V,U- V) color-magnitude diagram. 

We have inspected the position vs. time data for each of the 
high-velocity stars in Table|8]and Figure|^individually, and 
separated them into four categories: 

(1) Marginal: 26 stars have total speeds just barely above 
2.6 mas yr"', and/or velocity uncertainties that could bring 
them below fitot < 2. 65 m as yr"'. Table|S]and the lower right- 
hand panel of Figure|^show that the large majority of these 
stars are fairly faint, with V > 19, suggesting that their (mod- 
estly) high apparent speeds exceed the v^ax limit simply be- 
cause of measurement error. 

(2) Low P(x^): Six stars with Pix') > 0.001 in both ve- 
locity components nevertheless have P(x^) < 0.05 in at least 
one component. Thus, it can't be said with better than 95% 
confidence that these stars follow genuinely linear motion, let 
alone whether their fitted velocities are completely reliable. 
Most of these stars are fairly bright, and additional data from 
future epochs could either confirm or remove them as high- 
velocity candidates. 

(3) SMC: Five stars are grouped together at = 4.3-5.3 
mas yr"' West and fig = 0.6-2.0 mas yr"' North in the upper 
left-hand panel of Figure |^but are not so closely associated 
with each other in the (/ifi,/ie) plane which refers their mo- 
tion to the center of 47 Tuc. The Small Magellanic Cloud is 
in the background of our proper-motion field, and its average 
motion relative to 47 Tuc is (p.a,iis) — (-4. 7, +1.3) mas yr"' 
(Anderson & King 2003a). Thus, we identify these five stars 
as red giants in the SMC. Their faint V magnitudes and rel- 
atively red (U - V) colors (see the lower right-hand panel of 
Figurel2ll are consistent with this interpretation. 

(4) Certain: Nine stars remain which pass all of the first 
three tests, and which we therefore take to be genuinely high- 
velocity stars that could well be associated in some way with 
47 Tuc. 

In fact, all of the "certain" high-velocity stars have /itot > 
3.0 mas yr"' = 56.9(Z)/4kpc) km s"', which is obviously a 
lower limit to their total speed including motion along the 
line of sight. This makes it quite unlikely that any dynam- 
ical model could allow for them to be bound to the cluster 
currently. If they were ever members of 47 Tuc, they must 
have been produced by extreme events inside the cluster core. 
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Fig. 21. — Distribution of "high-velocity" proper-motion stars (ptot > 2.6 mas yr"', indicated by the dotted circles in the upper panels and corresponding 
roughly to the nominal escape velocity at the center of 47 Tuc), in velocity space, in position relative to the cluster center, and on a color-magnitude diagram. See 
the text and Table Islfor details, but note that these fast-moving stars are not confined to the very center of the cluster and show no preference for one particular 
direction of motion. None is a blue straggler. 



On the other hand, it may appear suspicious that six of the 
nine stars are gathered at t he b ottom of the main sequence 
in the lower right of Figure |^ By construction of our cat- 
egories, the relative velocity uncertainties on these particu- 
lar faint stars are not extreme, and their speeds seem secure 
enough; but this positioning on the CMD suggests the pos- 
sibility that at least some of them might be foreground field 
stars (of which we do expect to see a handful throughout our 
field; cf. the contamination estimate in Table[0. 

Of the three brighter "certain" high-velocity stars, two 
(M058939 and M053721, both at the main-sequence turn-off) 
were previously identified by McLau2hlin et al. (2003), who 
also pointed out that they have ^ |/ie| and thus appear 
to be moving almost purely radially, with respect to the clus- 
ter center, on the plane of the sky. This is true of only one 
other "certain" star (M054717), although it is much fainter 



and falls blueward of the main cluster sequence in the CMD 
of Figure]^ Even among the marginal and low-P(x^) stars 
in Table |8] only a handful show such a preference in their di- 
rection of motion. Since we do not have information on the 
line-of-sight velocity component for any of these stars, it is 
impossible to say exactly how their total space motions relate 
to the cluster center. But with this caveat, even if all of the 
non-SMC stars with apparent /itot > 2.6 mas yr"' are bona- 
fide high-velocity stars and genuinely associated with 47 Tuc, 
only a small minority show any obvious hint that they might 
have been ejected directly outward from the very center of the 
cluster (as could be expected to happen, for example, as a re- 
sult of interactions with centrally concentrated binaries or a 
large black hole). 

Our fastest stars are rather more extreme than the two 
"high-velocity" stars found in the radial-velocity survey by 
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TABLE 8 

Stars with Proper-Motion Speed fim > 2.6 mas yr"' 









^tot 


0^2000 


52000 


Ma 


^J■s 




Me 




ID 


y 


(71 — V\ 
i^U V ) 


[nicis yr ] 


[arcsec] 


[mas yr ] 


[mas yr ] 


comment 






('<\ 




(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


M 1 ^ 


M055901 


15.33 


1.97 


3.15 


-7.25 


-2.19 


-1.71 ±0.12 


-2.64 ±0.10 


2.40 ±0.12 


2.03 ±0.10 


low P(x^) 


M077141 


17.15 


1.49 


2.75 


4.86 


20.91 


-0.68 ±0.20 


2.66 ±0.20 


2.44 ±0.20 


1.27 ±0.20 


marginal 


M005828 


17.15 


0.85 


3.45 


17.39 


-71.22 


l.I8±0.27 


3.24 ±0.38 


-2.87 ±0.38 


1.91 ±0.28 


low P{x^) 




17.23 


0.87 


2.86 


-70.65 


5.85 


0.93 ±0.19 


-2.70 ±0.21 


-1.15±0.19 


2.61 ±0.21 




M027453 


17.46 


0.75 


2.65 


44.63 


-35.47 


-0.40 ±0.19 


2.62 ±0.37 


-1.94 ±0.27 


1.80±0.31 


marginal 


M058939 


17.66 


0.78 


3.19 


-12.25 


0.96 


-3.19±0.14 


-0.08 ±0.10 


3.17±0.14 


0.33 ±0.10 


certain 


M053721 


17.70 


0.85 


3.04 


-66.05 


-4.50 


3.04±0.18 


0.05 ±0.15 


-3.04±0.18 


0.16±0.15 


certain 


M057847 


17.81 


0.85 


2.64 


12.79 


-0.15 


-2.01 ±0.25 


-1.71 ±0.23 


-1.99 ±0.24 


-1.73 ±0.23 


marginal 


M049017 


17.93 


0.75 


3.45 


-8.75 


-9.45 


3.45 ±0.25 


0.08 ±0.16 


-2.40 ±0.21 


2.48 ±0.22 


low Pi-Y^) 


M073236 


18.20 


0.75 


2.71 


3.87 


16.45 


-2.59 ±0.29 


-0.80 ±0.22 


-1.37 ±0.23 


2.34 ±0.29 


marginal 


M047472 


18.28 


0.91 


3.07 


-7.18 


-11.18 


1.61 ±0.20 


-2.61 ±0.19 


1.33 ±0.20 


2.77 ±0.20 


certain 


M003793 


18.90 


2.17 


4.61 


-5.60 


-76.99 


-4.57 ±0.48 


0.59 ±0.54 


-0.25 ±0.54 


-4.60 ±0.48 


SMC 


M029657 


19.02 


1.09 


3.15 


63.47 


-32.54 


2. 16 ±0.45 


-2.29 ±0.60 


2.97 ±0.48 


-1.05 ±0.57 


marginal 


M067102 


19.06 


1.95 


4.77 


-38.13 


9.66 


-4.33 ±0.71 


2.01 ±0.81 


4.69 ±0.72 


-0.89 ±0.80 


SMC 


iVIU J / jUo 


1 Q OS 


1 AO 


J. J 1 


58.57 


-22.83 


1.49 ±0.57 


3.18±0.46 


0.23 ±0.55 


3.50 ±0.47 


low ) 




1 Q DQ 
1 y.vy 


1 24 




-34.14 


-19.69 


2.53 ±0.44 


-1.63 ±0.53 


-1.38 ±0.47 


2.68 ±0.51 


marginal 


M049738 


19.21 


1.27 


3.27 


60.06 


-8.68 


-0.90 ±0.55 


3.14±0.58 


-1.34 ±0.55 


2.98 ±0.58 


low P(x^) 


JVKMUiUy 


19.38 


1.42 


2.61 


-41.94 


-19.46 


2.57 ±0.71 


0.44 ±0.52 


-2.52 ±0.68 


0.68 ±0.56 


marginal 


MUliiUU 


19.40 


1.23 


3.19 


16.25 


-56.91 


-0.23 ±0.54 


-3.18±0.69 


2.99 ±0.68 


-1.10±0.55 


marginal 




19.42 


1.95 


2.65 


49.89 


-30.12 


2.45 ±0.83 


-1.02 ±0.90 


2.62 ±0.85 


0.39 ±0.88 


marginal 


JVlU/5jy5 


19.44 


1.83 


3.51 


1.45 


18.89 


2.47 ±0.73 


-2.50 ± 1.05 


-2.30 ± 1.05 


-2.65 ±0.73 


marginal 


M061 148 


19.47 


1.60 


3.83 


-68.89 


3.26 


1.92 ±0.69 


-3.31 ± 1.12 


-2.07 ±0.69 


3.22± 1.12 


certain 


JVl(JU34(jy 


19.49 


1.69 


3.14 


40.08 


-78.17 


3.05 ±0.97 


0.73 ±0.74 


0.74 ±0.79 


3.05 ±0.93 


marginal 


MO/z/OU 


19.50 


1.99 


2.72 


-56.36 


15.83 


2. 15 ±0.82 


1.67 ±0.70 


-1.62±0.81 


-2.19±0.71 


marginal 


TV jr^T /1 1 'in 

JVlUJ412y 


19.52 


1.45 


2.80 


40.97 


-26.76 


-1.27 ±0.78 


-2.50± 1.01 


0.30 ±0.85 


-2.79 ±0.95 


marginal 




19.55 


1.85 


4.79 


35.85 


-75.30 


-2.00 ±0.68 


-4.35 ±0.91 


3.07 ±0.87 


-3.68 ±0.73 


certain 


MU47 7 zi 


19.57 


1.34 


3.37 


57.87 


-10.89 


1.67±0.61 


-2.93 ±0.53 


2.18±0.61 


-2.57 ±0.54 


certain 


TV JA 1/1111 

JVlU14zlz 


19.57 


1.56 


3.36 


20.07 


-55.21 


2.07 ±0.63 


2.65 ±0.65 


-1.78 ±0.65 


2.85 ±0.63 


certain 


Muu/uyy 


19.58 


1.85 


3.39 


-0.08 


-68.26 


0.58 ±0.77 


-3.34 ±0.84 


3.34 ±0.84 


0.58 ±0.77 


marginal 


MU4Uo5 1 


19.65 


1.44 


4.81 


58.89 


-18.96 


-4.61 ±0.40 


1.37 ±0.39 


-4.81 ±0.40 


-0.11 ±0.39 


SMC 


MUzi/oj 


19.67 


1.67 


2.95 


12.94 


-40.61 


-1.09 ±0.52 


-2.74 ±0.61 


2.28 ±0.60 


-1.87 ±0.53 


marginal 


M()54717 


19.68 


1.19 


3.14 


14.08 


-3.44 


3.14±0.26 


-0.04 ± 0.36 


3.06 ±0.27 


0.71 ±0.36 


certain 


A ^ A /I ACT7 

JV1U4UjZ/ 


19.74 


1.55 


5.12 


-50.95 


-19.10 


-5.04 ±0.52 


0.92 ±0.63 


4.40 ±0.53 


-2.63 ±0.61 


SMC 


TV flA'7A"7jC;A 

MU/U7oU 


19.74 


1.49 


2.95 


-13.65 


13.66 


-1.96 ±0.50 


2.20 ±0.59 


2.94 ±0.55 


-0.17 ±0.55 


marginal 




19 75 


1 55 


2 74 


15.58 


-3.56 


-0.61 ±0.54 


-2.67 ±0.47 


-0.00 ± 0.53 


-2.74 ±0.48 


m tir'fTi ti n 1 
llldiglilill 


M051847 


19.77 


1.96 


3.55 


60.31 


-6.46 


3.31 ±0.99 


1.27 ±0.64 


3. 16 ±0.99 


1.62 ±0.65 


marginal 


M045617 


19.78 


1.50 


5.45 


-23.80 


-13.26 


-5.33 ±0.82 


1.13±0.86 


4.11 ±0.83 


-3.58 ±0.85 


SMC 


M074956 


19.78 


1.74 


3.07 


-54.67 


18.42 


-1.41 ±0.96 


-2.73 ±1.38 


0.46 ±1.01 


3.04 ±1.34 


marginal 


M019505 


19.80 


1.63 


3.12 


38.84 


-46.89 


-0.55 ±0.63 


3.07 ±0.78 


-2.72 ±0.72 


1.53 ±0.70 


marginal 


M052255 


19.80 


1.69 


3.05 


9.34 


-6.03 


2.29 ±0.62 


2.01 ±0.52 


0.84 ±0.59 


2.93 ±0.55 


marginal 


M091341 


19.81 


1.70 


2.72 


-50.32 


39.02 


-1.80 ±0.77 


2.04 ±0.93 


2.67 ±0.83 


-0.51 ±0.87 


marginal 


M064031 


19.82 


1.96 


4.03 


-2.89 


6.28 


-0.31 ±0.49 


4.02 ±0.96 


3.78 ±0.89 


-1.40 ±0.60 


certain 


M006034 


19.84 


1.57 


3.54 


-43.05 


-70.74 


-2.23 ±0.87 


2.75 ±1.06 


-1.19±1.01 


-3.33 ±0.93 


marginal 


M040685 


19.85 


1.64 


2.72 


57.48 


-18.91 


1.49 ±0.67 


2.28 ±0.70 


0.70 ±0.67 


2.63 ±0.69 


marginal 


M070501 


19.91 


1.80 


2.64 


1.00 


13.37 


0.19±0.58 


2.63 ±0.84 


2.64 ±0.84 


0.00 ±0.58 


marginal 


M063762 


19.96 


1.56 


2.83 


5.52 


5.99 


-1.41 ±0.34 


-2.45 ±0.41 


-2.76 ±0.38 


-0.62 ±0.37 


marginal 



iMevlan. Dubath. & Mavod J1991I) that prompted our current 
study. Unfortunately, these two likely cluster members are too 
bright (V < 14) to have been included in our proper-motion 
sample. However, it is worth noting that their line-of-sight 
velocities are |v,| = 32.4 and 36.7 km s"' relative to the clus- 
ter mean, which would correspond here to about 1.7 and 1.9 
mas yr"' (for D = 4 kpc) and would not even be considered as 
high velocities by our definition. Meylan et al. deemed these 
stars potentially interesting based in part on a lower estimate 
of the central escape velocity than we have now, and in part 
on the fact that they appeared in a much smaller sample of 
velocities. Now, however, it seems that they may simply be 
part of the normal population of the v elocity distribution at 
V 3(T() (see also Gebhardt et alJl995l) . 

It must be emphasized that there are an additional 1,392 
stars in Table [S] with measured proper motions that have 
probabilities less than 0.001 in at least one component. As 



we have already discussed, although many of these measure- 
ments are effectively noise, some could be associated with 
perfectly reliable data which sample motions that are nonlin- 
ear due to strong accelerations. Any such stars should prob- 
ably be considered alongside those in Table |S] in any future 
attempts to understand the extreme tail of the 47 Tuc velocity 
distribution in detail. 

But even ignoring any nonlinear motions, it is still signifi- 
cant that the fraction of stars with anomalously high straight- 
line proper motions is so small: the number is at most 
41 /1 2, 974 ~ 0.3% within < 5 core radii of the cluster cen- 
ter, and it could easily be lower than this by a factor of 4 or 
more. Together with the scattered distribution of these stars 
in the (/x«,/ie) plane and the highly regular appearance of 
the velocity distribution of all other stars below Vmax — 2.6 
mas yr~', this should serve as a useful constraint on the- 
oretical analyses of any processes that could produce very 
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fast-moving stars in dense globular cluster cores (such as 
the mere presence of a central black hole, or close encoun- 
ters between various combinations of single stars, binaries, 
stellar- or intermediat e-mass black holes, and black-hole bina- 
ries; see, e. g., Davies & Benz 1995; Drukier & Bailvn 2003 
iManelli eTa l. 2005). In this context, note that for a Kinp 
Jl966h model with Wo = 8.6, roughly 75-80% of all stars pro- 
jected to within R < 5ro on the plane of the sky are in fact 
expected to lie in a sphere of radius 5ro around the cluster 
center. 

6. VELOCITY DISPERSIONS 

Having seen that there is nothing particularly untoward in 
the full proper-motion velocity distribution at the center of 
47 Tuc, the last analyses we undertake involve examining its 
second moment (velocity dispersion) in particular, as a func- 
tion of stellar magnitude and clustercentric radius in various 
subsamples of our data. We do not attempt to model the ef- 
fects of the cluster rotation in any detail — even though in prin- 
ciple this affects both the radial and the azimuthal compo- 
nents of proper-motion dispersion — ^both because the patchy 
azimuthal coverage of our velocity sample (e.g.. Figure O 
makes it difficult to constrain this well using only our data, 
and because the rotat ion in the inner few co re radii is known 
to be small anyway ( Mev lan & Mav or 1986). 

In ^6.11 we consider the v eloc ity dispersion of blue strag- 
glers in our sample. Section l6^ next looks at the dispersion 
of normal main-sequence stars as a function of magnitude 
and radius, and obtain s qua ntitative measures of the veloc- 
ity anisotropy. Section l63l compares the proper-motion and 
line-of-sight dispersions for turn-off mass stars to derive a 
kinematic distance to 47 Tuc. Having finally estimated the 
line-of-sight velocity dispersion in the cluster core, we draw 
on it in ^6. 41 to argue that the presence of ordered stellar mo- 
tions (rotation) will affect our estimates of the kinematics on 
these small radial scales at the level of ^ 2% or even less — a 
source of error that we do not correct for. Last, ^6.5l fits the 

f oper-motion a(R) profile with very simple modifications of 
in2 ( 1966) models to allow for the possibility of a central 
point mass. 

Throughout this section, we report all velocity dispersions 
after correction for both velocity-measurement errors and lo- 
cal transformation errors, as calculated from equation (I16> . 
For convenience, however, we drop the subscript "true" on 
cr from now on. In addition, we prevent the highest-velocity 
stars of 85.3l from inflating the estimated dispersions by only 
using stars with 

( S.Omasyr-i , V < 19 

/i,ot < < 3.2masyr-i , 19 < V < 19.5 (21) 
[ 3.6masyr-i , 19.5 < V < 20 

in the calculations — as we did earlier, in ^5.21 All uncertain- 
ties are estimated through nu meri cal bootstrap experiments, 
which were also described in ^5.21 

6.1. Blue Stragglers 

Figure [TO] above clearly shows the presence of a number 
of blue stragglers in our "good" proper-motion sample, and 
it is worth noting that non e of these appear as high- velocity 
stars (cf. Figure ml . In fact. 'K ing & AndersonH 2001') used an 
earlier, more restricted version of the current data set to show 
qualitatively that the blue stragglers in 47 Tuc have a lower 
velocity dispersion than the giants. We can now quantify this 
result. 



For our purposes here, we define blue stragglers as stars 
with 15.25 < V < 17.25 (recall that the main-sequence turn- 
off is at y = 17.65; Tabled and (U -V) < 0.7. Throughout 
our entire field, there are 23 proper-motion stars which satisfy 
these criteria, and 18 of these have R < 16'.'5. By compari- 
son, there are 441 proper-motion stars in the same magnitude 
range but on the red-giant branch of the cluster [which we 
take to be 1.5 < (t/-y) < 2.5], and of these 191 lie within 
R < 20". Although our velocity sample does not provide an 
unbiased sampling of the true radial distribution of any stellar 
population, this already indicates that the blue stragglers are 
more centrally concentrated than normal cluster members of 
the same magnitude. 

Figure |22] and Table |9] compare in more detail the relative 
spatial distributions and kinematics of the blue-straggler and 
red-giant branch stars in the inner R < 20" of 47 Tuc. The up- 
per right-hand panel of Figure |22] illustrates the CMD-based 
definitions of these two subsamples, as explained just above 
and also written in Table|5] The left-hand panels of the figure 
then show their one-dimensional distributions of /ig, fiQ, and 
speed /itot = (Pr + A'e)'''^^ very similarly to Figure [TTl above. 
The curves in each panel are Gaussians with dispersions equal 
to the measurement- and transformation-error corrected val- 
ues given by equation ( I16> , which are listed in Table|9j for the 
/itot distribution, we use the average value cr^ = (cr| + CTq) /2 
in equation (0. The lower right-hand panel then shows the 
cumulative radial distributions, out to R = 20", of the blue 




M/mas yr-1 R/arcsec 

Fig. 22. — Kinematics of blue stragglers vs. giant-branch stars of the same 
magnitude, in the inner R < 20" (roughly one core radius) of 47 Tuc. Only 
stars plotted as open stars (the blue stragglers) or large filled circles (red gi- 
ants) are used to construct the velocity and speed distributions in the left-hand 
panels; see Tablelglfor the definition of the two subsamples. In these panels, 
the open histograms pertain to the blue stragglers, and the bold solid curves 
are Gaussians with error-corrected velocity dispersions given in Table|9l The 
shaded histograms are for the red-giant stars, and the bold dashed curves are 
Gaussians with intrinsic dispersions also given in Tablel9l The lower ct^ for 
the blue stragglers is consistent with them being on average twice as mas- 
sive as main-sequence turn-off stars. The cumulative radial distributions of 
the blue stragglers, red giants, and all stars in our proper-motion sample (the 
curves from top to bottom in the lower right-hand panel) support this qualita- 
tively, although quantitatively the spatial distribution of stars in the velocity 
sample does not accurately reflect the true density profile. 



28 



D. E. McLaughlin et al. 



TABLE 9 

Kinematics of Stars with 15.25 <V < 17.25 and 
R < 20" 



Blue Stragglers Red Giant Branch 
0<({/-y)<0.7 1.5 < (f;-V)< 2.5 



AT 




18 


190 


OR (mas yr" 
OQ (mas yr" 


') 


367+"' "" 
-9.<M 


0.621+<'"'" 

U.OJl 


iTiot (mas yr" 


') 


AT.n+o'Sia 


(/3> 




fj C97+0.I86 
"•■^^'-0.427 


-0.046:«;i?? 



Stragglers (upper curve), the red-giant branch stars (middle 
curve), and all stars in the full sample of good proper motions 
(lower curve). The differences between the three are clearly 
significant to a high level of confidence (although they are also 
strongly affected by complex selection biases in these velocity 
samples). 

The interesting point here is the ratio of the blue-straggler 
and red-giant velocity dispersions: (t^(BS)/(t^(RGB) = 
0.69 (+0.07, -0.14), which is, to well within the errors, in- 
distinguishable from 1 / V2. This is just what we would ex- 
pect if blue stragglers are on average about twice as mas- 
sive as the normal main-sequence turn-off stars in 47 Tuc 
and the two populations are in energy equipartition. This is 
the first evidence of its kind for a relatively high mass for 
blue stragglers, and it is consistent with the idea that these pe- 
culiar stars represent some combination of tight (contact) or 
even coalesced equal-mass binaries, and the remnants of col- 
lisional mergers between main-sequence turn-off stars (see, 
e.g.,Eavies, Piotto, & De Anaeli 2004; MapelH et al. 2004). 

A typical mass ratio of 2 would also imply that thermal- 
ized blue stragglers should be spatially distributed with a core 
radius roughly 2""^^ = 0.67 times that of the turn-off and 
giant-branch stars (see eq. ||6l). For 47 Tuc, then, ro(BS) ~ 
0.67 X 20"84 w 14", which is indeed co nsistent with an ob- 
served ro(BS) = 10"-18", as estimated bv lGuhathakurta et al.1 
©93). 

As well as they mesh with each other, these argu- 
ments based on kinematics and spatial distributions are 
still indirect. Spectroscopic mass measure ments have been 
made for five blue stragglers in 47 Tuc by iDe Marco et al.l 
( 12QQ51 m = 0.29-1. 88M fn), and for one other in the clus- 
ter bv lShara. Saffer. & Liv io ( 1997. m =1.1 ± O.4M0). The 
weighted average mass of these 6 stars ( usin g weights based 
on the relative errors as described by iDe Marco et al.) is 
(m) ~ 1 .4 ± Q.5Mq, which is to be compared with the main- 
sequence turn-off mass of about O.85M0 (Table [7). While 
perhaps roughly consistent with the indirect suggestions that 
m(BS) = 2 X m(TO) on average, the number of direct mass 
measurements is still very small, and the spread and uncer- 
tainties in the individual values are quite large. More work is 
needed in this area. 

The last line of Table |9] gives rough estimates of the ve- 
locity anisotropy {(3) obtained by comparing the radial and 
azimuthal proper-motion dispersions of both the blue strag- 
glers and the red giants (see eq. 1231 below). While the giant- 
branch stars are clearly consistent with the expected isotropy 
((/3) = 0), it is perhaps noteworthy that the formal > for 
the blue stragglers implies an apparent radial bias in their or- 
bits. However, given the small number of stars, this result is 



not significant at even the one-sigma level. 

6.2. Dispersion as a Function of Stellar Magnitude and 
Projected Radius 

In order to define the run of proper-motion dispersion vs. 
clustercentric radius, we first break our good velocity sam- 
ple of 1 2,974 stars into four broad magnitude/mass bins as 
in V < 18.5, corresponding to 0.77 < < 0.9Mq; 

18.5 < y < 19, meaning 0.73 < < 0.77 Mg; 19 < V < 
19.5, or 0.69 < < O.73M0; and 19.5 < V < 20, for 
0.65 < < 0.69Mq. We also apply the color selection spec- 
ified in equation ( I18> . to exclude blue stragglers in particular 
from the analysis. In each of these ranges we use stars sur- 
viving the velocity cuts in equation ( I2U to compute the mean 
velocities in the RA and Dec directions, and the correspond- 
ing error-corrected dispersions, in a series of narrow, concen- 
tric annuli. From these we confirm that (fia) = {l^s) = to 
within the uncertainties at all radii in general — as expected 
from the relative nature of our proper motions — and then set 
= (/^e) = to compute the true radial and azimuthal dis- 
persions (Tr and (Te, and the average = [(u^ + ctq)/!]^/^, 
in the same annuli (with A/lot - 1 replaced by A/Jot in eq- 1161 ). 
Finally, we find the relative difference (cr^ - CTq) / cr^. 

Table [Tol presents the results of all these calculations per- 
formed in a series of narrow annuli in each magnitude bins. 
(Table [Tolcan be found at the end of this preprint). Note that 
we increase the annulus size as the stellar magnitudes grow 
fainter, in order to avoid drastic declines in the numbers of 
stars per annulus. As a result, in the bins with V > 18.5, we 
actually use two sets of overlapping annuli to aid in evaluat- 
ing the extent to which the choice of bin positioning affects 
the measurements. All of the data required to repeat this ex- 
ercise for any other choice of magnitude or radial binning (or 
any alternate definition of the velocity sample) are contained 
in Table |5l which is available in full in the online edition of 
the Astro phy sical Journal. 

Fi sure I23I shows cts, (jq, and (o's-cr||)/cr^j as functions of 
R from TablellOl Note that we have only plotted every second 
number from the table for the fainter magnitudes V > 18.5, 
so that every point shown is statistically independent of the 
others. For clarity, we have attached errorbars only to the data 
in the first and third magnitude bins. 

In each of the top two panels of Figure |23l we draw 
the predicted one-dimensional, projected velocity-dispersion 
profile for a King ( 1 966) model with Wo = 8.6 and ro = 
20"84 (Figure [151 in E)- and with cro = 0.633 mas yr"' = 
12.0(D/4kpc) kms"' as in equation ^ of EJl These 
curv es w ere calculated as described in Appendix [E] (see 
eqs. IIB8I and IB9I ). and then smoothed by averaging cr^aA^R) 
over a 3"-wide annulus at each point, to match the binning 
of the brightest data. In this specific model, the observable 
central dispersion in any component of velocity is a(R = 0) = 
0.982 0-0 = 0.622 mas yr"', or 1 1.8 km s"' forD = 4kpc. Note 
that the model predicts a decrease of about 15% in erg = o-q 
from the center Xo R= 100". Although this amounts to just 
under 0.1 mas yr~', it is nicely resolved by the measured dis- 
persions for the brighter (more numerous and less uncertain) 
stars. The clear tendency for a decreasing cr as a function of R 
suggests — ^rather more strongly than do the weakly populated 
tails of the full, two-dimensional velocity distribution at any 
radius — that a King-model distribution function (with a finite 
escape velocity that decreases monotonically towards larger 
radii) is indeed a closer approximation to the cluster dynam- 
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TABLE 11 

Aperture Dispersions and Average Anisotropy 



Stellar Magnitude 




[mas yr"' ] 


[mas yr"' ] 


[mas yr"' ] 




(/3) 


V < 18.5 


6336 


0.589 ±0.006 


0.577 ±0.006 


0.583 ±0.004 


0.04 ±0.03 


0.06 ±0.04 


18.5 < y < 19 


2327 


0.587 ±0.010 


0.565 ±0.011 


0.576 ±0.008 


0.08 ±0.05 


0.11 ±0.07 


19 < y < 19.5 


2421 


0.591 ±0.013 


0.586 ±0.013 


0.589 ±0.009 


0.02 ±0.06 


0.03 ±0.09 


19.5 < y < 20 


1693 


0.479 ±0.026 


0.495 ± 0.024 


0.487 ±0.018 


-0.07 ±0.14 


-0 11+0-20 



ics than is a truly isothermal sphere with a perfectly Gaussian 
velocity distribution. 

Also in the top two panels of Fig. |23] note the slightly 
higher dispersion of the fainter stars inside R < 20", which has 
already been discussed in connection with Figure [19] above. 
Apparent too are the low values and large scatter in the error- 
corrected (jR and (Te at larger radii, for stars in the range 
19.5 < y < 20. As was also discussed above, this is due to 
the generally large velocity uncertainties at these magnitudes 
in these areas of the field. Our faint stars might provide useful 
kinematics constraints on better dynamical models for 47 Tuc 
if the definition of the velocity sample is revised to exclude 
stars with uncertainties close to or larger than the true velocity 
dispersion; but, as Figure |9l makes clear, this will drastically 
decrease the sample size. 

The bottom panel of Figure |23] shows the relative differ- 
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Fig. 23. — Proper-motion dispersions in the projected radial and azimuthal 
directions relative to the cluster center, and their relative difference, as func- 
tions of radius in four magnitude bins. Curves in the upper two panels are 
one-dimensional, projected cf(R) profiles for a Wq = 8.6 Kina >1966 ) model 
with ro = 20" 84 and (Tq = (1 + 1/45)'/^ X 0.633 mas yr"' (see text). Note 
that the stellar kinematics at y > 19.5 and outside of R> 20" in particular 
are error-dominated. The relative difference o f var iances in the bottom panel 
is related to the velocity anisotropy; see TablelTTIfor average values over the 
entire field. 



naturally the most uncertain kinematical parameter that we 
calculate, and its noisy profile reflects this. The fact that it 
rarely differs significantly from 0, and the absence of any sus- 
tained trend in clustercentric radius, suggest that cr| = ctq on 
the whole. Again, then, in keeping with the earlier modeling 
of Mevlan ( 1988, 1989) and our discussion in ^ the velocity 
distribution in the inner several core radii of 47 Tuc appears 
to be essentially isotropic. 

This impression can be better quantified by considering 
how the observable velocity dispersions in the R and Q di- 
rections and along the line of sight, z, are related to the intrin- 
sic, unprojected ar, (Jq, and ct^. For a spherically symmet- 
ric star cl uster with an ellipsoidal v eloc ity-dispersion tensor 
(ae = o-aI. Leonard & MerrittHl989h and lGenzel etanj2000h 
show that 

<7^(r)cos^e + aj(r) sin^ 6 
a^ir) sin^6l + CT^(r) cos^9 



(22) 



where r is the unprojected clustercentric radius and 9 is the an- 
gle between the unit vectors r and z. With velocity anisotropy 
parametrized by the usual function (3(r) = l-aj/a^ (e.g., 
iBinnev & Tremaine 1987), averaging the last two of equa- 
tions ( I22t over all 9 and all r leads to jLeonard & MerritJ 
fT989.) 



1- 



3(^i)-3(c7c: 
3(<7|)-(a| 



(23) 



Thus, for an isotropic velocity ellipsoid, (/3) = and we have 
{'^r) - {^%)- For purely radial orbits, (/3) = 1 and ((t|,) = 0, 
the same as the unprojected ag. For purely circular orbits, 
the unprojected (jy = 0, so {(3) = -co and the averaged proper- 
motion dispersions satisfy ((Tq) = 3(0-1). 

Equation ( I23t is strictly only valid as an estimator of the 
globally averaged anisotropy, and stars spread over an entire 
cluster should be used to compute it; it is not a particularly 
meaningful quantity in narrow annuli such as those consid- 
ered in Table [Tol fo r exa mple. However, as was also men- 
tioned at the end of ^5.31 the majority of stars with projected 
radius /? < 5 ro in a Kina -model cluster similar to 47 Tuc do 
actually lie in a sphere of the same radius. Thus, calculat- 
ing a single {(3) for all the stars in our proper-motion sample 
will still yield a useful approximation to the volume-averaged 
{(3) over the inner 100" ~ 4.8 core radii of this cluster (with 
a small contamination from stars at larger distances along the 
line of sight). 

Table [Tn reports the R and 8 aperture dispersions for our 
entire field, as well as their relative difference and the corre- 
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spending averages (aj^)^^^ and (/3), with the sample broken 
into the same magnitude bins as before. Note that the dis- 
persions in the fifth column of the table are just those used 
in Figure ^] above and that, as was also discussed in ^5.21 
the results at the faintest magnitudes are dominated by mea- 
surement errors and are thus unreliable. Apart from this, 
there is an apparent preference for a slight radial anisotropy 
({<^i) > {^h) and {13) > 0) in the average. Although the de- 
parture from isotropy is not highly significant statistically, the 
sense of it is just what we expect for a dynamical structure in 
which the stellar orbits are isotropic at the center of the cluster 
and become gradually more radially biased toward very large 
radii (see iMevlannlQS S). 

6.3. A Kinematic Distance 

Section 13.21 presented details of a set of radial velocities 
we have obtained for stars in the inner R < 105" of 47 Tuc. 
Here we compare the velocity dispersion of these stars to that 
of our proper-motion stars to estimate a distance to the clus- 
ter Recall that we have useful radial velocities only for 419 
bright giants with magnitudes 11 < V < 14, (corresponding 
to masses ~ Q.9Mq), whereas all of our proper-motion 
stars are fainter than this and can have masses as low as 
OT* ~ O.65M0. In order not to bias the derived distance by 
comparing the kinematics of stars with too widely different 
masses, we work only with rather bright proper-motion stars: 
V < 18.5, and thus O.TSM© < < 0.9Mq according to the 
formula in Table Applying the color cuts of equation ( I18> 
to exclude blue stragglers, and the velocity cut in equation 
illi . we are left with 6,336 proper-motion stars having an av- 
erage magnitude (V) = 17.74 and hence an average mass of 
O.84M0. 

Since the stars in the core of 47 Tuc are expected to be in 
energy equipartition, we might naturally worry that a proper- 
motion sample with V < 18.5 could have a dispersion some 
(0.9/0.84)'/^ = 1.035 times larger than a sample of stars with 
masses strictly equal to those in our radial-velocity sample — 
in which case, the true distance to the cluster might be longer 
than what we derive in this Section, by about 3.5% (roughly 
0.15 kpc). While this is a valid concern in principle, in prac- 
tice we have found that carrying out the analysis below using 
a smaller proper- motion sample of 1,895 stars with V < 17.5 
(for a mean mass of O.88M0) formally changes our inferred 
distance in the opposite way to what we would expect, but 
by a statistically insignificant amount: D{V < 17 .5) / D{V < 
18.5) ~ 0.99 ±0.05. Given our inability to resolve directly 
the signature of energy equipartition in the kinematics of the 
brightest proper-motion stars, we simply present the results of 
our work with the larger, V < 18.5 sample and do not attempt 
to compensate for any anticipated effect. 

Another potential complication is that in using these data to 
estimate the distance to 47 Tuc, we assume that the cluster is 
not rotating. As we have ment ioned alre ady, it is kno wn that 
47 Tuc does rotate (|Mevlan & Mavodl98 6: Anderson & Kind 
|2903a'), but we are saved from having to apply much more so- 
phisticated dynamical modeling (with more free parameters, 
constrained by additional data with further sources of uncer- 
tainty) by the fact that the level of rotation in the central-most 
regions of the cluster presently under consideration is low rel- 
ative to the random motions. In ^6.41 we give a brief but more 
quantitative a posteriori justification for our neglect of this is- 
sue. 

Given equations ( I22t and the definition (3{r) = 1 -a^/a^, 
it is straightforward to show that the velocity dispersions pro- 



jected along the line of sight and into the radial and azimuthal 
directions on the plane of the sky are connected by 



2-1 
1-/3 



(24) 



at any radius r in a spherical and non-rotating cluster, for all 
three dispersions measured in the sa me un its. If we then per- 
form the same spatial averaging as in ^6.21 substitute equation 
( I23> for j3, and use equation M5\ to convert between proper 
motions measured in mas yr"' and radial velocities in km s"', 
the distance follows as 



D 
kpc 




t2\1/2 



4.74(^2)1/2 



(25) 



For non-zero anisotropy, this formula is really applicable 
only to a global sample of stars covering the entire cluster, 
although — as with the expression above for {(3) — its use for 
our R <5r(, field is still justified. And in any case, if /3 = it 
is more generally valid for any local sample; (ct^ ) is then sim- 
ply an improved estimate of a single proper-motion variance 
given two independent estimates from the R and 8 compo- 
nents. Thus, we have estimated the distance to 47 Tuc in two 
different (though interdependent) ways, both of which are il- 
lustrated in Figure |24] 

The top panel of Figure|23plots cr^ against R from the first 
section of Ta bleFTOI The curves running through the points are 
models with two different core radii, which we 
will discuss below. The middle panel of the figure then shows 
the error-corrected line-of-sight dispersion CT; as a function of 
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Fig. 24. — Comparison of proper- moti on dispersion profile crp(S) for stars 
with V < 18.5 (first section of Table lTol against the line-of-sight dispersion 
profile cr^(R) for stars with V < 14, and the resulting kinematic distance es- 
timate in a number of circular annuli. Only every second point in the lower 
two pane ls is statistically independent, as they refer to overlapping annuli; 
see Table ll2i Curves in the upper two panels are Kins 1 1966) model profiles 
for Wo = 8.6. with two different core radii rg and two different ctq parameters; 
see text. The bold horizontal line in the bottom panel is at the median of 
the independent annulus-by-annulus distance estimates, and the lighter lines 
denote its 68% confidence interval (see eq. 1261 ). Taking the ratio of either 
pair of m odel curves in the upper panels yields essentially the same distance 
(eq. |271). 
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TABLE 12 

LiNE-OF-SlGHT VS . PROPER-MOTION KINEMATICS 



Radii 
[arcsec] 
(1) 



LOS: 11 < y < 14 

"TV 

[kms-'] 
(2) (3) 



PM: V < 18.5 
"TV 5> 

[mas yr"'] 
(4) (5) 



D 
[kpc] 
(6) 



0.0-15.0 
7.5-22.5 
15.0-30.0 
22.5-37.5 
30.0-45.0 
37.5-52.5 
45.0-60.0 
52.5-67.5 
60.0-75.0 
67.5-82.5 
75.0-90.0 
82.5-97.5 
90.0-105.0 



35 
76 
104 
96 
79 
77 
68 
53 
54 
49 
37 
37 
42 



9 76+'^* 
11.82+'05 

12 01+*-*'' 
1 2 Q4+0.89 



12.66+' 



5+0:91 

10 AS+O Sfi 

9 63+''-^8 
-0.94 

10.84!»;g 

11 01+""' 
ll.ui , 03 

10.87:!?' 

8.79:[g 

9.48:?;™ 



1283 
1314 
569 
498 
824 
1086 
1279 
1349 
1260 
1125 
952 
550 
169 



0.608:»;«i« 
0.604:0:«0' 
600+" ''" 
0.624ti 
0.608^;«|1 

A <:o<:+0.010 

n ';sn+" 008 
0.5SO_oQm 

570+0 "'''' 
-QQQ9 

<;67+fl-010 

U.JO/_Q0QC, 

546-ni "!" 
A 54o5fl.0t4 
0.563^^5 



a ^q+0.56 
-^■^ -0.73 

4.13 



4.22; 
4.37 



m 

-0.34 



4.39+ 



0+0.32 
'-0.39 
3 76+"-32 

-0.32 

3 50+" " 

3.58:"3' 

4 n4+"-^^ 

4 10+"3't 
4.20+"-1" 



3.44; 
3.55: 



radius (and'King'-model profile curves again), for stars binned 
in 15"-wide annuli. Unlike the proper motions, we have plot- 
ted these results for both of two slightly offset, overlapping 
sets of annuli, so that only every second a, point is statistically 
independent. Tablell2lgives the radii of the annuli defined for 
the radial velocities, followed by the error-corrected and 
the true proper-motion cr^ in each, and then the distance in- 
ferred by assuming velocity isotropy to ap ply e quation ( I25t in 
each annulus. The bottom panel of Figure |24]plots these dis- 
tance estimates, again for all of the overlapping annuli. Only 
seven of them are independent, and their median value is 

£) = 4.04 (+0.1 8, -0.49) kpc , (26) 

as indicated by the horizontal lines in the graph. The uncer- 
tainties here delimit the 68% confidence interval on the me- 
dian for 1000 artificial data sets generated by bootstrapping 
with replacement from the measured distances. 

Clearly, it is the relatively small number of radial veloci- 
ties at our disposal which limits the precision of any distance 
estimate here, and the noisiness of the profile a,{R) which 
limits its accuracy in any single radial bin. It would thus be 
preferable to use all of the velocities over our entire field at 
once to obtain a single D. This is nontrivial, however, since 
neither component of the dispersion is spatially constant and 
the radial-velocity and proper-motion stars have significantly 
different (and highly irregular) distributions on the sky; see 
Figure ^5 above. The total cr^ and CT; for our full samples 
therefore represent differently weighted averages of a single 
underlying velocity-dispersion profile, and if we were to put 
them directly into equation ( I25t a biased distance could re- 
sult in principle. It is possible, thoug h, first to cor rect for the 
different spatial samplings using our lKing I model de- 

scription of 47 Tuc. 

The bold, dash-dot curve in the top panel of Figure|23is the 
projected and smoothed one-dimensional velocity-dis persi on 
profile plotted in Figure |23l As was also discussed in ^5.21 to 
normalize this model to the data we first computed the dimen- 
sionless profile (tI^o^/ctq vs. (R/ro), as described in Appendix 
151 for Wo = 8.6. Then, given any fixed value of the King core 
radius ro, the average model dispersion for a set of stars with 
measured clustercentric radii {Rj} is given by equation ( I19> . 



To repeat: 

(^mod) ^ 1 ^mod(^//^0) 

Withro = 20"84 from our fit to the surface -brightness profile 
in ^.11 the measured positions of the 6,336 proper-motion 
stars being used here yield (cr,^jjjj)/crQ = 0.849±0.015. Equat- 
ing (cr^d)'''^ to the true (cr^)i/2 = 0.583 ± 0.004 mas yr"' 
from Table [TT] then impHes ctq = 0.633 ± 0.010 mas yr , 
just as in equation ( I20> above. Turning now to the radial- 
velocity sample, the 419 observed stellar positions correspond 
to (cr,nod)/'''o - 0-851 ±0.059 (only fortuitously close to the 
result for the proper-motion stars) if ro = 20."84 still, while the 
observed {a}) '/^ = 1 1 . 1 ± 0.4 km s"' after correction for mea- 
surement errors. Thus, cto = 12.05 ±0.85 km s"'. The model 
cr,(/?) profile for Wo = 8.6 and ro = 20."84, normalized by this 
value of (To and smoothed with an annular filter of width 15" 
to mimic the binning of the data, is drawn as the heavy dash- 
dot curve in the middle panel of Figure |24] An estimate of 
the distance to 47 Tuc, which effectively smooths over the 
fluctuations of a, and in narrow annuli and does not re- 
quire an assumption of strict velocity isotropy, then follows 
from using the proper-motion and radial- velocity <to values in 
equation M5V . 

D = 4.02 ±0.35 kpc , (27) 

in good agreement with our first result. Note that the ~9% 
relative uncertainty in D is twice the amount caused by num- 
ber statistics alone, due to the necessity of correcting for the 
different spatial distributions of the two velocity samples. 

As we discussed in connection with Figure l20l above, and 
can also be seen from the top panel of Figure|24] the observed 
(Tf^(R) profile in the inner w 40" (~ 2 core radii) appears rather 
flatter than expected for a Wo = 8.6 King (1966) model with 
the core radius appropriate to stars brighter than the main- 
sequence turn-off. This point will also appear again in ^6.51 
Thus, we have repeated the exercise just described using a 
W) = 8.6 King model with ro = 27 ."5 rather than 20."84. In this 
case, we find do = 0.620 ±0.010 mas yr"' from the proper- 
motion sample, and ao = 11.8 ±0.8 km s"' from the radial 
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velocities. The dotted curves in the top two panels of Figure 
I24I show the model profiles with this larger ro and the lower 
(Jo values, again smoothed by amounts matching the binnings 
of the data. Taking the ratio of the normalizations gives a 
distance identical to that in equation ( I27> . demonstrating that 
the result is insensitive to fine details of the model used to 
smooth over the spatial distributions. 

Recent determinations of the distance to 47 Tuc by standard 
CMD fitting range from D = 4 .45 ±0.15 kp c (Percival et al. 
I2002h to D = 4.85 ± 0. 18 kpc ( Gratton et al. 2003), while a fit 
to the white-dwarf cooling sequence by Zoccali et al. (2001) 
returned D = 4.15 ±0.27 kpc. Our result clearly supports 
the shorter white-dwarf value (which in turn implies a cluster 
age of 13 Gyr; see Zoccali et al.), although the lowest of the 
CMD-based distance estimates may still be consistent within 
the uncertainties. 

6.4. Neglecting Rotation 

It would be surprising if more sophisticated modeling did 
not uncover slight biases in our distance estimate due to our 
simplifying assumptions of spherical symmetry and zero ro- 
tation. But it seems very doubtful that the value of D could 
change by more than the 9% margin of error already al- 
lowed by the formal uncertainties. In fact, we expect that 
any rotation-related biases in the kinematics we have derived 
will be present only at the level of ^2%, which is comparable 
to the minimal, 1 / statistical errorbars that are unavoid- 
able even in the brightest and lowest-uncertainty subsamples 
of proper-motion stars. 

The signature of rotation in the radial velocities appears as 
a sinusoidal dependence of on the projected azimuth Q, 
and simply fitting a sine curve to the 419 reliable line-of-sight 
velocities available at R < 105" yields a rotation amplitude 
Vrot — 3 km s~'. The azimuthally averaged velocity variance 
due to rotation alone is therefore {V^f sin^ Q) = 4.5 (km s~')^. 
Meanwhile, the observed, error-corrected variance for all 419 
radial-velocity stars together is (ct?) = (11.1 kms~')^ (see 
^6.3> . The velocity dispersion over our field that is due to 
truly random stellar motions, rather than rotation, is then 
(1 1.12-4.5)'/^ = 10.9 km s"'— just over 98% of the observed 
value. 

The velocity distribution on these spatial scales in 47 Tuc is 
essentially isotropic ( ^6.2> . and there is ev idence that the in- 
clination angle of the rotation axis is ^45° J Anderson & Kind 
|2()03a), so any signature of rotation in our proper-motion dis- 
persions can also be expected to appear at roughly ^2% lev- 
els, which is too small to discern accurately once all sources 
of statistical uncertainty (number statistics; measurement and 
transformation errors; irregular spatial distributions) are taken 
into account. In any event, our distance estimate particularly 
should be rather robust, since it relies on the ratio (J-/(J^. 

6.5. Does 47 Tuc Harbor a Central Black Hole? 

Our final analysis looks at the kinematics of the innermost 
^3" of 47 Tuc, the area in which the effects of a compact 
central mass concentration (if one is present) should be most 
clearly observable. 

Returning for a moment to the results of ^6.31 it is a sim- 
ple matter to calculate the observable velocity dispersion (for 
stars at the turn-off mass) at the exact center of 47 Tuc from 
our estimates of the King-model scale parameter, ctq. With 
Wo = 8.6, the projected one-dimensional UmodiR = 0)/cto is 
0.982. For the cto values associated with ro = 27 ."5 (arguably 



a better description of the inner one core radius of the proper- 
motion velocity-dispersion profile in Fig. l24> . we thus have 



<j^,(R = 0) = 0.609 ± 0.0 10 mas yr" 
a-(R = 0) = 11.6±0.8kms"' . 



(m, ~O.85M0) 
(28) 

Obviously this assumes that velocity dispersion is nearly con- 
stant as R approaches 0. On the other hand, the innermost 
point in the top panel of Figure |23 appears to fall above this 
extrapolation. Even though the discrepancy is hardly signifi- 
cant at even the one-sigma level, it is nevertheless of potential 
interest in light of recent claims for evidence of massive black 
holes at the centers of globular clusters. 

It is now well established that the masses of super- 
massive black holes at the centers of galaxies correlate 
tightly with the central velocity dispersions of the stel- 
lar bulges (Ferrai-ese & Merritt 2000; Gebhardt et al. 2000). 
Recently, claims have been made, on the basis of radial- 
velocity studies, th at the Galactic globular cluster M15 
(Gerssen et al. 2002) and the very massive cluster Gl in 
Andromeda (Gebhai'dt. Rich, & Ho 2002, 2005) may hai-- 
bor large central black holes with masses M, ^ 10^- 
10^ Mq. If so, such black holes could lie on a simple 
extension of the M,-<t relation for galaxies. For M15, 
the data show less than a 1-cr significan ce for a black 
hole dOerssen et a l. 2003; McNamara, Harris on'. & AndersonI 
2003), and theoretical m odels show there is no need to 
invoke one (^Baumaardt et al. 2003a). Given the data of 
Gebhardt, Rich, & Ho (2002), Baumaardt et al. (2003b) con- 
cluded there is no need for a black hole in Gl either, although 
Gebhardt. Rich. & Ho (2005) use more recent data to argue 
for an improved significance of detection. It is thus unclear 
whether globular clusters do in fact contain "intermediate- 
mass" central black holes, and it is worthwhile asking whether 
our proper motions in 47 Tuc can add anything to this discus- 
sion. 

The galactic M,-a relation given bv lTremaine et al. ( 200^ 
uses the bulge velocity dispersion averaged over an effective 
radius, which in a Wo = 8.6 Kins ( 1966) model is a(R < R^) = 
0.91 a(R = 0). Given equation (|28ll, then, cr(R < R,) ~ 10.6 
km s"' in 47 Tuc, and extrapolating the Tremaine et al. rela- 
tion gives a "predicted" black-hole mass of roughly IOOOM0. 
This corresponds to a sphere of influence GM,/a{Q)^ « 
0.032 pc = l."6(D/4kpc) and suggests a potential observa- 
tional signature in our proper-motion data. On the other 
hand, there is no a priori justification for extrapolating the 
empirical scaling for galaxies down to the globular cluster 
range, and the results of doing so are highly dependent on 
the adopted parametrization of the relation. Thus, for exam- 
ple, the steeper dependence of black-hole mass on (central 
galactic velocity dispersion advocated by Ferrarese (2004) 
leads to M, w 360Mq in 47 Tuc, which would have a much 
smaller effect on the observable kinematics. Conversely, the 
log-quadratic relation of Wvithe (2006) implies an untenable 
(though highly uncertain) M, -^2-3 x lO^^^M©. 

Aside from these rou gh expectations, it is also worth not- 
ing that iGrindlav et a l7 (2001) use the positions and disper- 
sion measures of millisecond pulsars in the core of 47 Tuc, 
together with a Chandra upper limit on the X-ray luminosity 
of a central source, to estimate an upper limit of M, < 47OM0 
for the mass of any central black hole. This corresponds to 
a sphere of influence, r < 0.015 pc ~ 0"8(Z)/4kpc), which 
would be very difficult to probe with the current kinemati- 
cal data. On the other hand, an X-ray upper mass limit is 
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fairly sensitive to the assumed efficiency and spectral tem- 
perature for radiation from gas accreting onto a black hole 
(M, ^ e'^I^T) — uncertainties that are, of course, avoided 
with kinematical mass estimates. 

6.5.1 . Simple Models for the Velocity Dispersion Profile 

T o address this issue, we have constructed a series of mod- 
ified King (1966) models allowing for the presence of a cen- 
tral point mass, and fitted these to the proper-motion dis- 
persion profile (Jfj^(R). We first find the dimensionless den- 
sity profile pK(r) / po and enclosed mass profile GMxCr) /cTgro, 
both as functions of F = r/ro, for a regular King model with 
given W(). Then we add a black hole with dimensionless 
mass GM^/a^rQ, assume that the stellar density profile is un- 
changed, and solve the isotropic Jeans equation. 



d 

d7 



1 Pk(T) 
r^ po 



GMk{T) ^ GM, 
alro 



crlro 



(29) 



for a new cr,^(?)/(TQ profile (by integrating inward from the 
original tidal radius, where pa^ = 0). The projected, one- 
dimensional velocity-dispersion profile to be compared with 
the observed (J^{R) follows fro m th e usual averaging of a} 
along the line of sight (e.g., eq. IB9I ). 

Note that these models do not have the density cusp that 
is usually associated intuitively with the presence of a central 
black hole. This is because they do not attempt to address the 
origin of the black hole or the dynamical evolution of the stel- 
lar distribution function in its presence. Rather, they simply 
describe the self-consistent kinematics of stars around a dark 
point mass ^iven that the density structure is well matched by 
a normal King model. Baumgardt, Maki no. & Hut (2005) use 
detailed N-body simulations to show that the visible stars in 
globular clusters with intermediate-mass central black holes 
can indeed be reasonably well described by'King'-type density 
profiles (the expected density cusp is essentially confined to 
heavy stellar remnants), and that the stellar kinematics are ac- 
curately reproduced by isotropic Jeans modeling. The great- 
est limitation to our simple approach is the assumption of a 
single stellar mass, since in principle neutron stars or heavy 
white dwarfs that are strongly concentrated to the center by 
mass segregation can affect the kinematics of observable stars 
in ways similar to a true point mass. Although we have al- 
ready argued t hat th e total mass of such remnants in 47 Tuc 
is quite small (" ^5.21 see the discussion around Figurel20t. ig- 
noring the issue altogether means that fitting the projected so- 
lutions of equation ( I29> to an observed (J^{R) profile is likely 
to overestimate the mas s of an y black hole in general. 

To fit the black-hole iKind models, we first calculated the 
afj_{R) profile again for our "good" proper-motion stars with 
V < 18.5 and restricted colors as in equation dlSt . but binned 
into l"-wide annuli in order t o res olve the small-radius be- 
havior better than we did in ^6.21 and ^6.31 (However, we 
did check that our results are insensitive to the binning.) 
We then solved equation \29\ for the dimensionless profile 
(Jr/<JQ VS. r/ro, for a number of models with fixed Wo in the 
range Wb = 8.6 ± 0.5 and dimensionless black-hole masses in 
the range < GMt/a^ro < 0.50. Every model profile was 
projected along the line of sight, smoothed with an annular 
filter of width AR = 1" to match our treatment of the data, 
and fit to the observed iTp(/?) profile by finding the normal- 
izations ro and cto that minimize the usual error-weighted 
statistic. We used only points with R < 90" when computing 
X^, since the number statistics at larger radii are very poor. 



TABLE 13 

Central Black-Hole Model Fits to V < 18.5 Kinematics ^ 



GM./alro 




O"0 


'■o 










[mas yr~' ] 


[arcsec] 


[Mq] 












i.6 






0.00 


98.47 


0.614 


26.34 





best fit 


0.01 


98.57 


0.612 


26.94 


164 




0.02 


98.75 


0.610 


27.56 


334 




0.03 


99.03 


0.608 


28.19 


509 




0.04 


99.40 


0.606 


28.85 


690 


< 68% C.I. 


0.05 


99.89 


0.605 


29.51 


877 




0.06 


100.49 


0.603 


30.17 


1070 




0.07 


101.21 


0.601 


30.86 


1270 




08 


102 04 


\J.3yy 


J i . JO 


1470 




0.09 


103.01 


0.597 


32.25 


1680 




0.10 


104.11 


0.595 


32.86 


1890 




0.1 1 


105.34 


0.593 


33.48 


2110 


> 99% C.I. 


0.12 


106.72 


0.592 


33.97 


2320 




0.13 


108.24 


0.590 


34.37 


2530 




0.14 


109.91 


0.588 


34.68 


2730 




0.15 


1 1 1 72 


587 


34 92 


2930 




0.16 


113.66 


0.586 


35.08 


3130 








Wo = i 


i.l 






0.00 


98.25 


0.615 


30.72 





best fit 


0.04 


99.49 


0.608 


33.44 


805 


> 68% C.I. 


0.08 


102.68 


0.600 


36.80 


1720 


> 95% C.I. 


0.10 


105.15 


0.596 


38.47 


2220 


> 99% C.I. 








).l 






0.00 


98.76 


0.613 


22.39 





best fit 


0.05 


99.77 


0.603 


25.29 


749 


68% C.I. 


0.09 


102.30 


0.597 


27.50 


1430 


< 95% C.I. 


0.12 


105.37 


0.592 


28.78 


1970 


99% C.I. 



^Fits are to the proper-motion velocity-dispersion profile of stars with R < 90", 
V < 18.5, and P(x") ^ 0.t)OI in both RA and Dec components of proper motion. 
Stars are binned in 1 -wide annuli, resulting in 87 points fitted. 



Dimensional black-hole masses assume a distance of D = 
any other distance, multiply by (D/4.0kpc)"^. 



4.0 kpc to 47 Tuc. For 



For any value of Wq, the best-fit black-hole mass is, of 
course, that which yields the lowest x^ in the model grid 
of GM,/alro values that we defined. The 68% (l-cr) con- 
fidence interval is defined by the values of GM,/alr() for 
which x^ < Xmin"*" 1' '^^c 95% confidence interval is set by 
^ Xmin + 3.84; and the 99% confidence interval, by x^ < 
xL + 6-63. 

The first part of Table II 31 presents the details of several of 
the fits to (7^{R) obtained for Wo = 8.6, which best matches 
the spatial structure of 47 Tuc for R < 1000" (see The 
columns in this table are the dimensionless black-hole mass 
specified for each fit, followed by the the fitted model 
scale velocity and radius (both allowed to vary freely in the 
fitting), the physical value of the black-hole mass, and com- 
ments identifying the formal best fit and approximate 68%, 
95%, and 99% confidence intervals. Note that the dimen- 
sional black hole mass is obtained as 



M. GM, 

-^=25.40^ 



D 
kpc 



CTO 



mas yr 



arcsec 



(30) 



which clearly is quite sensitive to the assumed distance to the 
cluster Naturally, we take D = 4.0 kpc, from ^6.31 

Assuming Wo = 8.6, the best fit to (j^{R) for this sample 
of stars with V < 18.5 is achieved with M, = 0, although 
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R/ arcsec 

Fig. 25. — Upper panel: Number-density profile of stars brighter 
than the main-sequence turn-off (Vto = 17.65, or F475W ~ 17.8), as de- 
rived from our master list in Table |4| and compared to similar results from 
[Howell. Guhathakurta. & Gilliland (2000). To obtain our numbers (filled cir- 
cles), stars were binned into l"-wide circular annuli at R < 20" and into 
3"-wide annuli at 20" <R< 90". Bold, dash-dot curve is a Wq = 8.6 KinJ 
model curve with ro = 20" 84; the dotted curve has a larger core radius of 
tq = 21'.' 5. Bottom panel: Proper-motion velocity-dispersion profile for stars 
with V < 18.5, binned in radius as for the density profile. Curves represent 
fits of Wo = 8.6 Kina models with central black holes of the masses indicated, 
which correspond to the formal best fit (lower, sohd curve) and the 68%, 95%, 
and 99% (uppermost dashed c urve) upper hmits. Normalizations rg and oq 
for each case are given in Table[T3l 

M, < 700 Mq is allowed at the 1-cr level, and M, < 1500 Mq 
at the 95% (< 2-a) confidence level. This conclusion is es- 
sentially independent of the assumed value of Wo, as Table 
[nialso shows: for Wo = 8.1 and Wq = 9.1, the best-fit black- 
hole mass is also zero, and the 68%, 95%, and 99% upper 
limits are very similar to those found with Wo = 8.6. All in 
all, the proper-motion dispersions cannot be used to argue for 
the presence of an intermediate-mass black hole at the center 
of 47 Tuc — although neither do they strongly disallow such a 
possibility. 

What is perhaps more convincingly seen in Tablell3lis the 
fact that the velocity-dispersion profile of these stars (with 
(m») ~ 0.84Mq) always prefers a "core" radius, ro, which 
is rather larger than the w 2 1 " implied by the density profile 
of the same populatio n. Th is is yet another demonstration that 
(as was discussed in ^5. 21 and ^6.3> the bright stars in 47 Tuc 
appear to be moving in response to a more extended mass dis- 
tribution that is reminiscent of what we naively expect for the 
dominant, but unseen (in this cluster), population of stars with 
~0.5Mq. 

Figure l25l illustrates these results. In the upper panel, we 
show the number density of turn-off mass stars (F475 W mag- 
nitudes < 17.8), constructed from our master star list (Table0} 
as described in ^4.11 but now binned in annuli with AR =1" 
for R < 20" and AR = 3" for R > 20". The bold, dash-dot 
curve is the Wo = 8.6 Kine-model fit from Figure^] while the 
dotted curve has the larger ro = 27 "5 from Figure l24l above. 
The bottom panel of Figure |25] then shows ci^ vs. R for our 



TABLE 14 

Central Black-Hole Model Fits to V < 20 Kinematics " 



GM./alro 




O"0 


'■o 










[mas yr~' ] 


[arcsec] 


[Mq] 










Wo = i 


i.6 






0.00 


97.74 


0.613 


23.06 







0.01 


97.24 


0.612 


23.03 


140 


<68% C.I. 


0.02 


96.89 


0.611 


23.21 


282 




0.03 


96.65 


0.610 


23.53 


427 




0.04 


96.49 


0.609 


23.92 


576 




0.05 


96.42 


0.607 


24.36 


729 


best fit 


0.06 


96.44 


0.605 


24.83 


887 




0.07 


96.56 


0.603 


25.32 


1050 




0.08 


96.78 


0.602 


25.85 


1220 




0.09 


97.12 


0.600 


26.39 


1390 




0.10 


97.55 


0.598 


26.95 


1570 


> 68% C.I. 


0.11 


98.11 


0.596 


27.52 


1750 




12 




594 


28 08 






0.13 


99.63 


0.592 


28.61 


2120 




0.14 


100.61 


0.590 


29.12 


2310 


>95% C.I. 


0.15 


101.73 


0.588 


29.57 


2490 




0.16 


103.00 


0.586 


30.01 


2680 


>99% C.I. 


0.17 


104.41 


0.584 


30.37 


2870 




0.18 


106.00 


0.583 


30.68 


3050 








Wo = i 


i.l 






0.00 


97.79 


0.613 


21.61 







0.01 


97.32 


0.612 


21.59 


131 


68% C.I. 


0.05 


96.33 


0.609 


28.45 


856 


best fit 


0.09 


97.44 


0.600 


31.06 


1640 


> 68% C.I. 


0.13 


100.80 


0.591 


33.95 


2510 


> 95% C.I. 


0.15 


103.49 


0.587 


35.16 


2950 


> 99% C.I. 






Wo=' 


).l 






0.00 


97.87 


0.612 


19.60 







0.01 


97.44 


0.611 


19.63 


119 


< 68% C.I. 


0.06 


96.53 


0.604 


21.09 


752 


best fit 


0.11 


97.60 


0.596 


23.20 


1470 


68% C.I. 


0.15 


100.29 


0.590 


24.75 


2100 


95% C.I. 


0.18 


103.58 


0.585 


25.61 


2560 


> 99% C.I. 



^Fits are to the proper-motion velocity-dispersion profile of stars with R < 90", 
V < 20, P(x^) ^ 0.001 in both RA and Dec components of proper motion, and 
uncertainties of < 0.25 mas yr"' in both components. Stars are binned in l"-wide 
annuli, resulting in 87 points fitted. 

^Dimensional black-hole masses assume a distance of D = 4.0 kpc to 47 Tuc. For 
any other distance, multiply by (D/4.0kpc)"^. 



proper-motion stars with V < 18.5, calculated in annuli with 
the same AR inside and outside R = 20" as in the upper panel. 
The model curves drawn in the bottom panel of Figurel25lcor- 
respond to the best-fit M, = and the rough 68%, 95%, and 
99% upper limits on M, for the Wo = 8.6 models detailed in 
Table [T^ As was mentioned above, these curves were ob- 
tained by fitting the cr^ profile calculated in l"-wide annuli at 
all radii R < 90"; the coarser binning at large R in this figure 
is applied simply to make the plot clearer. 

We have also defined an alternative velocity sample by se- 
lecting all stars from Table |5l which have V < 20, P(x^) > 
0.001 in both components of proper motion, colors obeying 
equation ( I18> . and velocity errorbars and As both less 
than 0.25 mas yr"'. The last of these selection criteria allows 
for more robust determination of the error-corrected veloc- 
ity dispersions in this sample. We then fit our modified Kin^ 
models with central black holes to the true (Tf_i{R) in l"-wide 
annuli, exactly as we did before. 

Table [T4l which has the same format as Table [T3l shows 
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Fig. 26. — As for Figure l25l but now for stars with V < 20. Note the better 
description provided for the density profile of this sample (which includes 
stellar masses down to m* ^ 0.65 ; Table[7) by the ro = 27;'5|King-model 
curve in the upper panel. The kinematics sample in the lower panel is further 
restricted by placing upper limits of 0.25 mas yr"' on the uncertainty in both 
velocity components [as well as the usual limits on P(x^) and the ( U — V) 
colors]. Curves in the bottom panel are again fits of Wq = 8.6 black-hole lKind 
models to the observed (Tp(i?). The bold, solid curve is the formal best fit, the 
lowermost curve has M, = 0, and the other dashed curves represent the 68%, 
95% , and 99% upper limits on M, . Details of all these fits are given in Table 

El 

that the best-fit black-hole mass inferred from fitting this 
different sample of proper-motion stars is now M, ^ 700- 
800Mq (depending slightly on the assumed value of Wq). 
The 68% confidence interval is about (IOOMq, 1500 Mq) and 
the 95% confidence interval is roughly (O,23OOM0). At best 
this might suggest a ~l-cr consistency with the hypothesis 
of an intermediate-mass black hole, but the fact that different 
velocity-sample definitions yield formal best-fit masses that 
differ by nearly lOOOM© clearly emphasizes the weak nature 
of any black-hole signal in our data. It also has to be stressed 
that the sample of stars with V < 20 covers a non-negligible 
range of masses (m* ~ O.65-O.9M0; Table0, and moreover 
that the mixture of mass classes varies significantly with ra- 
dius as a result of selection on the basis of errors (see Fig.|9j — 
effects that are ignored in our simple models. 

Figure |26l shows, in its upper panel, the completeness- 
corrected density profile of stars with estimated V < 20 in our 
master star list, Table0] As in Figure|25] King ( 1966) models 
with Wo = 8 .6 and ro = 20"84 or ro = 27 "5 are drawn here. The 
latter now provides a better description of the spatial distribu- 
tion in this closer-to-average mass range. The lower panel 
of Figure |26l displays the true crp(/?) profile and the best-fit 
black-hole King model with Wo = 8.6 from Table [T4l along 
with the 68%, 95%, and 99% limits on M.. 

Whatever small signature of a central mass concentration 
might appear in Figures |25l and |26l obviously comes from 
the innermost R < 3" — or 0.06 pc, which would contain the 
sphere of influence of any < 2000-Mq black hole — and in 
particular from the high apparent cr^ in the central 1". How- 
ever, there are only 8 stars at R < 1" in the V < 18.5 veloc- 



ity sample and 1 1 stars in the error-selected V < 20 sample, 
which is why the significance of the rise is unavoidably low. 

Sample size aside, it is also clear that much more sophis- 
ticated modeling is required if these proper-motion data are 
to be used either to prove or to disprove the presence of an 
intermediate-mass black hole in 47 Tuc. Allowances must be 
made for a nontrivial stellar mass spectrum; careful consid- 
eration must be given to sample selection effects; and, most 
likely, higher-order moments of the velocity distribution will 
need to be examined in addition to the dispersion. On this last 
note, we now construct the full proper-motion distributions in 
the central 3" of the cluster. 

6.5.2. Velocity Distributions in the Innermost 3" 

Figure]^ shows the two-dimensional velocity distribution 
N(fitot) for stars with R < 3" in the two subsamples. These 
are o btained from the data in the same way as described in 
^5.21 The solid curves in both panels are the error-convolved 
distributions for regular Wq = 8.6 Ki ng m odels with no cen- 
tral black hole, also computed as in ^5. 21 but with ro and ctq 
now set by our fits to the velocity-dispersion profile alone (i.e., 
taken from the first line of Table II 31 for the upper panel, and 
the first line of Tablell4lfor the lower panel). The dotted lines 
are the one-sigma uncertainties on the predicted A^(/itot)- By 
and large, the data and the model distributions with no black 
hole agree reasonably well. The absence of a handful of ex- 
pected stars with /itot ^0.2 mas yr"' in the upper panel may be 
of some interest, but it necessarily has low significance given 
the small sample size. In any event, note that none of the 
high-velocity stars discussed in ^5.3l fall in the zone R < 3". 

Finally, Figure|28]shows the one-dimensional distributions 
of the separate radial and tangential components of proper 
motion for stars with R < 3" in the same two samples as be- 
fore. The observed distributions are obtained in the same way 
as the Niiitoi) distributions we have focused on to this point; 
in particular, the errorbars in Fig. l28lare obtained by the same 
bootstrap procedure described in ^5.21 For a given veloc- 
ity sample, the error-convolved model curves for N(ij,r) (left- 
hand panels of Fig.|28} are identical to those for A^(/ie) (right- 
hand panel s), due to the assumed isotropy. They are calcu- 
lated as in ^4.2.11 only using the appropriate one-dimensional 
Gaussian distribution of velocity errors rather than eq uation 
ilO\ and with the function N](ij,_^\R / ro) from equation ( IB14> 
taking the place of A^2(AiA , Aiy|^/'"o) in equations il 1> and J12l i. 
As in Fig.|^ we use do and ro from the first lines of Tables 
ElandOto compare the models and the data. 

The excess of th e observed velocity dispersion over that 
expected in a iKind ri966.) model with no black hole is per- 
haps more in evidence in these plots than in Fig. |^ The 
broader wings of the tangential N(^iq) distributions for both 
samples are particularly apparent, although it is also clear that 
the discrepancies are at the ^l-cr level. The distributions of 
radial proper motion in the left-hand panels may show more 
interesting features, but again these have low significance. 
Specifically, there appears to be an asymmetry in N(ij,r) in 
both velocity samples. However, in the brighter sample with 
V < 18.5 this is due to the presence of just two stars moving 
radially inwards with jiR < -1.3 mas yr"', vs. none moving 
outwards with the corresponding /i^ > +1.3 mas yr"'. In the 
fainter sample, the asymmetry reflects the presence of four 
stars with /ij? < -1 .3 mas yr"' vs. two with /ir > 1 .3 mas yr"' . 
Similarly, the very high observed dispersion ctr in the second 
sample (which drives the non-zero best fit ofM. in Tabled 
and Fig.l26> is due in large part to a single star with /ir ~ -1 .9 
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Fig. 27. — Two-dimensional velocity distributions within 3" of tlie clus- 
ter center for t he prop er-motion samples to which black-hole models are fit 
in Figures l25landl2^ Solid curves are en'or-convolved Kins 1 1966) model 
N{fiiot) distributions for the Wo = 8.6 fits with no black hole in the two pre- 
vious figures (top lines of Tables ll3l and [T4l . The dash-dot curves are the 
±1-(T uncertainties on the model. The observed (error-corrected) ct^ indi- 
cated overlap with the model expectations for no black hole at the one-sigma 
level, and more generally the observed full distributions appear consistent 
with the model curves at the same level (except for a paucity of slow-moving 
stars in the upper panel). 
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Fig. 28. — One-dimensional velocity distributions within 3" of the cluster 
center for t he prop er-motion samples to which black-hole models are fit in 
Figures|25land l26l Left-hand panels show ^(/ir) for stars in the two samples 
of Fig. 1271 while right-hand panels display N(p0) for the same samples. 



mas yr ' and one with mas yr ' — neither of which 

is expected in a regular lKind (or Gaussian) model for a sam- 



ple of this size, but either of which could simply be a sta- 
tistical oddity. It will be challenging indeed to build more 
comprehensive velocity data sets, in this or any other globular 
cluster, which might be able to overcome inherent statistical 
limitations such as these and discern unambiguously the sub- 
tle kinematical effects of any central mass concentration with 
M, ^ 10^ Mq or less. 

7. SUMMARY 

We have used multiple HST/WFPC2 and ACS images of 
the inner regions of 47 Tucanae (NGC 104), obtained be- 
tween 1995 and 2002, to derive proper motions for 14,366 
stars within R < 100" of the cluster center — the largest veloc- 
ity data set yet collected for any globular cluster At the same 
time, we have used a single set of dithered ACS images to 
construct an unbiased list of photometry and astrometrically 
calibrated positions for nearly 130,000 stars in a ~ 3' x 3' 
central field which contains the proper-motion sample. This 
"master" star list is presented in Table 0] above ( ^2.2> . while 
the data going into our proper-motion determinations are in 
Tablel5l( ^3.ll i. Both of these catalogues can be downloaded 
in their entirety through the electronic edition of the Astro- 
physical Journal. Supplementary online material is also avail- 
able: an image (4 7TucMaster . fits) of a circular region 
R < 150" around the cluster center, which contains the area 
covered by the master star list of Table |3 and an SM macro 
(pmdat .mon) for extracting and plotting position vs. time 
data from Table|5]for any star in the proper-motion catalogue. 

Section|2]described in detail the basic procedures and qual- 
ity checks which we used to obtain calibrated astrometry 
and relative position shifts (from local coordinate transfor- 
mations) for many thousands of stars with multi-epoch imag- 
ing from different HST instruments with different pointings, 
fields of view, and exposure times. Along the way, we used 
the unbiased master star list to estimate the coordinates of the 
cluster center, finding our result to be in good agreement with 
previous values from the literature (see Table[0. 

In |3lwe defined the subset of stars for which we derive rel- 
ative proper motions, in which the mean velocity of the entire 
cluster is zero by definition. We applied standard weighted, 
linear least-squares fitting to find the plane-of-sky velocities 
and uncertainties, and to estimate the quality of a straight- 
line fit to each component of motion. We found that some 
~10% of the 14,366 proper-motion vectors we obtained were 
either too uncertain, or not linear to high enough confidence, 
to include them in subsequent kinematical analyses. How- 
ever, these stars were retained in the velocity catalogue of Ta- 
ble |5] Although we have not examined the issue in this pa- 
per, in some cases the poor quality of linear proper-motion 
fits may reflect truly nonlinear motion (rather than simple 
measurement errors or noise) and thus be of physical inter- 
est. Section |3] also gave a brief description of a set of radial- 
velocity data — some already in the literature, and some as yet 
unpublished — which we eventually used only in order to esti- 
mate the distance to 47 Tuc. 

Section 13 used stars in our unbiased master list again (not 
the velocity sample, which has a highly irregular and biased 
spatial sampling) to define the surface-density profile of stars 
with the main-sequence turn-off mass in 47 Tuc (those with 
magnitudes V < 17.65, or F475W < 17.8, corresponding to 
TO, ~ O.85-O.9M0) from R ~ 70" to R = 0. After joining 
this profile on to ground-based V-band surface photometry 
from the literature we fit the cluster structure to /? w 1000" 
with a single-mass, isotropic iKind Jl966h model for a low- 
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ered isothermal sphere. Our best-fit core radius, central sur- 
face brightness, and central concentration are listed in Table 
[l] and are in good agreement with other estimates in the lit- 
erature. This model does not fit the density profile of 47 Tuc 
from R > 1000" to the tidal radius R ^ 3000"; a realistic 
stellar mass function and/or velocity anisotropy must be taken 
into account to explain the outer halo structure. However, it 
does accurately describe the R < 100" field covered by our 
proper-motion data, and thus we have used it as a convenient 
point of reference for interpretation of some of the observed 
kinematics. In particular, we have derived the projected one- 
and two-dimensional velocity distribut ions as a function of 
clustercentric position in any isotropic King! model ( ^4.2l and 
AppendixO, which could be profitably applied to other large 
velocity samples in star clusters. 

In ^we presented the one- and two-dimensional velocity 
distributions for our proper-motion sample and found them 
to be very regular overall. For speeds less than the nominal 
central escape velocity in 47 Tuc (/^max — 2.6 mas yr"', cor- 
responding to ~ 50 km s"' for a distance of 4 kpc), all distri- 
butions are well described by the single-mass, isotropic Kins 
model functions that we have derived, after convolution with 
the velocity measurement errors and the stars' spatial distri- 
bution. In turn, these model velocity distributions are nearly 
Gaussian in the core regions of 47 Tuc, since this massive 
cluster defines a rather deep potential well. We compared the 
observed and model distributions as a function of radius for 
the brightest stars with the smallest proper-motion uncertain- 
ties, showing that the velocity dispersion decreases slightly 
more slowly with clustercentric radius than expected for a 
single-mass King model ( ^5.21 Figure l20l. We argued that 
this implies that the mass distribution in the inner few core 
radii of 47 Tuc is dominated by stars less massive than the 
main-sequence turn-off mass, rather than by heavy remnants 
such as neutron stars or massive white dwarfs. This is consis- 
tent with other suggestions, from independent modeling and 
observations, that the total mass of such remnants is quite low 
in this cluster. 

In ^5.3l we took especial care to examine the properties of 
stars with apparent proper-motion speeds above the nominal 
central escape velocity of 2.6 mas yr~' . Only a few dozen such 
stars are found, amounting to < 0.3% of our total sample. Of 
these, fewer than ten, or < 0.1%, can be said unequivocally 
not to have been scattered to high velocities by measurement 
errors, and it is still possible that some of these may be in- 
terloping field stars. Thus, we have not found any evidence 
for a significant population of very fast-moving, "nonthermal" 
stars in the core of 47 Tuc, such as had been suggested on the 
basis of some earlier radial-velocity surveys involving much 
smaller samples. We emphasize, however, that we have not 
examined in any detail the ^ 5% of stars in our catalogue 
which are not fit well by straight-line proper motions. It is 
conceivable that some legitimate nonlinear motions among 
these stars could be of physical interest in connection with 
the close encounters and scattering processes expected in this 
dense stellar environment. 

Section|6llooked at a few aspects of the proper-motion ve- 
locity dispersion only. In ^6.11 we showed that the velocity 
dispersion of 18 blue stragglers within R < 20" (one core ra- 
dius) is smaller than the dispersion of comparably bright stars 
on the cluster red-giant branch, by a factor of about \/2. To- 
gether with other studies in the literature showing that blue 
stragglers are more centrally concentrated than the giants in 



47 Tuc, this is consistent with the blue stragglers being on av- 
erage about twice as massive as normal main-sequence turn- 
off sta rs. 

In ^6.21 we presented the run of proper-motion disper- 
sion with clustercentric radius for stars in various magnitude 
ranges, compared it with the expected behavior in a single- 
mass King ( 1966) model, and showed quantitatively that the 
average velocity anisotropy in the inner R < 100" ~ 5 core 
radii of 4 7 Tuc is — as expected — negligible. 

Section l63l compared the line-of-sight velocity-dispersion 
profile defined by 419 bright giants (V < 14) within R<105" 
with the proper-motion dispersion profile of 6,336 of our stars 
with V < 18.5 (and, thus, masses comparable to the giants). 
This led to a kinematic estimate of the distance to the cluster: 
D = 4.02 ± 0.35 kpc. While this result was obtained by ignor- 
ing the rotation of 47 Tuc, on the small spatial scales probed 
by our kinematics the rotation is quite weak and should not 
affect our estimate of D by more than 2-3 percent at most 
( ^6.4> . Our distance is consistent with results in the literature 
from fitting to the white-dwarf cooling sequence. We also 
used a King model to extrapolate the velocity dispersions to 
R = 0, giving improved estimates of the central and cr- (see 
TablelD. 

Finally, ^6.5l focused on the kinematics and density struc- 
ture in the innermost few arcseconds of 47 Tuc. We fit 
the proper-motion velocity-dispersion profile with simplistic 
models based on those of King (1966) but allowing for the 
presence of a central point mass. We found no clear need here 
for an "intermediate-mass" black hole of the type recently 
claimed to have been detected in other globular clusters. Nev- 
ertheless, one with a mass M, < 1000-1500 Mq could be ac- 
commodated by the data at the l-a level. It would thus seem 
difficult to use these data alone either to prove or to disprove 
whether the well known correlation, between supermassive 
black-hole mass and nuclear velocity dispersion in galaxies, 
extends to globular clusters. 

The overarching conclusion from this work is that the core 
of 47 Tucanae is remarkably unsurprising in its dynamical 
structure. The velocity distribution is essentially isotropic 
and close to Gaussian, with radial and (secondarily) mass 
dependences which — aside from a few relatively small de- 
tails which can be at least qualitatively understood — ^roughly 
match the expectations for standard King-type models of dy- 
namically relaxed, lowered-isothermal spheres. Our work was 
aimed expressly at an investigation of the core of this cluster, 
and our results have emphasized that any truly new insight 
into star-cluster dynamics from large proper-motion samples 
may be more likely to come from surveys of more off-center 
fields, where the potential is weaker, relaxation may not be so 
advanced, and velocity anisotropy in particular could be much 
more pronounced. 
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APPENDIX 

A. CORRECTING VELOCITY DISPERSION FOR MEASUREMENT AND TRANSFORMATION ERRORS 

This Appendix addresses the need to interpret our observed velocity dispersions in terms of the true dispersion (of the intrinsic 
stellar velocity distribution) and the additional contribution made by our measurement errors. There are two sources of error. 
The first, which is random, relates to the effect that simple errors in position measurements have on the observed proper motions. 
The second, which is systematic, comes from our use of local transformations to bring a star's coordinates at any epoch into our 
master-frame system (see ^2.3.2t . 

Recall that by applying standard linear regression to infer proper motions from our measurements of displacement vs. time, we 
have assumed that the position measurement error for each star at each epoch is Gaussian distributed with a mean of 0, in which 
case the errorbars on the fitted /i^ and jig represent the dispersions of Gaussian distributions of RA- and Dec-velocity measure- 
ment error The errorbars on the two velocity components need not be the same for any one star and are in general different from 
one star to the next, but the mean measurement error is always assumed to be 0. Thus, the distribution of measurement error in 
either component of velocity for a set of any M stars from our sample is given by the normalized probability density 



2Ai 



(Al) 



where S/j, represents the error in either the RA or the Dec component of the observed proper motion and A; is the uncertainty in 
that component for the /* star in the sample, as given by our least-squares fitting. 
Any observed velocity /iobs can then be viewed as the sum 

Mobs = /^coi- + <5Ai , (A2) 

where ^coi is a random variable drawn from the (a priori unknown) error-free velocity distribution as a function of stellar position 
and magnitude in 47 Tuc, and the measurement error (5/i is a random variable drawn from the distribution (IA1> . If the measurement 
errors are uncorrelated with the velocities themselves (as the lack of correlation between the observed fia and the least-squares 
uncertainties Aq, in Figure|5]of ^3.1.2l implies is the case), then the expectation value of the observed RA or Dec velocity for a 
set of J\f stars is simply 

E(pobs) = E(p^or)+EiSfJ,) = Eipcor) , (A3) 

since E{5ij) vanishes by virtue of equation (lAU . Thus, an unbiased estimate of the mean velocity with no measurement error is 
provided by the usual 

1 ^ 

(/^)coi- = (M)obs = ^ X! ^'>''^ ■ ^^^^ 

For uncorrelated random variables, the variances are similarly additive, regardless of the form of the underlying distribution(s) 
of the variables. Thus, 

Vai-(Aiobs) = Var(/icor) + Var((5Ai) , ( A5) 

or 

£[(Mobs-(M>obs)']=£[(Mcor-(A'>cor)']+£[('5Ai)'] • (A6) 

Switching now to more standard notation, we write (j^^^ = Var(/^obs) and al^^ = Var(/i(;oi-); use equation (lAU to evaluate Var((5/x); 
and finally use the sample standard deviation as an unbiased estimate of ciobs to obtain 



2 _ 1 

N- 1 



■ j\r 

y~^(M',obs~ (A^)obs)^ 
. 1=1 



JzE^r (A7) 



1=1 



That is, the proper estimator for the error-corrected velocity dispersion (Tcor of any M stars with unequal velocity uncertainties 
is obtained by subtracting the rms errorbar in quadrature from the usual standard deviation of the observed velocities. The only 
assumption behind equation (IA7> is that the velocity errors are normally distributed, and thus it holds true for any form of the 
actual velocity distribution. In particular, the latter need not be Gaussian — indeed, it need not be specified at all. 

There remai ns one further correction to be applied to (Tcor in order to reach an estimate of the true velocity dispersion: As was 
mentioned in ^2.3.21 the use of local transformations to map stellar positions at our ten epochs into the master-frame coordinate 
system introduces an unavoidable error in the observed velocity dispersion, which is completely separate from the inflation of 
dispersion due to measurement error. 

Suppose that, as in ^2.3.21 we use the positions of the A/Jrans nearest neighbors of a star at some epoch to define the local 
transformation of that star into the reference-frame coordinate system. Suppose further that there are no position-measurement 
errors in any frame at any epoch. Then the formal error in all local transformations is zero, and the mean and dispersion of the 
observed velocities are automatically (/i)cor and dcor in the notation of the preceding discussion. How do these relate to the true 
moments of the intrinsic velocity distribution? 
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To address this question we look, without loss of generality, at a single component (any component) of position in the master- 
frame system, measured for a given target star and for its nearest A/trans neighbors at any two epochs separated by At in time. Call 
these positions xq for the target star and {xj; j =1,... , Attans} for the neighbor stars. By applying local coordinate transformations, 
the proper motion /icor we derive for a star is really measuring the change in its position relative to the centroid of its nearest 
neighbors at the two epochs. That is, if 

1 

A/trans ^ 

denotes the centroid of the neighbor stars, then 

/icor At = {x'Q- x'^^J - (Xo - .Xcen) = (x'q-Xq)- {x[^^ - .Xcen) , ( A9) 

where the primes identify the second-epoch positions. Then, noting that the expectation value of the posit ion c hang e for any star 
is E(x' -x) = AtE(^irue) (where /itnje refers to the true distribution of absolute proper motions), equations ( IA9t and ( IA8> together 
imply that 

£(Mcor) = . (AlO) 

This only says, quite sensibly, that on average the velocity of a star relative to its nearest neighbors will be 0. This is, of course, 
exactly what the true situation is, so that even though the local transformation scheme leads to an error in every observed stellar 
velocity, the average error is 0. 

By contrast, the relationshi p be tween the observed dispersion dcor and the second moment crtrue of the true velocity distribution 
is nontrivial. From equation (IA9> we have 

(Af)^Var(//cor) = Var(x(,-xo) + Var(x^g„-Xcen) 

= (Af)2Vai-(/i,„e) + (Af)'Var(/i,e„), ^ ' 

where /icen = (1 /^timis)jyj^'i"ix'j—Xj)/ At from equation ( IA8> . Clearly, Var(/itrue) = ftmc while the variance of /i^en follows from 

its definition as the sample mean for the A/Jrans neighbors of our target star: Var(/i(;en) = ^t^rueZ-'^trans- Finally, then, equation 
yields 

<^L=(^ + J—)<^L- (A12) 



Aft. 

Thus, from the combination of equations (IA4> and jA10> we expect to find 

1 ^ 

(M)obs = Jf'^^^Uohs = (A 13) 

for the mean relative proper motion (in any direction) of any set of stars; and by combining equations ( IA7> and (IA12> we calculate 

AT 



crL. = ( 1 + ■ 



A/irans / A/"" 1 



L i=l 



^(A*i>bs- (A*)obs) 



AA^ ' 

/=1 



(A14) 



to estimate t heir t rue one-dimensional velocity dispersion. 

Equation (IA14t . with A/Imns = 45 following the discussion in ^2.3.21 is used repeatedly in our analyses of ^and ^ to find the 
intrinsic RA and Dec or radial and azimuthal velocity dispersions for various subsets of our proper-motion sample. 

B. PROJECTED VELOCITY DISTRIBUTIONS IN KING (1966) MODELS 

It is well known that the velocity distribution at any radius in a lKind (119661) model cluster takes the form of a "lowered 
Maxwellian," i.e., a Gaussian distribution in v = (vj + v^ + v^)'/^ minus a position-dependent constant related to the local escape 
speed. When comparing such models to data, however, it is necessary to integrate this simple distribution over all phase-space 
coordinates for which we have no direct information. This will always include at least a projection along the line of sight, which 
is already nontrivial. In a case where only proper-motion data are being considered, an additional integration over one component 
of velocity is also required; while if only radial-velocity data are being modeled, yet another integration over a second velocity 
component is necessary. In this Appendix, we develop general expressions and outline a computational scheme for these various 
marginalized velocity distributions. 

If £■ = v^/2 + (j){r) is the specific energy in a self-gravitating cluster of equal-mass stars, with (l)( r) a relative potential defined 
so that = at some tidal radius to be determined, then the phase-space distribution function of a lKind Jl966t model (see also 
iBinnev & Tremainell987l Section 4.4) is defined as 

U(2woV/'(e-^/-o-l) , E<0 ^gj^ 
\ , £>0, 

where ao is a velocity scale and g is a normalization factor related to the central space density. The "lowering" term (-1) in the 
top line allows for the possibility of a finite escape velocity: this distribution function only allows for stars with E <0, and thus 



V < y™ax(r) = V^2^ . (B2) 
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The usual space density p(r) then follows from integrating f(E) over all allowed velocities: 

P(r) = /(£)47rv2t/y=gi|e'*'Werf[VWW]-\/4W(r)/7r [l+2W(r)/3]| , (B3) 

where W{r) = -4>{ r) /ctq is the dimensionless relative potential. If Wo is the (positive, but otherwise arbitrary) value of W at r = 0, 
then equation IB 3 1 defines g in terms of Wo and the central density po- The full density and potential profiles are obtained from 
Poisson's equation in dimensionless form: 

— {r^—)=-9^r\ (B4) 
dr \ dr J po 

in which 

7= r/ro , ri = 9al/(4TrGpo) . (B5) 

Equation ( IB3> is used to write p(F)/po in terms of W(T) and Wo, after which equation ( IB4> is integrated numerically [subject to the 
specification of Wo and the additional boundary condition (dW /dTyf^o = 0] to find W and p / po as functions of r/ro . The radius at 
which p{r) vanishes defines the tidal radius of the clu ster, and the concentration parameter c = log(r, / ro) is a one-to-one function 
of Wo (the iKing model fit to 47 Tuc in Figure^lof ^4.11 with c = 2.01, has Wo = 8.6). Subsequent projection of p{r) along the 
line of sight z is straightforward, yielding the surface density as a function of projected radius R = (r^ — z^'fl'^: 

m = 2r"'^ W)dz. (B6) 
Pora Jo 

where R = R/r^, z = z/ro, p = p/po, and 7= r/ro is evaluated as (R^ + z^y^^ in the integrand. While p(r) and S(/?) are actually 
mass densities, for a single-mass cluster they are directl y pr oportional to the stellar number-density profiles, and to the luminosity 
density or surface-brightness profiles. It is equation iB6\ that was fit to the surface-brightness profile constructed in ^4. II to 
constrain Wo, ro, and po (and thus to derive c and pv.o) for 47 Tuc. 

All of the above is completely standard and has been repeated here simply to give some basic definitions and to establish 
our notation. The usual next step in kinematical applications of these models is to obtain the dimensionless velocity-dispersion 
profile, a,-/(Jo vs. r/ro, from the definition 

3p(r) Jo 

which impUes (using eqs. IB II and IB2I . and eq. IB3I evaluated at r = 0) 

a2(?)_ 1 1 fV^y exp[W(?)-y72a2]-l 



^'('•)=TT\ / v^fiE)4nv^dv , (B7) 
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crl ^p(r) Jo exp(Wo)erf (VW)) - ^4Wo/7r [1 +2Wo/3] cTq cto ' 

After finding the profile cr^ = /ctq' it is again straightforward to integrate it along the line of sight to obtain the projected velocity 
dispersion 

p(7) aj(7) dl (B9) 

for direct comparison with radial-velocity observations. More generally, because of the assumption of velocity isotropy in these 
models, is als o the dispersion o f any component of proper motion as a function of projected clustercentric radius. Note that 
neither equation ( IBSt nor equation ( IB9> requires that the parameter cto correspond exactly to the actual velocity dispersion at the 
center of the cluster; and indeed, cr, (r = 0) ^ (To and <7~(R = 0) ^ tio in general for these mode ls. ( It is only in the limit of infinite 
W) that the central dispersions tend to cto). Section |6l compares the predictions of equation ( IB 91 . for models appropriate to the 
spatial structure of 47 Tuc, to the observed dispersion profiles of our proper-motion and line-of-sight velocity samples. 

We are additionally interested here in the full distribution of projected velocities, not just the second moment; and we have 
data for the joint distribution of velocities in two orthogonal directions on the pla ne of the sky, not just a single component. Thus, 
consider first the distribution of th ree-d imensional velocity implied by equation (IBU . At any fixed radius r, this can be written 
(using the definition of Vmax in eq. IB2h as 



ge 



<J2< r 



n3(v.,v,,v,|r)= ;; ,-(.'H-^/2-5_,-L/2.„^ vi+vi+vi<vL. (bio) 

with v^2x/2co = W(r) already known as a function of radius from the procedure outlined just above. Recall that 
/// nT,(v;t,Vy,v~)dvjcdvydv. = p(r) by d efinition of the phase-space density /. 

Marginalizing the distribution JB IOI over the two velocity components Vx and Vy (for example) gives the velocity distribution 
at r in the v- direction alone: 

ni{vz\r)= / / n3(vx,vy,v,\r)dvxdvy 

J J-VmiJf) 
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whic h is the n umber of stars per unit z-velocity and per unit spatial volume, at radius r in the cluster. Alternatively, equations 
(IB2> and (IB3> can be used to write v^^^ and g in terms of W(r) and p{r), so that 

, I, P(r) exp[W(r) - v^/2ag] - 1 - [W(r) - v?/2ag] 

n\{y.\r) = — = , , . (B12) 

%/2^CTo exp[VK(r)]erf[VWX^]- V4W(^[l+2VK(r)/3] 



By virtue of the assumed isotropy, the distributions n\{vx\r) and ni(vy\r) are formally identical to equation JB12t . Clearly, 
ni(vmaxk) = always. 

For comparison with observations, equation (IB12> must be averaged along the line of sight. Denoting this coordinate by 
z = (r^-R^yl^ again, and normalizing to a density of one star per unit area at projected radius R on the plane of the sky, we define 

2 Jo' p{r)dz 

or, in terms of the model's dimensionless radius, potential and surface density, 

^^^^^ ^^ = _^ r--'^'" exp[H^(^-v^/2ag]-l-[H^(7)-v,?/2ag] ^ ^^^^^ 

\/2^(7o ^0 exp(Wo) erf ( ^1%) - VWo/tt [ 1 + 2W)/3] 

where E is the dimensionless 'S{R)/poro given by equation (IB6> and, as also in that definition, R = R/ro, z = z/ro, and 7= r/ro = 
(R^ + z^y^^. Note that this distribution has dimensions of inverse velocity and is normalized such that jNi(v^\R) dv- = 1 at any R. 

Thus, A^i (v,\R) is the probability that a single star at projected clustercentric radius R/ro will have a line-of-sight velocity between 
and (vj + c/vj). As with the volume probability distribution ni{v-\r ), the corresponding distribution for any one component of 

velocity on the plane of the sky has the same form as equation JB14> . 
The two-dimensional distribution of velocities projected onto any spatial plane can similarly be derived from equation JBim . 

Obviousl y prop er motions are the motivation for looking at this, so let us now write /i^ = Vj and /i, = Vy. Then integrating 

equation (IB10> over all velocities (with |v,| < Vmax as usual) and defining the auxiliary variable 

U(r) = = W(r)-^\^ (B15) 

2CTg 2CTq 2CTg 

gives 

,^ p(r) exp[t/(r)]erf[Vr7(7)]-v^t7WA 

«2(Mx,Mv k") = T ^= , : (BlO) 

2nal exp[W(r)]evnVWi^]-V^W0^[l+2W(r)/3] 

for the number of stars per unit velocity element dp^dpy, per unit spatial volume, at radius r in the cluster. Through U(r), this 
distribution is a function only of the total proper-motion speed, /itot = ip^ + fi^,y^^ (i.e., it is circularly symmetric in the Px-fJ-y 
plane of phase space), again because of the assumed velocity isotropy. It vanishes at any radius r if p,^^j2<jl > W{r). Averaging 
along the line of sight as in equation ( IB13I I — including the normalization there to one star per unit area on the sky — then gives 
the observable 

2t:<jI^(R)Jq exp(iyo)erf(VW^)-V^HV7^[l+2Wo/3] 

for the dimensionless model quantities E, R, z, and Fall defined as above. Note that Nxip-x, Py\R) has units of inverse velocity 
squared (it is analogous to a spatial surface density), and it is normalized so that its integral over all allowed /ii and py (at any 
fixed R/ro) is always unity. Finally, note that p^ and /i,. could, in fact, be any two orthogonal components of velocity without 

changing any of the results here. In particula r, the joint distribution of radial and tangential proper motions, N2{pR, pelR), has 
exactly the same functional form as equation ( IB17> . 

Unfortunately, equations ( IB14> and (IB17> can only be evaluated by numerical integration. To do this, for any given Wo we 
first compute look-up tables for the dimensionless potential W{7^ and surface-density profile E(/?). Then we define a series of 
dimensionless projected radii ranging over < R / ro < r, / ro, where the upper limit r, / ro is determined by the value of Wo. Next, at 
each point in this radius sequence we define a series of dimensionless velocity values covering the interval < u^/2ao < W{R/ro), 
where can be associated either with or with (/i^ + /i^). Finally, we perform the integrals JB14> and ( IB17> at each pair of fixed 

(R/ro,u-/2<Tl), ultimately obtaining two-dimensional look-up tables for both A^i(Vj;|/?) and A^2(Mj:iMvI^)- 

Application to a real velocity sample further requires convolution both with the observed spatial distribution of the stars in 
the sample and with the distribution of their velocity-measurement errors. Section l4!2l describes how we have done this for 
N2{fJ-x,P'y\R) in particular, and § jSl and lOl compare the results to our observed proper-motion distributions. 
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TABLE 4 

129,733 Stars in the Master Frame 



ARA 


ADec 


F475W 


RA 


Dec 


RA 


Dec 




OK 


c 


ID 


.^meta 


ymeta 




>'ref 


[arcsec] 




[hours] 


[deg] 


[hh:mm:ss] 


[dd:mm:ss] 






[%] 




[pixe 


s] 


[pixels] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


88.1030 


-100.8748 


21.51 


0.4068860 


-72.109304 


00:24:24.789 


-72:06:33.49 


25.39 




97.51 


MOOOOOl 


1238.94 


983.50 


4078.27 


4167.41 


85.7765 


-100.8247 


21.75 


0.4067457 


-72.109290 


00:24:24.284 


-72:06:33.44 


26.90 




95.67 


M000002 


1285.47 


984.51 


4031.80 


4172.30 


84.9804 


-100.7607 


23.83 


0.4066977 


-72.109273 


00:24:24.111 


-72:06:33.38 


29.44 




86.33 


M000003 


1301.39 


985.79 


4015.78 


4173.04 


87.6399 


-100.6749 


19.11 


0.4068580 


-72.109249 


00:24:24.688 


-72:06:33.29 


27.79 




99.56 


M000004 


1248.20 


987.50 


4068.54 


4164.60 


86.0665 


-100.6557 


23.00 


0.4067632 


-72.109243 


00:24:24.347 


-72:06:33.27 


25.31 




92.87 


M000005 


1279.67 


987.89 


4037.15 


4168.19 


86.9218 


-100.6459 


23.29 


0.4068147 


-72.109241 


00:24:24.533 


-72:06:33.26 


26.41 




93.39 


M000006 


1262.56 


988.08 


4054.17 


4165.84 


82.3194 


-100.3598 


19.76 


0.4065373 


-72.109161 


00:24:23.534 


-72:06:32.98 


27.28 




98.35 


M000007 


1354.61 


993.80 


3961.76 


4171.78 


86.2647 


-100.3566 


19.76 


0.4067751 


-72.109160 


00:24:24.390 


-72:06:32.97 


24.29 




98.95 


M000008 


1275.71 


993.87 


4040.34 


4161.74 


87.9017 


-100.2980 


23.65 


0.4068738 


-72.109144 


00:24:24.745 


-72:06:32.91 


25.16 




92.68 


M000009 


1242.97 


995.04 


4072.81 


4156.43 


81.5615 


-100.2525 


19.87 


0.4064916 


-72.109131 


00:24:23.369 


-72:06:32.87 


24.61 




97.97 


MOOOOlO 


1369.77 


995.95 


3946.40 


4171.56 



Note. — Tablel4lis published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here to illustrate its form and content. 
Key to columns: 

Column (1) — RA offset from the cluster center in Table|^ measured in arcsec at date 2002.26, from ACS images taken as part of program GO-9028. 
Column (2) — Dec offset from the cluster center in Table|3l measured in arcsec at date 2002.26, from ACS images taken as part of program GO-9028. 
Column (3)— Calibrated F475W magnitude. 

Column (4) — Calibrated absolute right ascension (J2000) in fractional hours. 

Column (5) — Calibrated absolute declination (J2000) in fractional degrees. 

Column (6) — Calibrated absolute right ascension (J2000) in (hour: minute: second) format. 

Column (7) — Calibrated absolute declination (J2000) in (degree:arcmin:arcsec) format. 

Column (8) — Faintest magnitude a star could have at the given position to be counted as "OK.'" See text for details. 

Column (9) — OK flag. Either 0K=1 (for a bona fide stellar detection) or OK=0 (for a likely PSF artifact). Set by the discrimination criterion illustrated in Fig. [71 

Column (10) — Percentage completeness: the fraction of the local neighborhood (within 100 WFC pixels, or ::± 5") where the star could have fallen and still be found. Each entry in this table thus 
con'esponds to a completeness-con'ected total "number" of (100/r) stars at the given position and magnitude in the master frame. 
Column (11)— Sequential ID number, i-unning from MOOOOOl to M129733. 

Column (12) — Star's x position, in pixels, in our meta-image of the master frame (4 7TucMaster .fits, available in the electronic edition of the Astrophysical Journal 
Column (13) — Star's v position, in pixels, in our meta-image of the master frame (4 7TucMaster . fits, available in the electronic edition of the Astrophysical Journal 
Column (14) — Star's a" position, in pixels, in the original ACS/GO-9028 image j8cd01a9q_w. fits. 
Column (15) — Star's v position, in pixels, in the original ACS/GO-9028 image j 8cd01a9q_w .fits. 



TABLES 

Displacement AND Proper-Motion Data for 14,366 Stars in the Core of 47 Tucanae 



DATASET 


MEYLANel 


MEYLANe2 


GILLILUl 


MEYLANe3 


GILLILU2 


WFC-MEUR 


HRC-MEUR 


HRC-BOHL 


WFC-KING 


HRC-KING 


DATE 


1995.82 


1997.84 


1999.51 


1999.82 


2001.53 


2002.26 


2002.26 


2002.28 


2002.52 


2002.56 


M000872 


8.3886 


-88.5393 


18.317 


19.165 


18.44 


0.4020806 


-72.105878 








M000872 


7.18E-03 


n/a 


1.21E-03 


3.08E-03 


n/a 


6.96E-05 


n/a 


n/a 


-6.46E-05 


n/a 


M000872 


3.11E-03 


n/a 


2.12E-03 


1.81E-03 


n/a 


3.48E-04 


n/a 


n/a 


2.98E-04 


n/a 


M000872 


9.46E-03 


n/a 


l.llE-03 


1.78E-03 


n/a 


-4.97E-06 


n/a 


n/a 


-1.69E-04 


n/a 


M000872 


2.91E-03 


n/a 


1.63E-03 


2.46E-03 


n/a 


3.23E-04 


n/a 


n/a 


1.44E-04 


n/a 


M000872 


-9.50E-04 


3.43E-04 


0.264 


0.8515 


-9.89E-04 


3.13E-04 


0.764 


0.5141 






M000979 


4.5367 


-87.9300 


18.649 


19.517 


18.61 


0.4018484 


-72.105709 








M000979 


5.68E-03 


n/a 


-6.36E-04 


-2.59E-04 


-1.77E-03 


3.98E-05 


n/a 


n/a 


2.14E-04 


n/a 


M000979 


3.55E-03 


n/a 


1.84E-03 


3.71E-03 


2.33E-03 


1.99E-04 


n/a 


n/a 


2.78E-04 


n/a 


M000979 


3.48E-04 


n/a 


2.08E-03 


1.89E-04 


2.84E-03 


1.19E-04 


n/a 


n/a 


-1.59E-04 


n/a 


M000979 


3.42E-03 


n/a 


1.53E-03 


2.58E-03 


2.44E-03 


2.04E-04 


n/a 


n/a 


3.03E-04 


n/a 


M000979 


-2.27E-04 


3.82E-04 


0.799 


0.5254 


-4.19E-04 


3.43E-04 


0.544 


0.7034 







Note. — Tablelslis published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here to illustrate its form and content. The first two lines of the Table 
contain the names and dates of the 10 HST data sets we have used for position measurements (see Table Izl. Following are 6 lines each for a total of 14,366 stars which have calibrated V 
magnitudes brighter than V < 20.5, colors < ((7 - V) < 3, and proper motions obtained as described in j|3lof the text. The lines recorded for each star contain the following information: 
Line 1: Columns are: (1) the star's ID (coixesponding to those in Table|4l: (2) its right-ascension distance from the cluster center (in arcseconds, positive Eastward, measured on the master frame); 

(3) its distance in declination from the cluster center (in arcseconds, positive Northward, measured on the master frame); (4) its V magnitude; (5) its U magnitude; (6) its F475W magnitude; (7) 
its absolute RA position (J2000) at date 2002.26 (program GO-9028); and (8) its absolute Dec position (J2000) at date 2002.26 (program GO-9028). 

Line 2: Columns are: (1) the star's ID, followed by (2)-(l 1) its RA displacement (in arcseconds), relative to the position in Line 1, in each of the 10 HST data sets or "epochs" we have analyzed 

(following the order in the top two lines of this table). An entry of 'n/a' in any column indicates that the star's position in that data set was not used in the proper-motion measurement. 

Line 3: Columns are: (1) the star's ID, followed by (2)-(ll) the uncertainties (in arcesonds) in the conesponding RA displacements in Line 2. 

Line 4: Columns are (1) the star's ID, followed by (2)-(l 1) its Dec displacement (in arcseconds), relative to the position in Line 1, in each of the 10 HST data sets. 

Line 5: Columns are: (1) the star's ID, followed by (2)-(l 1) the uncertainties (in arcesonds) in the conesponding Dec displacements in Line 4. 

Line 6: Columns are: (1) the star's ID, followed by (2) its weighted least-squares proper-motion //q- in the RA direction (in units of arcsec yr~' , positive Eastward); (3) the uncertainty in ', 

(4) the of the weighted least-squares straight-line fit yielding fi^', (5) the probability P that this or a higher x'a could occur by chance if the true RA-motion of the star follows a straight 
line; (6) the star's weighted least-squares proper-motion fis in the Dec direction (in units of arcsec yr~' , positive Northward); (7) the uncertainty in fis \ (8) the x" of the weighted least-squares 
straight-hne fit yielding \ and (9) the probability P that this or a higher x's could occur by chance if the true Dec-motion of the star follows a straight line. 
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TABLE 10 

Proper-Motion Kinematics in the Core of 47 Tucanae 



W) = (Me) sO 



Aimulus 
[arcsec] 
(1) 


Retf A/" (jUa) (M5> 

[arcsec] [masyr"'] [masyr"'] 
(2) (3) (4) (5) 


[mas yr~'] 
(6) 


[mas yr~' ] [mas yr~' ] 
(7) (8) 


[mas yr" 
(9) 


[masyr '] 

(10) (11) 






14 < y < 16.75, 1.5 < (U-V) < 2.6, 


and 


16.75 <V < 18.5,0.6 < {U- 


-V) < IS 


);AR = 3" 



0.0-3.0 
3.0-6.0 
6.0-9.0 
9.0-12.0 
12.0-15.0 
15.0-18.0 
18.0-21.0 
21.0-24.0 
24.0-27.0 
27.0-30.0 
30.0-33.0 
33.0-36.0 
36.0-39.0 
39.0^2.0 
42.0^5.0 
45.0^8.0 
48.0-51.0 
51.0-54.0 
54.0-57.0 
57.0-60.0 
60.0-63.0 
63.0-66.0 
66.0-69.0 
69.0-72.0 
72.0-75.0 
75.0-78.0 
78.0-81.0 
81.0-84.0 
84.0-87.0 
87.0-90.0 
90.0-93.0 
93.0-96.0 
96.0-99.0 



2.12 
4.74 
7.65 
10.61 
13.58 
16.57 
19.56 
22.55 
25.54 
28.54 
31.54 
34.53 
37.53 
40.53 
43.53 
46.52 
49.52 
52.52 
55.52 
58.52 
61.52 
64.52 
67.52 
70.52 
73.52 
76.51 
79.51 
82.51 
85.51 
88.51 
91.51 
94.51 
97.51 



69 

176 
303 
382 
353 
245 



-0.087!0;0«^ 

I +0.044 
-0.043 

o.oo3+!'!!3^ 



-0.081+ 



■0.005^; 
■0.030^11 

033+0-039 

143 -0.025:^P' 
50 -0.003!^-^| 
47 0.008-*^-"" 



84 ^.022+ 
119 -0.103!^-^| 

154 -o.m^^l 

■0.006+ooS 



162 
180 
209 
229 
259 



-0 075+0 046 

-0 020+0'0'*2 
"■""^"-0.048 
-0 046+0 042 



■0 032+0 040 

244 -0.024:™0| 

278 -0.030^11 

269 -0.059^11 

277 -0.033^11 

277 -0.014^11 

260 -0.04o3|i 
222 
224 
219 
227 
212 

165 -0.022 „ 

129 ^.054^1 

93 116"«''»3 

36 -0.175^™ 

28 -0.0383-S^ 



-0.028:«;«« 

n n^s+0 043 

-0 004+0 037 
u.uu^ „ 

-0.020+*>**3* 



025+0->"i 
-0 045+0 044 
-0 074+O-034 
-0.035^; 
-0.048^;«g 
-0 064+0 034 
_n 1 0+0 053 

A 020+0.097 

-0.162;°"* 

-0 047+^-^''^ 

-O 049+O OSl 
^■"^^-0.062 

-0 022+""'' 

0.038+11 
0.004:0-040 

-0 022+0 044 
"■'^^^-0.045 
-0 027+0-038 

-0 013+O"* 
-0 062+0-041 
_n 030+0-035 

-0 046+'' o35 

041+«-035 

-0 OO5+0-036 

012+*04l 
'^•"'^-0.036 
-0 025+0-041 

036+0-046 
U.UJD_f,Q44 

-0 032+0-037 

105+*'''** 
"•""-0.041 

042+0-0"" 

"•"^ -0.044 

-0.012 
0.011 
-0.007 



0.052- 



"J- "^-0.071 
5 64+0-027 

611+O-027 

616+<' ''^'' 
"■""'-0.025 

0.601^«;i 

0.584:»;«30 

0.638:!J;!]i 
0.468:»;«|| 

0.666!»;| 
0.628ii 

631+0-029 
"■"-"-0.038 
n i:g7+0.036 

-Mii 
^97+0.033 

"■■'^ -0.039 
Q (,41+0.0.30 
"■"^'-0.035 
5 62+0-025 
"■•'" -0.029 

550+O'''" 

"■-'•'"-0.036 
n S7O+0.026 

535+<''°^* 

"■-'-'-'-0.023 

565+'''027 

"■-'"-'-0.030 

539+0'''^2 

"■-'-^ -0.026 
570+0-027 
"■-""-0.035 
n S/lR+O-024 

0.567:!';«i 
0.584:«;«30 
n 57S+0-027 

U.JtJ „„33 

540+O'''^^ 

"■-'^"-0.039 
n 570+0.039 

-»-fl59 
n S7O+0.051 

0.680!^'9ZS 
0.490 



604+0 036 
"•""^-0.051 
601+0 031 
"•""'-0.035 
n <;o<;+0.021 

62{)+''-''20 

"•"^"-0.023 
612+0-023 
"•"'^-0.024 
561+0-024 
"•^"'-0.026 
f;33+0.034 
"•"-'-'-0.042 
616+'' '05 
"•"'"-0.134 
624+'' 05l 
• -Q.Q^7 
0,557+8.834 



n 624+"-"+^ 
"■"^^-0.046 
600+0 036 
"■"""-0.039 
618+'''''-'' 

Mii 0.640+^-^^^ 

545+0:026 
n 505+0.033 



0.6383;«| 
571+0.033 

"•-"-0.039 

0.624+''-''3*' 

5oi^8:8;?i 

"•■'^ -0.040 
n 5oa+0.028 
"-^° -0.029 
(,71+0.030 
"•""-0.039 

554+0-024 

"•-'-' -Q.Q24 

577+°^^ 

+8'8S 

0.586 „„,. 
596+''-o" 

"-"^"-0.039 

569+''-0^5 

"•-'"^-0.030 

579+0.023 
605+^-^^^ 

"•""^-0.032 
571+0.032 
"•-^ ' '-0.042 

546+0-030 

"•^^"-0.034 
607+0-030 
"•""'-0.035 
543+0-027 

U.JtJ „ „3. 

551+*-''3'' 

"•-'-"-0.038 

0.487 
0.477: 
0.568 
0.482: 
0.494 



n 5O5+0.033 
U.JOJ_(,Q32 

0.625^;«| 

622+°'055 
-Q-Q53 
643+'''0*i 
"■"^-^-0.044 
656+'' 034 
"■"^"-0.034 
n 5^-7+0.035 



.Mil 



692+'' 053 
"•"^^-0.066 

o.569:«:!;i 

0.580:«;«| 

595+0-022 
"•■'^■'-0.022 
590+0-024 
"•■'^"-0.024 
601+0-030 
"•""'-0.034 
679+0-034 
"•"'^-0.035 
491+0.052 
"•^^'-0.056 
703+0 061 

0.599+ 

n 6')s+8:8f§ 

5 39+0.036 
"•-'-'='-0.039 
649+0-039 

623+0-''38 
u.uzj g Q4Q 

620+0-030 

-0.035 "^"^"-0.034 
i+!!-8?? 0.667+«-03^ 
0.539:0:021 
559+*:o26 
0.564 
0.524 
604+^:**^^ 

"•""^-0.029 

521+''-024 

A 55O+0.030 

556+« 036 
"•-'-'"-0.034 
532+fl-023 
"•■'-' -0.025 
6IO+'' 032 
"•"'"-0.032 
544+0-025 
„ „32 

562+0-028 

505+8:8^^ 

"•■'"-'-0.035 

456+0-029 

"•^-'"-0.029 

466+*'-0" 

"•^""-0.040 

0.599: 
0.438 



56i+o:o32 

"■^"'-0.033 

601+0-031 

"■""'-0.035 
5 65+0-029 
"■-'"-'-0.031 
599+0.030 

566+^'028 
0.601^; 
0.608|| 

f'-528_o:o24 
595+0-031 

614+0-032 
"■"'^-0.027 
5 61+0-029 
"■-■"'-0.030 
n 579+0.034 
"■-"^-0.035 
n '581+0-027 
^'■''^'-0.031 
526+0-027 
"■-^ -0.030 

542+0-030 

521+0-035 
-Q-Q33 
545+^-037 

667+^'0^' 
n 50£+8.867 



659"™-"^° 
n 505+8.021 
"■-"'-'-0.023 
599+8.017 
"■-'^'^-0.018 

n fti s+O-016 

-»-fiiS 
607+0-018 
"■""'-0.020 
n 57^+0.019 
"■-"■'-0.021 
634+0-023 
"■"-■ -0.024 
Q 541+0.065 
"■-'^^-0.075 
665+0 042 
"■""-'-0.051 

611+0-03' 
-0-836 
0.635!«;i 
0.600:0.025 

609+0 026 

"■""^-0.028 
593+0.055 
"■-''-'-0.026 
gi 1+0.020 
"■""-0.023 
n ftl S+0-023 

570+0 018 
-QfilS 
563+0-OM 
"■-^"^-0.021 

5o^+0.0T9 

-8-01.3 
567+0'022 

-Q-Q23 
567+0'0l9 

-Q-Q18 
559+0-019 

3{'8i§ 
^■587_o:oi9 

559+0.025 

55(:+0.020 
"■-'-'"-0.020 
595+0.021 

535+0-019 
u.jjj „„22 

5 52+0-021 
"■-^-^ -0.024 

5 1 4+0-025 

"-^' -0.024 
5 02+0-026 
"-^" -0.028 

575+0-039 

0.593^: 
0.483«:i 



_0 21+0.22 
"■^'-0.19 

o-iolJis 

13+0^12 

15+0-10 

• ;8:1? 
"•"-0.12 

-0 20+0-14 
"■^"-0.14 

-o-3o:0:l^ 

0.35:0-25 

-0 23+''-30 
■ -031 
0.08™'| 
05+^:23 
"•"-'-0.21 
0.39:«'l' 

-0.27+0^'^ 

-0 21+^:'^ 
"■^'-0.18 

-0 06+0-16 
"■""-0.16 
-0 33+0-15 

021:013 
002+8: B 

0.13+^''^ 

30+^:'^ 
-0-13 
-0 27+^- '2 

9^8:1? 

-0.13 
1+0.16 
-0.16 
,+0.15 
-0.16 
,+0.14 
'-0.16 
-1-16 

.11? 



0.26:| 

o.i9:_^ 
02+0:16 

"•"^-0.16 
17+0-15 
"•"-0.16 

-o.io:o-|4 



"■-"-0.16 

-0.07:0- 16 
-0.07 

0.06 

0.35 

0.69:| 
-0.04 

0.36: 



^8:1^ 



'-0.40 



0.0-5.0 


3.54 


53 


2.5-7.5 


5.59 


119 


5.0-10.0 


7.91 


182 


7.5-12.5 


10.31 


209 


10.0-15.0 


12.75 


213 


12.5-17.5 


15.21 


176 


15.0-20.0 


17.68 


129 


17.5-22.5 


20.16 


83 


20.0-25.0 


22.64 


35 


22.5-27.5 


25.12 


20 


25.0-30.0 


27.61 


32 


27.5-32.5 


30.10 


42 


30.0-35.0 


32.60 


52 


32.5-37.5 


35.09 


75 


35.0-40.0 


37.58 


103 


37.5^2.5 


40.08 


125 


40.0-45.0 


42.57 


110 


42.5^7.5 


45.07 


124 


45.0-50.0 


47.57 


158 


47.5-52.5 


50.06 


145 


50.0-55.0 


52.56 


147 


52.5-57.5 


55.06 


151 


55.0-60.0 


57.55 


164 


57.5-62.5 


60.05 


173 


60.0-65.0 


62.55 


161 



015+'' 086 
-0.034^«'0g -0.126+«'0^i 

-0 059+^:052 

0.0285g 

"•"""-0.043 

-0.094+0-051 
— n 1 Q9+8:8^^ 

-0 1 23+0-106 
"■'^■^-0.116 

-0.30i:0;li« 

-0.185:0;103 

-0 044+0-100 
"■"^^-0.103 
204+0 086 

-O.lQl 
097+0 076 

-0.077 

_Q 110+^-06^ 

A 066+^:^^ 

017+''^'"0 
-0^020|i 

-0.038:0;057 

-0 003+0-051 
"■""-'-0.058 
-0 01 1+0-052 

-0 071+0-0'i'' 

"■""-0.055 

0.021+0-054 
o n75+8:8i? 



18.5 < y < 19, 0.6 < (U-V) < 1.9 ; AR 
0.638: 



:5" 



O+0.048 
°-0.Q70 
636+0-039 
"•"■^"-0.044 

604+0-028 

"•""^-0.035 

0.587+0-0 

0^598:«|| 
594+0-026 
-8-029 
O 5S4+0-031 
"-^° -Q.Q39 
599+0.044 

650+0-052 
"•"-'"-0.082 
A 505+0.053 
U.5S5_o 

0.548:0:0« 
0.537:0:0^7 

A ,1 / 1 +0.048 



0.687^;048 



■0.168:0.086 



A 132+0 "43 

^■'^^-0.045 
_0 137+0.043 

^.089+11 
^.018:0-046 

-0.048:0-054 

-0.157:; I 

-0.180:»;(g 

-0 130+0 143 
"■'^"-0.141 

0.038:0; |07 



666+0-043 

"■"DO_0,053 

632+0-040 

622+0-031 
"■"^^-0.039 
639+0-027 
"■"^^-0.032 
632+0-028 
"■"-' -0.033 

0.581+0-028 



0.548 
0.57O 
0.614 
0.605+ 



+8:818 



-0.157 

0.571+0-05? 



0.075 



0.126t?g^ 0.595: 
099+0 075 

A 055+0-064 
^■"^-'-0.066 
A 047+0-072 
-0-076 
015+0 063 
-0-06* 
032+0-069 

O23+0-056 
"•"^^-0-057 
-0 013+0-049 
"■"'-'-0.053 
001+«057 
"•""'-0.062 
037+0-051 
-0-051 
0.07i:0-0« 

050+'':^^'^ 
"•"-'"-Q.046 
025+0-038 
"•"^-'-0.050 
A 001+0.056 
^•""'-0.050 



0.640:0'0| 
0.620:o:o§i 

0.598:^;O| 
640+0-037 
-0-04^ 
632+0-0'*'' 
-0-065 
62 1 +0-045 
"■""^'-0.053 
A 59(:+0.031 
"■-'^"-0.041 
A 5(,4+O.036 
"■-'"^-0.043 
545+0.036 
"■^^^-0-040 
0.583: 

O.6I2: 

0.533: 

0.557+ 



m 
m 



m 
m 

'-0.04?. 



"•-'•"-0.075 
441+0.048 
"•^^'-0.067 

0.482:«^o« 
0.679-' 

0.648: , 
0.618:; 

A 57-1+0.043 

"•■"-'-0.051 
564+0.034 
"•^"^-0.042 
625+'' 037 
"•"^^-0.042 
544+0.035 

-0-043 
O 59 7+0.034 
"•-^^ -0.040 
542+0-040 
"•-^^ -0.047 
5 5 5+0-039 
u.jjj „ „4, 

n 5fis+8.833 



-U.032 

0.589:0.036 

0.608^ 
0.648:; 

631+0-025 
"■"^'-0.030 
5 66+0-025 
"■-'""-0.025 
5 1 5+0-032 

U.JIJ n „34 

590+0-046 
"■-'^"-0.049 

7O4+0-066 
"■'"^-0.082 
651+0 097 
"■"^'-0.117 
A 536+O-098 
"■-'-'"-0.107 
A 512+0-077 
"■•"^-0.096 

A 605+<'-070 

"■""-'-0.075 

0.633^;O|^ 

0.661^^o^« 
640+0-042 
"■"^"-0-05? 
634+0-056 
"■"^^-0-055 
635+0-041 
"■"-'-'-0.052 
(,] 1+0.036 
"■""-0.039 

0.582:0.040 
0.506:0;029 
5 75+0-037 

"■-^'-^-0.039 
614+0.042 

577+8:8t^ 
574+fi:8t^ 

"•-' '^-0.037 



0.624:0.062 

0.684:O;039 

631+*'-o3^ 
-0-034 
590+*^029 
"•-'^"-0.031 
604+0-030 

"•""^-0.033 

605+0 030 
"•""-'-0.030 
0(,i8+!i.«« 

0.606: 
0.571+ 



0.589+ 



0+0.084 
-0.101 

590+O.062 

"•^^"-0.076 

615+0070 

"•"^^-0.085 
A 577+0.061 
"•-'^^-0.083 

0467+0.043 

621+^:'^^^ 

646:0.833 

610+^-^^^ 
"•"'"-0.059 
553+0.044 
"•■'-'^-0.055 
544+0-035 
"•-'^ -0.034 
603+0-038 
"•""^-0.041 
A 577+0.041 
U.J / ' n (jig 

O 535+0.04o 
u.jjj p j,43 

n 5ao+0.040 



0.547 



-0.047 



656+0 038 

"•"-'"-0.034 

A £00+0.028 

-0 029 
620+0-024 

-0-025 
620+0-020 

-0-021 
A 617+o.OM 

"■"' '-0.024 

''■^^^-o-oi 

5 69+0-025 
"-^" -0-027 

A 5AO+0.033 

-0-034 
A 641+8-841 

"■"^'-0.048 
6 20+0 061 
"■"^"-0-071 
A 564+0-048 
"■-'"^-0-053 
A 566+0-048 
"■-'""-0.054 
A 565+0-049 

556-10-041 
"■-'-'"-0.042 

641+0.835 
■ -0-034 
643+O-028 

622+0-041 

A 595+0-031 
"■-''-'-0-040 
A 579+0-024 
"■-"'-0.026 
A 593+0.029 
"■-''-'-0.028 
A 543+0.028 

5 56+0-026 

5 77+0-029 

-0-032 
A 545+0-030 

561+^-030 

"•^"'-0031 



19+0-25 
"•'^-0.26 

-0 30+0- '6 
-0)5 

-0 07+0-1'' 

"■"'-Q-I5 

19+0-14 

-0-13 
O9+0 11 
"•"^-0-13 

-0 13+0-12 
■ -0-13 
-0 36+0- '7 



-0.05 
0.41 



;8:i: 
;8:i 



0.20+ 



r,+0.47 
"-0-49 
-0 19+0 50 
"■'^-0.49 
-0 37+0-44 
"■■"-0.45 
30+0-37 
"•-'"-0-33 
59+0.26 
"•^^-0.25 
1 3+0.20 
"•'-'-0.22 
-0 02+0-21 
"•"^-0.22 

0.07+0-| 
0.27+o:| 



0.23+ 



i+o:i8 

-0.18 
A 07+0- 19 
"■"'-0.18 

A 26+0.28 
J4+8.M 

26+^:^^ 
"•^"-0-20 
0.24+0-20 



0.10- 
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D. E. McLaughlin et al. 
TABLE 10 — Continued 



Aimulus 




..•V 


[arcsec] 


[arcsec] 




(1) 


(2) 


(3) 


62.5-67.5 


65.05 


162 


65.0-70.0 


67.55 


159 


67.5-72.5 


70.04 


162 


70.0-75.0 


72.54 


157 


72.5-77.5 


75.04 


153 


75.0-80.0 


77.54 


143 


77.5-82.5 


80.04 


138 


80.0-85.0 


82.54 


145 


82.5-87.5 


85.04 


115 


85.0-90.0 


87.54 


95 


87.5-92.5 


90.03 


85 


90.0-95.0 


92.53 


66 


92.5-97.5 


95.03 


46 


95.0-100.0 


97.53 


23 



[ma.syr"'] [masyr"'] [masyr"'] [masyr"'] 
(4) (5) (6) (7) 



[mas yr"' ] 
(8) 



[mas yr"'] 
(9) 



[mas yr ' ] 
(10) 



(11) 



-0 042+''-'''** 
006+''*''* 

"•""-0.052 

-0 063+'"'54 
0.040+''| 



n 544+0.032 
-Q.Q29 
529^-°^^ 

n sofi+0.025 

^70+0.025 
"■•^'^-0.026 

"■■^■^ -0.027 
Q 5<;7+0.029 

0.543tO;«g 
0.508:0;«30 

550+'' ®' 
520+'' ''^^' 
0-475!»;g2 

0.486^1" 
0.5363U^ 



0.058;*' 



029+"°^'' 
137+''''*' 

-^.004^11 

_f) (,5^+0.048 

-0 077+""**' 

-0 035+'' ''''' 

135+*'*'* 

063+''-''''' 
"•""^-0.072 



0.113: 

23 -0.022 



-0.013+ 

-0 094+' 
"■'"^-0.064 

"■"^"-0.058 
-0.078^1 

A 14C4^.I17 



569+*'"^'* 
0.513+ 
0.556 



0.538+ 



i+0.027 
-0.037 
Q c 54+0.033 
"■■^■^ -0.045 
0.545:!;;«|« 

0.490:!];!]|? 

Q4g4+0.042 

Q 47Q+8:8t^ 
559+S!:8tA 
570+^-^^^ 

n 4fin+'''''*'^ 

-Q.074 
619'™ "* 



520+" °^' 

'^•°^^-0.041 
A £10+0.037 
"•°^^-0.040 
551+0.028 
"•■^■^ -0.042 
5f,f;+0.043 

595+0.046 



0.546: 
0.540: 
0.524 
0.459 
0.493 



-'-0.059 
+0.042 

'rim 



i^o 

45o+S:o§i 



0.545+ 



0.508+ 
0.514; 
0.600 



0.588+ 



j+0.038 
-0.038 
n coc+0.038 
"■="^'-0.041 

"■""•'-0.044 
567+O.042 
"•-■"'-0.043 
Q 5on+''"32 



0.6 12^"- 
0.589^1 
0.482!^' 
0.427^- 



l?3 

"•""-'-0.277 



mi 

-0.044 
1+0.038 



n C7O+0.047 

0.543+ 
0.571 

"•569Lo.o43 

5|f;+0.031 
"•^^"-0.033 
Q 504+0037 
"•^"^-0.037 
519+0.043 
"•-"^ =^-0.048 
5,5+0.044 

503+''-''^' 
"•■'"^-0.041 
n AS 1 +0.032 

4^9+0.036 
"•^^ -0.044 

0.468!«:»i 
0.4563?!? 



-0 25^-^° 
-Q.19 

-0 ii+<'-i^ 

"•'"-0.19 

06+<' 2'' 
"•""-0.20 
25+0.19 
"•^•'-0.19 
36+0.21 
"•^"-0.20 
n 1 O+0.22 
"•^^-0.20 
-0.06+0-20 

('•Ol^oi 
0.47; 

0.57; 

0.06; 
-0.46; 



3f 



0.55^1 



19 < y < 19.5, 0.6 < (U-V) < 1.9 ; AR= 10" 



0.0-10.0 


7.07 


184 


5.0-15.0 


11.18 


353 


10.0-20.0 


15.81 


337 


15.0-25.0 


20.62 


156 


20.0-30.0 


25.50 


74 


25.0-35.0 


30.41 


109 


30.0-40.0 


35.36 


166 


35.0-45.0 


40.31 


219 


40.0-50.0 


45.28 


276 


45.0-55.0 


50.25 


303 


50.0-60.0 


55.23 


320 


55.0-65.0 


60.21 


365 


60.0-70.0 


65.19 


383 


65.0-75.0 


70.18 


358 


70.0-80.0 


75.17 


324 


75.0-85.0 


80.16 


315 


80.0-90.0 


85.15 


265 


85.0-95.0 


90.14 


175 


90.0-100.0 


95.13 


91 


95.0-10.'^.0 


100.12 


24 



n 79S-tO.033 

0.666li 
0.668:0;0i2 

706+0 038 
"•'""-0.047 

0.617:0;0fi 

613+0 043 
"•"'-'-0.052 

653+O-043 

n fi9S+0 037 

5 7 5+0.030 
"-■"■'-0.035 
n cco+0.035 

5 52+0.035 
"■•'■' -0.037 
n sco+0.033 

543+0.029 
"■■'^-'-0.036 
522+0 031 
"■■'-0.037 

513+O.ra 

"--^'-'-0.035 
524+0.040 
"--^ -0.049 
600+0 054 
"-"""-0.061 

676+''-06'' 

"-"'"-0,081 

0.710!;; I'l 



722+0 041 

690+^-^^ 
"■"^"-0.030 
674+fl 027 

"•"'^-0.030 
703+0.041 
"•'"•'-0.042 
541+0.058 
"•^^'-0.064 
561+0 054 
"•-'"'-0.063 

0.625!0:0« 
0.608!0-045 

n cir9+0.040 
"•-'•' -0.043 
527+0-039 
"•-^^'-0.038 
556+0.035 
"•-'-'"-0.040 
576+<>-030 
"•^'"-0.032 
n cao+0.029 

497+0.032 
"•^"-0.033 

520+^-030 

"•-'^"-0.033 
A +0.036 



0.053^;|^ 
-0.018ii 
-0 088+^-042 
_0 120+''-'"2 

_0 120+*'-''^2 
-Q.QS2 
-u.uys „„7j 

040+^-060 

"•"^"-0.060 

062+<'-05* 

"•""^-0.056 
017+fl-047 
"•""-0.046 
034+0-044 
"•"^^-0.041 
052+0-038 
"•"■'^-0.044 
026+0-037 
"•"^"-0.041 
046+0-039 
"•"^"-0.038 

0.080^0-039 
0.010?-®* 



0.033 



-0.193 



-o.o.'^y- 



=■-'-0.088 
■ro+O 177 



_0 201+0-051 

"■^"'-0.056 
-0 1 1 3+O-037 

"■"-'-0.037 
_0 051+0.035 

-0 042+0 076 

^■"^^-0.065 
_0 05 1+0 055 
^■"-"-0.055 

055+O 047 

O67+0-041 
"•""'-0.043 
054+0-043 
"•"■'^-0.039 
075+0-044 
"■"'^-0.041 
062+0-039 
"■""^-0.038 
011+0 038 
"■""-0.039 
1 3+0.037 
"•"'^-0.036 
A 011+0 041 
-ft!!?? 

-0 026+''-0*' 
030+0 

"•"^"-0.040 

-0.042+'i-046 



673+0 034 
647+0.025 

"•"^'-0.031 
619+0.026 
"•"'=^-0.026 
n fi3S+0.037 

A COC+0.049 
"-"'■'-0.072 

559+0.050 

"•-^-^ -0.069 

563+0 041 

"•-^"-^-0.048 

504+O 038 

"•-'"^-0.045 
A C9o+0.o40 
"•■'^ -0.038 

0.576+0-0| 
o.628!o:of 
"■enS? 

550+0 031 
"■^-■"-0.033 
5 30+0-033 
"■■'-'"-0.038 
541+0.031 
"■■'^'-0.032 

570+0-''35 

"•-""-0.040 
A m+0.040 



705+O-040 

"•'"■'-0.044 

0.630liil 
0.618^;025 

n fi/lR+0 040 
'^•"4S_oo40 
657+0-072 
"•"■"-0.082 

6 07+0 068 
547+0-049 

0.560!; 
0.608:li;i 
621+0 037 

"■"^'-0.041 
634+0-034 

A 597+0.036 
"■-■^'-0.041 
569+0-032 
"■-■"^-0.030 
A c-T.^+0.035 



-0.066 

0,1 9.r 



-0,072 



"■3fy:s:ii3S 

576+0-035 
"■^'"-0.038 
541+0-035 

4» 

'-'-""-0.040 "-'^'-0.038 
A 51 3+0.040 A 404+0.038 
-0.039 

404+O-039 4 5 8+0-043 
0.507!;i 0.599!ii;iii 

o,586:;i'«^ 0.614!;;:;!^ 



'-'-^^■-0.036 
549+0.044 
"-^^ -0.052 

564+0-058 

"-'"^-0.064 

f,i9+H.0f!3 

"•"'^-0 067 

'■"' '-II 1 12 



714+O.029 

-Q-!!?2 
661+O-018 
"•""'-0.020 

0.646:0;0|| 

0.676:0;0i? 

0.602!0;042 

0.584!O.04i 

0.587^0.03 

0.584!0.03 

A CO 1+0,030 
"•=''* i-0,031 
576+0-025 
"-""-0,027 
5 96+0-022 
"-' ^"-0.029 
A sofi+0.022 

5 54+0.022 
"-'-' -0.024 
A 1^30+0.026 
"-'-'°-Q.024 
531+0-024 
-QQ26 
541+0-027 
"-'^'-0.029 

0.518!0;032 
514+0.038 
-Q-B39 
609+0-0*5 
"-""^-0,042 

0.646+':?«5 



-0.05!0.|6 

-o-isiii 

-0 17+0-13 
"•"-0.12 

-0.16!0;1« 
0.38:0;^* 

16+0-34 
"•'"-0.37 

-0 26+0-27 

"•^"-0.25 
_0 16+0.23 

Im 

-0.20 
-,+0.19 
-0.22 

26+0-17 
"•^"-0.20 
7+0.16 
-0,18 
0.11+0,17 



0.19 

0.32+ 



0.07+ 



'-0,16 

0.29:0;lf 
0.08:0;j| 
00+0- M 

"•""-0.16 

_0 25+0-23 

"• -0.20 
-0 41+0-29 

-0.06+O-g 

U.I ^_,|^4g 



0.0-10.0 


7.07 


149 


5.0-15.0 


11.18 


237 


10.0-20.0 


15.81 


207 


15.0-25.0 


20.62 


89 


25.0-35.0 


30.41 


69 


30.0^0.0 


35.36 


105 


35.0^5.0 


40.31 


150 


40.0-50.0 


45.28 


195 


45.0-55.0 


50.25 


221 


50.0-60.0 


55.23 


243 


55.0-65.0 


60.21 


266 


60.0-70.0 


65.19 


262 


65.0-75.0 


70.18 


236 


70.0-80.0 


75.17 


221 


75.0-85.0 


80.16 


234 


80.0-90.0 


85.15 


211 


85.0-95.0 


90.14 


134 



043+0-060 
"•"^^-0.063 

036+^-0^1 

"•"-'"-0.055 
-0 050+0-059 
"•"^"-0.055 

-o^2i5i:i 

-0.208:00^* 
n n7s+0-087 

~"-""*-0,079 
040+0-066 
"•"^"-0.071 
063+0-056 
"•""^-0.056 
035+0-052 

0.035!«:«1 
0.048:0;045 
024+0-047 

-0.Q49 

0.094+«.056 

0.073; 
0.045; 

0.034: 

"•"-"^-0.060 



19.5 <V < 20, 0.6 < (U-V) < 1.9 ; A«= 10" 



670+O-040 

"•"'"-Q.Q41 

697+0-032 
703^0.038 
"•'"-'-0.044 
61 5+0-055 
"•"'-'-0.064 
3 99+0.073 
"•-'' -0,092 

0.498:0;0i 
0.468!0;0ff 
490+0-052 

"•^^"-0.063 
A 409+0.040 
'^•482_o,050 
47 I +0,046 
"•^"-0,049 

496+0-050 

"•^^"-0,052 
466+0-046 
"•^""-0.054 

406+0-051 

"•^""-0.053 

347+0-052 

"-'^'-0.053 

345+0-057 

3 29+0-074 
"•^^ -0,053 

202+0-109 

"•^"^-0.102 



mi 



022+0-063 
004+0-0** 

"•""^-0.053 

001+0-052 
"•""'-0.057 
5 5+0-082 
"•"■'■'-0.070 

-0 040+0 097 
"■"^"-0,097 

-0 091+0 076 
"■"^'-0,083 
053+0-069 
"■"■'^-0.064 
1 3f,+0.067 
"■'^"-0.062 
092+0-058 
"■"^^-0.060 

0.03 3:0-058 

0.089:0;0| 
094+0-049 

"•"^^-0.056 

0.112+0050 
036+^-^^^ 

"•"■'"-Q.Q51 
A A3O+0.050 

-mz 

009+0-''*' 
"■""^-0.053 

-0 044+0'063 

"■"^^-0.065 



0.671+O-048 
7i9+8:8^§ 

726+^-^^ 

"■'^"-0.058 

A .ro9+0.050 

"■="^^-0.068 
3 99+0-095 
"■^^^-0,157 
0.548:0;0j« 
0.463:0079 
n ^8ft+0,079 

41 9+0,063 
"■^'^-0,075 
A 377+0,063 
"■- ' '-0,072 
0.3 66:0.053 

0.418:0;056 

0.23T 
0.36O 
0.341 
0.16O 



Am 



673+0.057 

0.682:||;«i 
f; 1 7+0.080 
"■""-0.105 
3A5+0.099 
"■-'"■'-0.180 
442+0.084 
"■^^ -0.122 
479+0.088 
"■^'^-0.107 
553+0.068 
"■■'"-'-0.090 
535+0.053 
"■■'-'■'-0.067 
A CC9+0.066 
"-^- -0.077 
5 9 5+0.068 
"-^'-^-0.079 

o-506:0:0™ 

0^458j:g 

310+8:8?^ 

0.323+ 
0.246 



0704+O.059 

A j^aAM 

-Q-P53 
725^0.057 

"- '^-'-0,060 
A 599+0,085 
"■■'^^-0,104 
343+0,110 
"■-'^-'-0,183 
5 1 9+0,090 
"■■"^-0,106 
5 1 7+0.079 
"■■"'-0.096 
5 20+0-066 
"■■'^"-0.073 
4 60+0 061 
"■^""-0,073 
A 379+0,067 
"■^ ' -0,072 



-0.076 
.+O.077 



'-0.246 



o.46o:o:o68 

0.507:j|- 

0.438:°; 

307+0'"°+ 

ooo+o-"' 

"-"""-0,000 

0.08 1:"-'"^ 
246+''-"2 

"■^^"-0.246 



635+0.053 

o.64o:^:| 

0.680:»;»i 

630+0 079 
"■"^"-0,083 
449+0-104 
"■^^^-0.135 
475+0.085 
"■^'-"-0.101 
414+0.086 
"•^'^-0.109 
459+O.080 
"•^■''^-0.090 
503+0-057 
"•^"^-0.076 
552+0-058 
"-^- -0.068 
5 31+0-062 
"-^■"-0.072 
A 49 2+0.061 
"•^^^-0.072 
A 379+0.066 
"-" -0.077 

A 309+O.O65 

"-'° -0.077 
n 4RR+0.069 

ol5+8:8?^ 

"•^-'°-Q-Q87 
146+0^134 
"•'™-0.146 



0.21+O-24 

31^8:5^ 
^3? 

"•'•'-0.23 
-0 10+0-41 

"•'"-0.46 
_0 53+0.88 

"•■'-'-0.96 

0.18:0.5* 

o.43:o| 

0.25:O;*5 
-0.18:0;*o 
_0 75+0.37 

"•'■'-0.34 
-0 29+0.35 

"■^^-0.32 

0.36+'!^38 



0.32; 
-0.43; 
-2.00 
-1. 

0.96 



Mi 



